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CHAPTER I 


GENERAL INTRODUCTION 


1.1 INTRODUCTION 

Environmental Biology deals with the study of evolution 
and co-existence of Biological species in the surrounding 
environment. Generally, the evolution of the species depends 
on many factors such as over-crowding, age structure, past 
population si z e, sources of food supply, interactions with 
other species, topographical, ecological and environmental 
conditions in the habitat including seasonal and climatic 
variations. Due to rapid industrialisation and deforestation, 
the environment and ecology of the earth has been d eteriorating 
continuously leading to the extinction of rarer species belonging 
to the animal and plant kingdom and creating health hazards 
in the form of pollution to various living beings. Because | 

of various undesirable factors which arise due to ecological | 

and environmental gradients in natural habitats, the tendency | 

of the species, therefore, is to migrate from unfavourable to j 

| 

favourable regions for their survival. In general, migration j 

I 

of the species in the habitat can be studied by identifying j 

it with convective and dispersive processes. One of the ; 

important aspects in theoretical biology, therefore, is to 
investigate how convective and dispersive migration affect ! 
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the stability and persistence of model communities of intera- 
cting species* 

Keeping the above in view, the following problems have 
been studied in this thesis. 

(i) Pollutant dispersion from time dependent point and 
area sources. 

(ii) Effects of connective and (or) dispersive migration on 
the linear and nonlinear stability of systems of two 
interacting species with functional response in finite 
and infinite habitats with e cological and environmental 
gradients in and around it. 

(iii) Effects of dispersal on the stability of gonorrhea 
epidemic model. 

We may note that the available literature regarding the 
above is rather extensive wherein a large number of important 
ideas have been brought to focus. We shall, therefore, 
review only the most pertinent researches from the point of 
view of the problems taken up in this thesis. 

1.2 ENVIRONMEhm POLLUTION 

Air pollutants, such as sulphur dioxide, carbon dioxide, 
etc., are dispersed in the environment by the process of 
molecular diffusion which arises due to changes in concen- 
tration and depends upon various factors such as types and 
number of sources, stack heights, meteorological conditions 
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and the topography of the terrain. A great deal of attention 
has been devoted to study the molecular diffusion process by 
using the well known Pick's law of diffusion and these have 
been well documented by Crank [[22][j, Deiningor L 24 , 

Pasquill [[92[], Scorer [[108[], Stern []124[] and Sutton 
C 125^ • Hoffert [[56]]]] has presented an excellent review 
which deals with scientific and engineering aspect of air 
pollution by taking into account the environmental factors 
such as wind turbulence, stability of atmosphere and atmospher 
chemistry. 

Using Gaussian plume and other models, several investi- 
gators have studied the dispersion of air pollutants due to 
a point source under various conditions, Cleary et al. £ 17 [] , 
Pay and Rosenzweig [[32]], Gupta and Gupta [[46]]], Horst 
[[ 59 3 > Lamb and Seinfeld [[ 67 [] , Lubini and Tirabassi [[ 76 [] y 
Pasquill []92[], Pattle [[ 93 [] , Peters and Klin zing []94[], 
Peters and Richards [[95]]], Ragland and Dennis [] 99 [] , 

Ragland []lOu[], Smith [[119]]. Although the Gaussian plume 
model is advantageous due to its simplicity in calculations, 
it does not represent the real system in a complex terrain 
and changing meteorological conditions as it ignores the 
boundary layer effect near the ground where topography plays 
an important role, Hoffert []56[], Saffman [][ 106 [] • Pay 
and Rosenzweig []32[] have discussed a steady state two 
dimensional diffusion model suitable for predicting ambient 
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air pollutant concentrations averaged over a long period 
of time and resulting from the long distance transport of 
pollutants. In reality, the dispersion from a point source 
should he studied b;y considering a three dimensional 
diffusion equation with variable diffusivity and wind velocity, 
Lamb and Seinfeld [] 67 ]] . Analytical solutions of such 
equations for elevated point source have been obtained under 
various conditions, Ermak [] 30]] » Peters and Klinzing [1 94]]]. 
Using numerical techniques, Ragland and Dennis []99]]j have 
studied the dispersion of a pollutant from an elevated point 
source by taking topographical effects into account. They 
have considered various forms of wind speed and eddy diffu- 
sivity in determining the concentration profiles. Recently, 
Lupini and Tirabassi []76]] have proposed a modified Gaussian 
model that gives the geometric features of the solutions of a 
two dimensional steady-state advective-dif fusion equation 
under a height structured atmospheric boundary layer condition. 
Peters and Richards css: have developed a procedure to j 

t'- .-Ji 

study dispersion from point sources with simultaneous reversibly 

I 

chemical reactions. The validity of the procedure necessi- j 
tates the chemical reactions to be rapid enough relative to j 
the transport so that local chemical equilibrium is achieved | 
at receptor site* j 

i 

In recent years, the dispersion of pollutant from a 
line source has been investigated by many workers, Ghock []] 14= 3 i 


f 
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Dilley and Yen [] 27 ^ , Liu and Seinfeld C 75 H> Snarma and 
Myrup 113 ^[ , Walters L 130 H • In particular, Walters 
H 130 Z] has studied the dispersion by considering uniform 
air velocity and taking eddy-dif fusivity as a linear function 
of height. Dilley and Yen £ 27 Z| have performed a more 
general study of dispersion from a cross wind line source by 
taking power law models of horizontal and vertical velocities 
where diffusivity also varies as a power law profile. A 
comparison between this and trajectory model has also been 
made, Liu and Seinfeld [Z 7 5 ^] > Sharma and Myrup £ 113 ZH • 

Chock (Z 14 ZJ has proposed a simple line source model to 
describe the down-wind dispersion of pollutants near the 
roadway* It is potentially more accurate than the existing 
Gaussian models especially under adverse meteorological 
conditions when the plume becomes iirportant* An attempt has 
also been made by C alder £ 8 Zj "to examine the general structure 
of multiple source plume models for urban air pollution. 

Investigations relating to the study of dispersion from 
area sources have also been carried out by several workers, 
Draxler [Z 28 Z] » Egan and Mahoney [Z 29 Z] » Martin [Z 78 Z] * 

Omstedt and Rodhe [Z9lZ» Ragland [Z1 9 8 ZD > Stern £Z 124= ZD • 

Liu and Seinfeld [Z 75 Z] have compared the validity of grid 
and trajectory models of urban air pollution due to area 

Using similarity solutions, Gifford and Hanna jZ 35 3 


sources. 
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have analysed the diffusion from an area source for power 
law variations of velocity and eddy diffusivity. i. study 
of atmospheric transport of pollutants from an area source 
has been carried out by using integral metnod, lebedeff and 
Hameed [[68]. lunge and Sankar Subramanian £88] have 
studied the atmospheric dispersion of gaseous pollutants 
from a continuous area source using generalized dispersion 
theory, Sankar Subramanian and Gill £ 107 ] . Draxler £28] 
has considered a simple sector-average Gaussian dispersion 
model to develop long term average concentrations. The model 
accounts for day/night variations in mixing depth as well as 
stability changes with time, wind speed adjustment to stack 
height and a more realistic treatment of calm winds. 

Chu and Seinfeld £16] have derived a comprehensive 
mathematical model to study the dynamic behavior of urban 
aerosols by considering the effects of factors such as pollutant 
reaction, nucleation, coagulation, condensation, settling and 
washout* Barrie £6] has presented an improved method of 
modeling reversible sulphur dioxide exchange between a falling 
rain drop and air by taking the phenomena associated with 
ion-enhanced diffusion as well as the microphysical effects of 
internal circulation. More specific problems dealing with 
long range transport of airborne materials and/or the effects 
of deposition, reversible washout have also been discussed, 

Scriven and Fisher £109, 110], Slinn £117, 118], Smith £120] . 
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It may be noted here that problems such as pollutant 
dispersion from a time dependent source or effects of foggy 
environment on dispersion process have not been investigated. 
Keeping this in view, in Chapter II, the environmental 
pollution from a time dependent point source with irreversible 
chemical reaction is studied by using unsteady state three 
dimensional diffusion equation and in Chapter III, the rever- 
sible absorption of a pollutant from an area source in a 
stagnant fog layer is discussed. 

1.3 POPULATION DYNAMIC S 

Population Dynamics deals with the formulation and study 
of models governing the various types of interactions in 
biological associations. A successful model is one that 
explains what is currently happening and predicts what will 
happen in future. It can be noted that a simple model may 
not always be realistic and realism might demand more compli- 
cated models for biological and ecological systems, though 
it may be difficult to draw sufficient conclusions from them. 

The dynamical system describing the evolution of two 
interacting species can, in general, be written as 


II 

|tJ 

f^) + 

♦ (M 1 . 

V 

dN 

2 _ 
cFE~ " 

g(tfg) + 

V H 1> 

V 


(1.3.1) 
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where and Ug represent the population distributions 

of two species at time t. Uote that the expressions 
measuring the rate of growth are simply divided into two 
parts, viz., autogrowth and interaction terms. The terms 
f(M^), g (Ip) stand for the autogrowth of the species and 
<i> (IT^Ug), ) stand for the interactions between the 

two species. Since (11^, Ug) and ^(U^jNg) represent 
the interactions between the species which affect the rate 
of growth of each species in different ways, their ratio is 
a constant (different from unity). If we assume, therefore, 
that ^(E^,Eg) = p <J>(N^,Ng ), p being a constant, then (1.3.1) 
be comes 


dJSL 

dlT = + V 

dh 

W~ = g ( N 2^ + P 4> ^1» 


(1.3.2) 


The sign of f,g, <J> and (3 depend on the mode of interactions 
between the species. Depending upon various forms of , 
(1.3,2) represents different hinds of interactions. The 
classical Lotka-Vol terra equations for system of interacting 
species is a particular case of the system (1.3.2). Deaklin 
£23^1 has established restrictions on the applicability of 
the conservative Yolterra equations in ecological theory. 

Using Volterra’s and other models, Goel et al. [339]]] have 
investigated various aspects of the populations of interacting 
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species exhaustively. Vandermeer and Boucher [[129]] have 
analysed tne possible types of mutualistic interactions in 
population models. 

Experiments by Boiling [] 57 J and a number of other 
researchers, summarized in Murdoch and Oaten []85]], show 
that the type II response is the most common one and is now 
widely used in ecological modeling, Evans csin. Iisu H 603 , 
Kazarinoff and Driessche [[63J, Lin and Kahn [] 73 , 74]], 

Oaten and Murdoch [[89, 90]] . However, it should be noted 
that functional response of this form has not been considered 
in the systems of two species with mutualistic or competitive 
interactions* In this thesis, therefore, we study the 
following Boilings type of functional response 

P'KX 

V = Trisf U.3.3) 

for mutualistic and competitive types of interactions between 
two species, where a represents the strength of the functional- 
response and (3', a constant. It is noted that when a = 0, 
(1.3.3) corresponds to Volterra’s type interaction. 

As pointed out earlier, migration of the species plays 
an important role in the evolution of the species and as such 
its effects should be taken into account in making the 
mathematical models by introducing convective and dispersive 
terms. In contrast to convection, dispersal consists of the 
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apparently aimless, undirected movements of animals that 
seem to be wnolly random. There is, however, enough 
evidence in some cases that migration- to avoid crowding 
is the primary cause of dispersal, Gurney and liisbet 3 ^ H • 

An analytical treatment for the spatial diffusion of a 
biological population has been given by Skellam [U5] in 
1951, under the assumption that dispersal is due to random 
motion of individuals. In this context, a general model 
has been derived to describe the dynamics of dispersal of 
biological species whose motion can be viewed as a continuous 
process, Gurtin and Mac Camy L 48 J » Me Murtrie 3 84 3 , 

Skellam £ll6^. The effects of convective migration 
on the species evolution have also been investigated in various 
cases by making deterministic models, Gomins and Blatt £ 18 3 > 
Me Murtrie 3 84 U • 

Eventhough any realistic model of complex ecosystem 
is nonlinear and involves many species, mathematical ecolo- 
gists have largely restricted their nonlinear models to two 
or three species. In general, the system governing the 
evolution of the two interacting species with dispersive and 
convective migration can be written as (Segel and Levin 
3 112 3 3 Comins and Blatt 3 18 3 » Me Murtrie 3 84 3 ) 

^ + V r A H 1 = f(M 1 ) + «KH 1 ,K g ) + A. (D 1 aN 1 ) 

+ Vg.A H 2 = g(3Sf g ) + § + A. (1> 2 AN 2 ) 


(1.3.4) 
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defined in |_0, °°) * R, where R is a hounded domain with 
smooth hound ary | J^, Lg are the positive variable dispersal 
co-ef ficients of the species? Y^, Vg are convective velocities 
(constant vectors) of the species. 

As investigation of stability of system of the type 
(1.3.4) in its various forms is going to be an important 
aspect of our study, in the following, we give an account of 
the literature related to both the linear and nonlinear stabi- 
lity of such systems. 

(i) Linear stability 

Investigations relating to evolution, linear stability 
and continued co-existence of interacting species with 
temporal variations have been quite extensive, Brauer £ 7 [] , 
Deaklin 23 9 Rujii £ 33 ]] , Gilpin £ 37 [] , Goel et al. £ 39 [] , 

Kerner [[64]], Leou and lump son [] 69 [] , May []8l[], Maynard 
Smith [[ 83 [] , Murdoch and Oaten []85[], Murray [[ 87 [] , Post 
and Iravis [] 96 [] , Rescigno and Richardson []101[], Iravis 
and Post []128[], Vandermeer and Boucher []129]]. In 
particular, Post and Travis [[96]] have derived sufficient 
conditions for local stability of interacting communities. 
Murdoch and Oaten [[ 85 [] have presented an excellent review 
of the extensive work on factors that affect the stability 
of prey-predator interactions, but have limited their discu- 
ssion to linear stability properties. Gilpin [[ 37 [] has 
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investigated the feasibility and stability of prey-predator 
systems. May £ 81]]] has shown that essentially all models 
that have been proposed for prey-predator systems possess 
either a stable equilibrium state or a stable limit cycle. 

Such a stable limit cycle provides a satisfactory explana- 
tion for those animal communities in which populations are 
observed to oscillate in a rather reproducible periodic manner. 
Leou and Tump son £69]] have established a necessary and 
sufficient condition for stable co-existence of two species 
competing for two complementary or substitutable resources. 
Travis and Post £ 128 3 have examined the sensitivity of 
mutualistic communities to perturbations from the equilibrium 
and have investigated the response of equilibrium population 
densities to the changes in environmental and biological 
parameters. 

May |™ 82]], in his monograph, has surveyed the work 
relating to stability versus complexity in mutispecios models. 
He £ 79 £ has proposed a simple mathematical model to caution 
against the belief that the increase in population stability 
is a mathematical consequence of the increase in multispecies 
complexity. He £80£ has also discussed the relation 
between stability and complexity in ecological systems with 
many interacting species and has shown that complexity 
decreases the population stability. Steele £ 123 £ has 
considered some ecological models and has pointed put that 
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complexity can induce stability. By carefully analysing a 
class of predator-prey models, Gilpin £ 37] has concluded 
that ecological systems do persist after enrichment. 

Though Skollam [[ 115_J has considered the dispersal of 
the species as early as 1951, extensive investigations 
relating to the study of effects of migration on the linear 
stability of the equilibrium state for systems of interacting 
species are only of recent origin, Oasten and Holland £ ll[] , 
Casten and Case [[ 12 3 » Chewing £ 13 3 , Comins and Blatt [[ 18 3 j 
G opalsamy [[ 44, 45[] , Gurney and Nisbet [[ 47 [] , Hadeler et al. 

[[ 49 [] , Hadeler and Rothe [[ 50 [] , Jorne [[6l[], Kerner [[65[], 
Levin [_ 70[] , Murdoch and Oaten [[85[], Murray [[ 87 [] , Rosen 
j_ 103 [j , Segel and Jackson [[ lll[] , Segel and Levin [_112[], 

Smith [[l2l[], Steele [[l22[]. A more complete review 
regarding the development of the mathematical theory of spa- 
tially distributed populations can be found in Levin [[7l[], 
McMurtrie [[ 84 [] . 

Comins and Blatt [[18[] have discussed the dynamics of 
prey-predator association on a spatially nonuniform habitat 
by means of a biased dispersal and observed that the effect 
of asymmetry in dispersal is to stabilize the otherwise 
neutrally stable Lotka-V olterra model. They concluded that 
the stabilizing effect is not only produced by dispersal but 
may also be caused by ecological variation. Gurney and 
lisbet [[47 3 have shown that dispersal produced by the wholly 
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random motion is incapable of exerting any stabilizing 
influence, but an introduction of a suitable nonlinearity 
into the dispersal behavior of a species wnose characteristics 
are otherwise wholly linear can lead to stabilization under 
a wide range of conditions. Casten and Case cisn have 
concluded that the linear stability of the equilibrium state 
will never be upset after dispersal terms are incorporated 
into the system for each component species. 

Environmental heterogeneity and interacting species 
systems have also been discussed by various authors, G-opalsamy 
[[[45[[], Levin [[71, 72[], Maynard Smith [[ 83 [] , Murdoch and 
Oaten [[85[[J, Me Murtrie [[ 84 [j , Segel and Jackson C iii I] » 
Steele [[ 122 [} . Cornins and Blatt [[ 18[[J , Levin [[70[[], 
Murdoch and Oaten [[[ 85[[j and Smith [[ 12l[[] have all studied 
models of systems showing that the stable co-existence is more 
probable in a spatially heterogeneous environment than in a 
homogeneous environment. 

Effects of nonhomogene ous boundary conditions on the 
dynamics of dispersive populations in a linear habitat have 
also been investigated in cases of competitive and prey- 
predator systems respectively by Gopalsamy [[ 44 [] and 
Hadeler et al • [[49[[]. 

It is noted here that the linear stability of systems 
such as (1.3.4) with functional response of the form (1.3.3) 
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has not been investigated and, therefore, we study some 
problems in this thesis along these lines. 

(ii) Nonlinear stability 

For nonlinear systems, there is unfortunately no general 
procedure for obtaining solutions in a closed form. In fact, 
when a physical problem leads to nonlinear differential 
equation, one is most often content with qualitative behavior 
of the solution. It is widely recognized that analysis of 
nonlinear stability and/or domains of attraction of the 
equilibrium states is essential for a full understanding of 
the stability and persistence of ecological systems, a view 
eloquently expressed by Holling £ 58 3 . In such cases, 
Liapunov’s Direct Method is an excellent tool not only to 
study the nonlinear stability of the equilibrium state for 
systems of interacting species, but also to investigate the 
effects of variable dispersal co-efficients on the stability 
of the equilibrium state generally existing in the real 
situations. In recent years, various authors have attempted 
to derive sufficient conditions for nonlinear stability of the 
equilibrium state for systems of interacting species with 
temporal variations, Alan Hastings CO. Demetrius £ 25 3 , 
Fuji! £333> Gat to and Rinaldi £34]], Gilpin and Case £ 38]]] , 
Goh £ 40,41,423 » Goh and Agnew 43 3 > Harrison 
Holling Jl 58 !]* Hsu - C 60 I]» Rosen £ 1043 » Siljak £ 1143 > 
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lakeuchi et al. £ 126 3 . Rujii £ 33 3 > in particular , has 
considered the system with two prey and one predator and has 
shown that, for a set of parametric values, tne three species 
system can he stable, whereas a system with two prey species 
without a predator is unstable. He has also shown that his 
model can have a nonlinearly stable limit cycle in three 
species even when the equilibrium state is locally unstable. 
G-atto and Rinaldi £34]]] have considered a generalized lotka- 
Volterra model and have studied the stability properties of 
its nontrivial equilibrium state by means of an energy function 
first proposed by Volterra in the context of conservative 
ecosystems. Hsu £6C>3 has established certain criteria for 
nonlinear stability of the locally stable equilibrium state. 

Goh £423 too has investigated the nonlinear stability in a 
class of prey-predator models including the model in which the 
predator has type II functional response. Gilpin and Case 
£38 3 have demonstrated the existence of multiple domains 
of attraction in competitive communities. It should be noted 
that a system with more than one stable equilibrium state could 
feasibly be driven from one stable state to another by a large 
perturbation, Holling £ 58 3 • 


Goh £413 has considered the generalized Lotka-Volterra 

system 


air. 

aw = K i < r i + 


n 

£ 

3=1 


a i3 


N . ) , i=l,2, . . . ,n 
3 


( 1 . 3 . 5 ) 



17 


with positive definite function 

n N. 

V(N) = E c, [Xl.-ISf* - N? In (-i)] (1.3.6) 

i=l 1 1 1 1 K* 

where c 1 , Cg,,.,, c n are positive constants and 
E? (i = l,2,.».,n) is the equilibrium state of the system 
(1.3.5), He has concluded that if the nontrivial equilibrium 
state of the system (1.3.5) is feasible and if there exists 
a constant positive diagonal matrix C with c^,c , 0 ^ 

as diagonal elements such that CA+A 0 is negative definite, 
where A is the matrix formed by interaction co-efficients 
a. ., the equilibrium state is nonlinearly stable. When (1.3.5) 

X J 

is considered as a two species system, sufficient conditions 
for nonlinear stability of the equilibrium state ares 

(i) the equilibrium state is feasible 

(ii) the equilibrium state is locally asymptotically stable 

(iii) both the species sustain density dependent mortalities 
i.e. a i;L < 0, i = 1,2. 

# 

Effects of dispersal on the nonlinear stability of the 
equilibrium state for systems of interacting species have also 
been studied by various workers, Alan Hastings Z 2 H» Chow 
and Tam £ 15] , Conway and Smoller £ 19 ] » Harada and Eukao 
IZ 51 3 ? Jorne and Carmi Murray £ 86 Z , Rosen 10S Z] > 

Rothe Z 1° 5 Z » Segel and Levin £ 112 Zj , Steele £ 122 Z • In 
particular, Steele £ 122 Z has considered a prey-predator 
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system in a finite one dimensional spatial domain with zero 
flux boundary conditions and from the analysis, he has 
conjectured that diffusional effects can never damp out all 
spatial population fluctuations if the nonlinear terms are 
included. Murray Q86]3 ^ as demonstrated the error in the 
above analysis and has shown that the effect of uniform 
diffusion on the Lotka-Volterra equations for any even number 
of interacting species is to damp out all spatial variations. 
However, he has not made any conclusions about the effect of 
different dispersal abilities of the species. Jorne and Carmi 
£ 62]]] have considered unequal positive dispersal co-efficients 
in the dispersive Lotka-Volterra system for a number of intera- 
cting species with zero flux boundary conditions and have 
shown that the role of dispersal is to damp out all spatial 
variations. 

Harada and Fukao Q 51^ have investigated the co- 
existence of two competing species over a linear habitat of 
finite length by considering the nonlinear effects and have 
derived approximate conditions for the co-existence of these 
species using the fundamental mode of spatial variation. 

Alan Hastings £ 2 3 has further examined the system 
(1.3.5) with variable dispersal co-efficients i.e. he has 
considered the system 
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defined in £ 0,°°) x R, where R is a hounded domain with 
smooth boundary and ih(i=l,2, . . . ,n) are measures of the 
dispersal rate which may depend on x,t or any of the N/ s, 
with boundary conditions 

5h.(x,t ) 

= U X e 3H for t > 0 (1.5.8) 

0ft 

i 1,2, • . . ,n 

n being the unit outward normal vector to the boundary. 

Using Liapunov's Direct Method, it has been established that 
if the equilibrium state is stable without dispersal, it 
remains so with dispersal. 

It should be noted here that in the above studies the 
nonlinear stability of the dispersive system such as (1.5.7) 
has only been investigated, but no attempt has been made to 
investigate such problems with Holling’s response which might 
tend the system towards instability even if it is stable 
without response. 

In this thesis, therefore, the linear and nonlinear 
stability of the dispersive system of the type (1.3.4) with 
or without convective migration has been studied by consi- 
dering Holling's response in linear one, two dimensional and 
circular habitats using non-homogene ous and homogeneous 
boundary conditions respectively. 
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1.4 SPREAD OP GOMQRRHEA 

In general, the spread of communicable epidemic diseases 
depend upon the number of Defectives and susceptibles in the 
population, modes of communication of the disease, social, 
cultural and economic conditions of the population, environ- 
mental and geographical factors, etc. A detailed account of 
the study of epidemic diseases can be found in the monographs 
of Bailey (3 3 ,4 and in the lecture notes by Waltman C • 

An interesting review of spread of epidemic diseases has also 
been presented by Dietz j_ 26 3] • In particular, spread of 
infectious disease has been studied by various investigators 
considering temporal variations, Bailey C 3 > 4 I1. Cooke [20], 
Hethcote [3 53 , 54 31 , Waltman [3 131 3] • 

Gonorrhea is one of the communicable diseases which 
spreads by sexual contacts, Bailey no. Cooke and Yorke 
[ 3213 ]} Hethcote [3 55 3] , Lajmanovich and Yorke [3663] , 

Wichmann 3 1323] • Cooke and Yorke 3 21 3 were the first to 
develop a mathematical model for gonorrhea. Lajmanovich and 
Yorke 3*56 3] have considered a generalized model for gono- 
rrhea in a nonhomogene o us population and studied asymptotic 
properties of its equilibrium. Hethcote 3 ^ as use( I 

a host-vector model to study the control of gonorrhea. Yorke 
et ali have pointed out that saturation in a sexually 

active core population limits the incidence of gonorrhea and 



21 


this core causes gonorrhea to remain endemic. Bailey Q 5 J 
has also reviewed the existing state of development in 
mathematical modeling of this disease. 


In general, the model describing gonorrhea is of host- 
vector type involving two sets of populations i .e . males 
and females (Bailey CO. Laomanovich and Yorke 11663, 

Yorke et al . [[ 133 [] , Wichmann £ 132 ”[ ) . If S^, B^ are 
number of susceptibles , infectives of male population and 
Sg, Bg are the corresponding quantities of female population, 
then the system describing the temporal variations of the 
disease when the two populations are mixing homogeneously 
can be written as (Bailey £ 5]] ) 


dS- 


ds ? 

JL 

dT“ “ 

- qSjBg +a i K 1 

dF"" ~ ~ b 2 ^2^1 + a 2 

dB, 


dB 

dF“ “ 

l i¥s -«1 N 1 

dF“ = b 2 S 2 ¥ l - ®2 W 2 


(1.4.1) 


where a^ } ag ,b 1 , b g are positive constants. In writing 
(1.4.1), it is assumed that the incubation period (the time 
elapsing between the receipt of infection and the appearance 
of symptoms) is assumed to be negligibly small so that' an 
individual becomes infectious to other susceptibles immedi- 
ately after being infected. 

Attempts have also been made to study geographical 
spread of infectious diseases by considering dispersive 
migration of susceptibles and infectives, Bailey £3]]]* 
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Capasso £ 9 ] , Capasso and Portunato C 10 H » Max cat i and 
Pozio £ 77 "2 , Radcliffe [[ 97 J . In particular, Capasso 
and Portunato £ 10_J have studied the asymptotic behavior of 
a class of reaction diffusion equations related to models 
for the spatial spread of infectious diseases. Mar cat i and 

P ozio on have investigated the global behavior of a 
vector disease model by considering spatial spread and here- 
ditary effects which is applicable to the spread of malaria. 

They have shown that if the recovery rate is less than or 
equal to a threshold value, the disease dies out, otherwise 
the infectious people density tends to a homogeneous distri- 
bution. A theoretical expression giving the velocity of 
propagation for geographical spread of host-vector and carrier 
borne epidemics have been d eveloped by Radcliffe £ 97 3 • 

Tbieme Q 127 3 has constructed a deterministic model for the 
spatial spread of an epidemic. 

It may be noted from the above that the effects of 
dispersal on the spread of gonorrhea has not been investi- 
gated. To study such effects, spatial spread of the population 
must be considered in the mathematical model. In such a 
case, the equations (1.4.1) governing the population distri- 
bution of susceptibles and infectives can be generalized as 
(Bailey £ 3, 5 3], Capasso £ 9 "2 , Marcati and Pozio £77^], 
Radcliffe £ 97 3 ) 
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where D^, h g are dispersal co-efficients of the male, 
female populations due to self diffusion; the constants 
a^, a g are cross dispersal co-efficients of the infectives 
and a^, a g , b 1 , h g are positive constants. 

As the total number of promiscuous males and females 
are constant 


i • 6 • 

S 1 + 

Ss5 

H* 

li 

O 

and 

S 2 + 

N 2 ~ C 2 
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the system (1.4.2) can be reduced to the following form 

aw. 
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In this thesis, the local (linear) and global (nonlinear) 
stability of the system (1.4.4) is investigated to see the 
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effects of spatial dispersal on the spread of gonorrhea 
in Chapter XII. 

1.5 SUMMARY 

In the present chapter, a general introduction for 
the work presented in the thesis is given. 

In Chapter II, the unsteady state problem of pollutant 
dispersion from a time dependent point source has been inves- 
tigated by considering the effects of irreversible chemical 
reaction. In the case of instantaneous point source, it has 
been shown that the concentration decreases as time increases 
at a particular location and the point of maximum in the 
concentration-distance profile moves away from the source as 
time passes. It has also been shown that the shape of the 
plume at any instant of time is spherical and its size 
increases as the diffusion co-efficient increases or chemical 
re act ion rate decreases. However in the case of constant 
flux, it has been pointed out that the concentration at a 
particular location increases with time and it decreases as 
the distance from the source increases at a given instant. 

In this case, the concentration-distance profile does not have 
any maximum. It has also been noted that the plume shape 
is ellipsoid of revolution for large times and its size 
increases as the diffusion co-efficient increases or chemical 
reaction rate decreases. 
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In Chapter III, the reversible absorption of a pollutant 
from an area source emitting with constant flux in a stagnant 
fog layer has been discussed. It nas been noted that the 
concentration of tile pollutant in the gaseous phase as well 
as in fog droplets increase with time and decrease as the 
height from the source increases. It has been further pointed 
out that the concentration of the pollutant in the gaseous 
phase increases with the backward reaction rate, but decreases 
as the forward reaction rate increases. However, reverse 
is the case for the concentration of the pollutant in fog 
droplets . 

In Chapter IY, effects of convective and dispersive 
migration on the linear stability of the equilibrium state for 
two species system with mutualistic interactions and functional 
response have been investigated. In both finite and semi- 
infinite habitats, it has been shown that the otherwise stable 
equilibrium state without dispersal remains so with, dispersal 
under flux and reservoir conditions. It has also been shown 
that the effects of convective and dispersive migration are to 
stabilize the equilibrium state further in both one and two 
dimensional finite habitats under reservoir conditions. 

In Chapter Y, effects of convective and dispersive 
migration on the linear stability of prey-predator system 
with functional response have been studied. In both the 
cases of finite and infinite habitats, it has been shown that 
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an otherwise stable equilibrium state without dispersal remains 
stable with dispersal as well under flux and reservoir condi- 
tions. In the case of finite habitat, it has been noted that 
the instability of the equilibrium state is not affected by 
dispersal under flux boundary conditions, whereas dispersal 
can stabilise an otherwise unstable equilibrium state without 
dispersal under reservoir conditions. further, under reser- 
voir conditions, it has also been pointed out that the effect 
of convective migration of the species is not only to stabilize 
the already stable equilibrium state, but also to stabilize 
even an otherwise unstable equilibrium state provided a certain 
condition involving convective and dispersive coefficients 
holds . 

Moreover, it has been pointed out that the equilibrium 
state which is unstable in one dimensional habitat may become 
stable in two dimensional finite habitat under reservoir 
conditions. 

In Chapter VI, similar effects have been investigated 
for a competitive system. Under both flux and reservoir 
conditions, it has been shown that an otherwise stable equili- 
brium state without dispersal remains so with dispersal in the 
case of finite and infinite habitats. However, in the case 
of finite habitat, dispersal may stabilize even an otherwise 
unstable equilibrium state under reservoir conditions. Also, 
the effect of convective migration is to stabilize an otherwise 
unstable equilibrium state with dispersal. 



27 


It has been further pointed out that an unstable equili- 
brium state in one dimensional habitat may become stable in two 
dimensional finite habitat under reservoir conditions, as in 
the case of prey-predator system, provided certain inequality 
involving convective velocities and dispersal coefficients is 
satisfied . 

In Chapters VII and VIII, studies similar to those in 
Chapters V and VI have been carried out in circular finite 
habitats for prey-predator and competitive systems with 
functional response by considering the effects of dispersal* 

In Chapter IX, effects of variable dispersal co-efficients 
on the linear said nonlinear stability of the equilibrium state 
for two species system with mutual is tic interactions and 
functional response have been investigated using Liapunov's 
Direct Method. It has been found that the equilibrium state 
is nonlinearly asymptotically stable in the entire positive 
quadrant of the phase-plane with or without dispersal. 

In Chapter X, similar effects on the linear and nonlinear 
stability of the equilibrium state for prey-predator system 
with functional response have been studied. The functional 
response on the equilibrium state has been shown to have a 
destabilizing effect. It has also been asserted that an 
otherwise stable equilibrium state without dispersal remains 
stable with dispersal as well. Further, if the equilibrium 
state is linearly stable without dispersal, then there exists 
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a subregion of the positive quadrant in the phase-plane 
where it is nonline arly asymptotically stable with or without 
dispersal. 

Various cases under which the equilibrium state is 
nonlinearly asymptotically stable are discussed and the 
corresponding regions of stability are also estimated* 

By considering the effects of dispersal, the linear and 
nonlinear stability of the equilibrium state for two competing 
species system with functional response have been examined in 
Chapter XI. It has been pointed out that the functional 
response possesses a destabilizing effect on the equilibrium 
state. Further, it has been observed that an otherwise stable 
equilibrium state without dispersal remains so with dispersal. 

The regions of attraction in various cases under which 
the equilibrium state is nonlinearly asymptotically stable 
have been estimated. 

In Chapter XII, effects of dispersal on the linear and 
nonlinear stability of the endemic equilibrium state for gonorr- 
hea epidemic model have been investigated. It has been noted 
that the equilibrium state which is nonlinearly stable in the 
feasible region of the positive quadrant in the phase-plane 
without dispersal remains so with self dispersal when cross 
dispersal co-efficients of the infectives are zero. Even in 
the presence of cross dispersal co-efficients, if the equilibrium 
state is linearly stable with dispersal (self and cross), there 
exists a subregion of the feasible region where it is non- 
linearly asymptotically stable. 
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CHAPTER II 


ENVIRONMENTAL POLLUTION PROM A TIME 
DEPENDENT POINT SOURCE 

2.1 INTRODUCTION 

The dispersion of air pollutant due to point sources 
has been investigated by several workers using Gaussian plume 
and other models under various conditions, Gill and Sankar 
Subramanian £4,5^, Hoffert £ 7 , Lamb and Seinfeld £83* 

Pasquill Q9U> Rattle cion. Saffman j~ ll^j , Smith Q 12,13^]. 
In particular, Lamb and Seinfeld £8^ have studied the 
dispersion of a pollutant from a point source by considering 
steady state three dimensional diffusion equation with 
variable diffusivity and wind velocity. Analytical solutions 
of such equations for an elevated point source have been 
obtained by Ermak £3^] under certain conditions. It may, 
however, be noted here that the dispersion of pollutants from 
a point source with time dependent flux has not been studied 
so far. 


Keeping these in view, in this chapter, the pollutant 
dispersion from a time dependent point source is investigated 
by considering time dependent three dimensional diffusion 
equation with first order irreversible chemical reaction. 
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2.2 BASIC EQUATION 

Consider the dispersion of a reactive air pollutant 
from an elevated (H being the height of the source) point 
source in a stable environment so that it is reasonable to 
use the unsteady state convective diffusive equation to 
describe its concentration distribution. By taking into 
account the effects of ciiemical reaction and wind velocity, 
the equation governing the concentration of the pollutant 
is written as follows i 



where V (assumed to be a constant) is the mean wind velocity 
along the horizontal x-axis, D is the diffusion coefficient 
and k is the rate of chemical reaction. 

Taking the source as the origin of the co-ordinate 
system, the initial and boundary conditions may be written as 

c(s,0) = 0 for all s = (x^ + y^ + z^)V^ > 0 

c(s,t ) — 0 as s — 00 , for t > 0 (2.2.2) 

-4tds^D |“ = w(t) ass-0, for t > 0 . 

0 s 


The last boundary condition implies that the point source 
has a prescribed time dependent flux W(t). It may be noted 
that the equations (2.2.1) and (2,2.2) are also applicable 
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approximately in a non homogeneous environment where diffusion 
co-efficients in the x,y,z directions might slightly vary 
from each other (see Appendix), 

To solve the system (2.2.1) with conditions (2.2.2), 
taking Laplace transform of (2.2.1) and keeping in view the 
initial condition, we get 


V 


5c 

9x 


- 2 - 

D( — rr + 
9x 




(2.2.3) 


where c is the Laplace transform of c, m = k+p and p is the 
Laplace’s variable. 


The boundary conditions become 
c( s,p ) — 0 as s — 00 

(2.2.4) 

-4as 2 D = W(p ) as s -* 0 
where W(p) is the Laplace transform of W(t). 

Assuming the solution of (2.2.3) in the form 

o = exp (|fj) f(s) (8.2.5) 

and which on using (2.2.3) and (2.2.4) gives 



( 2 . 2 . 6 ) 
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By finding inverse Laplace transform of (2.2.6), the 
solution of (2.2.1) satisfying the initial and "boundary 
conditions is written as 

1 -i Y + i°° - Y X 

c(x,y , z,t ) = f _ V(P) exp Qpt + 

Y — 1°° 

- (”■“£• + §) ly/ ' 2 dp} (2.2.7) 

where 7 is a real positive number such that all the 
singularities lie in the left hand side of the line, fiji(p) = Y , 
in the Bromwich contour, fig. 2.1. 

2.3 PARTICULAR GASES 

To study the practical situations, various forms of the 
flux W(t)have been considered. 

Case_I : 

When the point source is instantaneous, W(t) may "be 
taken as 

W(t) = W 0 6 (t ) 

where 6(t) is Dirac delta function. Tnis case may be 
applicable to an instantaneous exhaust of pollutant from a 
stack or to an explosion of a nuclear device, a gas storage. 
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In this case, using the above expression of W(t) in 
(2.2.7), the concentration c^(x,y,z,t) can be written as, 

°1 = OTI C (_ iM exp { pt + !§ - (E^) 1 / 2 S) dp]] 

(2.3.1) 

where 

2 

b = lyr + k . (2.3.2) 


The integrand in (2.3.1) has branch point at p = -b. 
Using the Bromwich contour of fig. 2,.l, this integral can be 
evaluated to give the concentration distribution as 


W n exp 


c., = 


(Si) 


— . 

4tc Ds 


J 

b 


exp (-ut ) sin 


c ) 1/z 


s 31 du. 


(2.3.3) 


It may be noted from (2.3.3) that o ± tends to zero asymptotically; 
as t f 


It is also noted that tne expression (2.3.3) can be 
simplified to the following form 


W. 


(4-jrDt 




exp Clf ' + k)t " 


(2.3.4) 


When V = 0, k = 0, it reduces to the same form as given by 
Carslaw and Jaeger j“ 1 3 * 




! 
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The shape of the plume at a particular instant of time 
t = t Q (say) is obtained by putting c ^ = c o (a constant) in 
(2.3.4), giving 

2 „2 „ 

(/[jJ t k) t Q — (a constant) . (2.3.5) 

2 2 2 2 

Using s = x + y + z , (2.3.5) is written as 


(x - Vt Q ) 2 + y 2 + z = R 2 (2.3.6) 

where 

R 2 = 4Dt 0 (Q o - kt Q ). (2.3.7) 

It is easily seen that (2.3.6) represents a sphere of radius 
R in (x,y,z) co-ordinates with centre at (Vt ,0,0) moving with 
velocity Y along the horizontal x-axis. It may be noted that 
the size of the plume increases as D increases or k decreases. 


Further, the radius of the spherical plume increases 
as time passes i.e. as t Q increases till t Q < Q Q /2k and the 
rate of increase of this radius increases as D increases or k, 
t decrease . 


The concentration at the ground level can be written, 
from (2.3.4), as 


= exp ■ & (x " vt) ■ 


.2 (v 2 +H 2 


(2.3.8) 
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showing that it decreases as the height of the source increases. 


The concentration along the central line is also obtained 
from (2.3.4) by putting y = 0 and z = 0 as 


W 


ll (x,Q,0,t ) 


o 


/ /I -rsd- 


-yrvs- exp r -kt 




The expression given by (2.3.9) may be written in the 
dimensionless form, by using the transformations 


t — fi.. -h x — 2 y — XP 

_ vr2 _ tt2 ____ tt3 

k = jp- k , b = b, - f“ C 1 


( 2 . 3 . 10 ) 


and dropping the bars for convenience, as follows 


c ll(x,0,0,-t) = 6XP C-W (x-vt) 2 :. (2.3.11) 

Tne expression (2.3.11) for c il^ x> o,0,t) ^ s calculated 
numerically and plotted in fig. Z»2 for different values of 
t and k = 0.5, V = 2.0. It is seen from the graph that the 
concentration decreases as time increases at a particular 
location and the point of maximum in the concentration- 
distance profile moves away from the source as time increases. 



e 



Fig. 2.2 Variation of C with X for V= 2.0 , k = 0.5 
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Case XI : 


When the exhaust (flux) from the point source is uniform 


i.e. W(t ) = W 


1 


the concentration Cg(x,y,z,t) of the pollutant can be written, 
from (2.2.7), as 


W. 


1 y + n°° 1 


Co — exp (gjj) 2tcl ^ - rT ex PLP t 


'2 ~ 4tiDs 


Y~l c 


( £i^)V2 s3 ap . 

(2.3.12) 


In the above expression the integrand has a simple pole at 
p = 0 and branch point at p = -b. Using the Bromwich contour 
of fig. 2.1, this integral can be evaluated as 

Vx 

W 1 exp ^ r ,Y 2 . k^/2 

c 2 = — OTT {exp C“ ^ + B J 

- ~ i exp (-ut) sin f (-^) 1//2 s|] du}. (2.3.13) 

71 ^ U 

The asymptotic form of (2.3.13) may be found by taking 
the limit as t 


= - Wl 
c 2 


OT5 exp C§J - (—£ + |) 1/2 o 


(2.3.14) 


It is noted from (2.3.14) that Cg decreases as k increases. 
The expression (2.3.14) is same as the one obtained by Gupta 
and Gupta £ 6 . 
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The shape of the plume near the source may he obtained 
by taking Cg = c Q (a constant; in (2.3.14), giving 


2 

(~"~2 + ^)"^ 2 s — (a constant) 

2 2 2 2 

which on using s = x + y + z , simplifies to 


k x 2 + (Jfj + k) (y 2 + z 2 ) - Q 1 Vx = DQ 2 


(2.3.15) 


It may be seen that (2.3.15) represents an ellipsoid of 

VQ 1 

revolution in (x,y,z) system with centre (~^g ,0,0) and semi 
axes (A,B,B) where 

Q 2 9 

.2 y l 


a = x l d + seU 

P DQ? 

B = - 7 -— . 

k 


(2.3.16) 


It is observed from (2.3.16) that the size of the plume 
decreases as k increases for fixed D. Further, the size of 
the ellipsoid increases with D. 

The concentration at the ground level is obtained from 
(2.3.13) as 

W 1 exp ^2D^ t r- ,V 2 , k^l/2 

° 2 l( X; y,-H,t) - {expt( ^ +5) 

(x^+yt-H 2 ) 1 / 2 ] 

- | i exp (-ut (^) V 2 (x V +H 2 ) nd u } . 

71 


(2.3.17 ) 
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It may be noted from (2.3,17) that the concentration at the 
ground level decreases as the height of the source increases. 


Ihe concentration along tne central line can also be 
obtained from (2.3.13) by putting y = 0 and z = 0 as 



(x,0,0,t ) 


exp 
" I’ldDx" 




exp 


r r + 

L*“v — 5f + 

4IT 


S) l/2 


X 


oo 

- - f I exp (-at) sin £ (HzS.) 1 / 2 x ^j du} . 

(2.3.18) 


Using the same dimensionless transformations as in the 
previous case and defining 



(2.3.19 ) 


the expression (2.3.18) is written in the dimensionless form 
as follows (the bars have been dropped for convenience) i 


c 2 


(xj 0,0, t) 


t 


exp (2f) 
Ircx 



- | jf - exp (~ut )sin£ (u-b 


k) 1 / 2 

) 1/2 X 3 du}> 

(2.3.20) 


The expression (2.3.20) for Cg is calculated numerically 
and plotted in fig. 3>.3 for different values of t and k = 0.5, 
V = 2.0. It is observed from the graph that the concentration 
of the pollutant at a particular location increases with t 



CE 

Ace 


r a5»*®~ 


l 


1 
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and decreases as x increases at a given instant. It is also 
noted that the concentration-distance profile does not have 
any maximum in this case. 

From the above two cases, the following results regarding 
concentration distribution can also be obtained as shown below j 

Case II I : W(t) = + W 0 6(t - t Q J , t > 0 

This case corresponds to constant flux for all time t 
and an instantaneous flux at time t = t . The concentration 
c 3 (x,y,z,tj of the pollutant in this case is given by 
c 3 (x,y,z,t) = c 2 (x,y,z,t j,o < t < t Q 

= c 2 (x,y,z,t) + c 1 (x,y,z,t-t Q ) , t > t Q 

where c 1 (x,y,z,t~t ) is obtained from (2.3.4) by replacing t 
by (t-t Q ). When t = 0, the concentration of the pollutant 
is given by 

c 3 (x,y,z,t) = c 1 (x,y,z,t) + c 2 (x,y,z,t). 

Case IT : W(t) = W Q 6(t), 0 < t < t Q 

= w 1 , t > t Q 

This case corresponds to an instantaneous flux at time 
"t = 0 and constant flux after an elapse of time t Q . In this 
v ‘^ase the concentration c 4 (x,y,z,t) is given by 
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c 4 (x,y,z,t) = c 1 (x,y,z,t) , 0 < t < t Q 

= c 2 (x,y,z,t-t o ) + c 1 (x,y,z,t) , t > t Q 

where c g (x,y, z,t-t Q ) is the expression obtained from (2.3.13) 
by replacing t with (t-t ). 

Ca se V : W(t) = W Q 6(t) , 0 < t < t Q 

= W 0 6(t-t 0 ),t > t o 

This case corresponds to an instantaneous flux at time 
t = 0 and another at t = t . In this case concentration 
c 5 (x,y,z,t) is given by 

j z jt ) = c^(x,y,z,t) , 0 K t < t^ 

= c 1 (x,y,z,t) + c 1 (x,y,z,t-t 0 ) , t > t Q . 


2.4 CON 01 US ION 

By considering the effects of irreversible chemical 
re act ion, the unsteady state problem of pollutant dispersion 
from a time dependent point source has been investigated. In 
the case of instantaneous point source, it has been shown 
that the concentration decreases as time increases at a 
particular location and the point of maximum in the concentration- 
distance profile moves away from the source as time increases. 

It has also been shown that the shape of the plume at any 
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instant of time is spherical and its size increases as the 
diffusion co-efficient increases or chemical reaction rate 
decreases. However, in the case of constant flux, it has 
been pointed out that the concentration at a particular 
location increases with time and it decreases as the distance 
from the source increases at a given instant. The concentration- 
distance profile does not have any maximum in this case. It 
has also been noted that the plume shape in this case is 
ellipsoid of revolution for large times and its size increases 
as the diffusion co-efficient increases or chemical reaction 
rate decreases. 

In both the cases, it has also been noticed that the 
concentration at the ground level decreases as the height of 
the source increases. 
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APPENDIX 


Consider the following partial differential equation, 

(i) 


2 

a c 


iL 2 4. v -d~ c - . _ -q , t) 

at + v ax 1 " ^1 + ^2 


a 2 c a 2 c 

5" + S 3 2 “ kc 

0Xg ax 3 


where D^, Dg, D^ are the diffusion coefficients in the directions 
of x^jXg and x 3 respectively. 


the flux Q at the source may be written as 


Q 


where 


s -I = 


•S.D. ||- = W(t), for s ± 


2 2 2 

x 1 x 9 x* 

+ D^ + ^ 


0 


(ii) 


(iii) 


/% /\ 

represents ellipsoid with semiaxes (f jp s^, V p- s^,s^) 

3 3 

enclosing the source (s^ -* 0), S is the surface area of the 
ellipsoid and D is the diffusion co-efficient along the normal 
to this surface. When D^, Dg, D^ are slightly different from 
each other, the surface area of the ellipsoid can be approximated 
as 

✓vr— Arr*“ 

I 2/3 


411 CV j) 


AT /» 2 _ ; 

S 1 • ^ d s l ,s lJ 


4tt;s. 


(iv) 


The diffusion co-efficient D along the normal to the 
surface having direction cosines (l,m,n) relative to the 
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principal axes of diffusion (see co. p . 7 ) is given by 


D = l 2 D x + m 2 D 2 + n 2 I> 3 

and hence in the present case 

2 

X„ 


D. 


D 


x 


= D. 


1 \2 




3 * 


(v) 


Prom (ii) to (v), we have 


4tc s 2 Dg ||— = W(t), when s-j_ - 0. 


(vi) 


Using the transformations 


D, 


x = C^) 1 / 2 x 1 , y = (^|) 1/2 x 2 , 


Z = Xr 


the equation (i) and the condition (vi ) reduce to 


3c 

aT 


a „ a 2 , 2 , ,2 

+ v I# = d( 2-£ + 2-g + - kc 


ax 


ax ay az 


(vii) 


- 47ts 2 D || = W(t), for s -* 0 


(viii ) 


P 2 2 2 

where s = x +y +z . 


Remark : It may be noted that a sphere in (x,y,z) system 
corresponds to an ellipsoid in (x^XgjX^) system. 



CHAPTER III 


REVERSIBLE ABSORPTION OF A POLLUTANT PROM AIM AREA 
SOURCE IN A STAGNANT FOG LAYER 


3.1 INTRODUCTION 


The pollutant dispersion from a time dependent point 
source has been investigated in Chapter II by considering 
first order irreversible chemical reaction. In the present 
chapter, the reversible absorption of a pollutant from an 
area source emitting with constant flux in a stagnant fog 
layer is discussed. 

Various attempts have been made to study the dispersion 
of the pollutants in the environment from area sources, Draxler 
[]] 13 > Egan and Mahoney £ 2 ^ , Lebedeff and Hameed £ 3 ]] , Liu 
and Seinfeld £ 43 * Martin L 5 ^ > Ragland Q 6 ^ , Stern [[ 10 ^ . 
Specific problems related to deposition and reversible washout 
have also been studied, Slinn Smith £ 9 ]] . 

It may be however noted that climatic effects such 
as presence of fog in the environment on the pollutant 
dispersion have not been studied. To see such effects in a 
simplified form, in the following, "t* 1 © pollutant dispersion 
from an area source in a stagnant fog layer with reversible 
absorption is studied. 
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3.2 BASIC EQUATIONS 

Consider the situation of a stagnant layer of air 
containing a uniform distribution of fog droplets into which 
a gaseous pollutant A is emitted at the ground commencing 
at t = 0. The pollutant is absorbed reversibly by the 
droplets according to 

k 

A -r— — > B 

where A denotes the species in the gas phase and B is its 
form dissolved in the droplet phase. The constants k^ and 
k^ are first order forward and backward reaction rates. The 
gas-phase concentration c(z,t) in the region 0 < z < H is 
g iven by 

B = K fv-V + Vp < 3 - 2 -n 

where c denotes the concentration of B, E is the diffusion 
P 

co-efficient of the pollutant in the gaseous phase and H is 
the thickness of the stagnant fog layer measured vertically 
from the ground surface z = 0. 

The initial and boundary conditions are written as 


c(z,Q) = 0 

for 

z > 

0 


- K I-” = Q 

■ az H 


at 

z = 0 . 

(3.2.2) 

= 0 

3z 


at 

z = H 
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As the fog is stationary, the local concentration 
c (z,t) is described by 

dc 

= k c - k c (3.2.3) 

d t a b p v 

with initial condition 


° p ( 2 , 0 ) = 0 . 


(3.2.4) 


Defining the following dimensionless variables 


cK 
QH » 


P 


c K 
QH~ 


k 

a 


tK 

H 2 

H 2 k 

a 


n - 

* — — -v T 


*b = 


E \ 

K 


(3.2.5) 


the equations (3.2.1) - (3.2.4) can be written in the 
dimensionless form as 


with 


and 


dc 

9t 




c(n,0) = 0 n > 0 

_ |1 = 1 n = 0 

on 


9c 

an 


= 0 n = 1 


(3.2.6) 


(3.2.7) 
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3 c _ _ 

jT = ~ Vp (3.2.8) 

with 

Cp ( r 'j0)=0. (3.2.9) 

'the solution of (3.2.8) subject to (3.2.9) is 
_ T 

Cp (n , t ) = k a J c(n,t») exp C T * ) 3^’ (3.2.10) 

which on substituting in (3.2.6) gives the final equation to 
govern c(n , T ) as follows: 

2- _ _ _ T 

|t = ^-§ - S a c + k a 5 b f c(n,t' ) exp GV T - t ' ) 3 dt ' • 

(3.2.11) 

On solving the equation (3.2.11) with conditions 
(3.2.7) by using Laplace and finite Fourier cosine transforms, 
we get (dropping the bars for convenience) 

c(n,r) = A 0 +B q exp (-p 1q t) + C 0 x 

oo 

+ 2 I C A n +B n exp (-P ln T ) + C n exp (- p 2n T )3 C0S ( n7tr 0 

(3.2.12) 

and which is used in (3.2.10) to give 



0.1 


0.3 


0.5 


0.7 


Fig. 3.1 Variation of C with r| for k a =1.0 , k b =10.0 




5.5 




Fig. 3.3 Variation of C with for k a =1.0 ,k b =1.0 


• • 



1.3 



Fig. 3.5 Variation of C with r^ for k a =10.0, k b = 1. 0 






5 



Fig. 3.7 Variation of C with for r=5.0 
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(i) 


(ii) 


(iii) 


(iv) 


(v) 


(Ti) 


(vii) 


For fixed t , both c and c decrease with 

P 

increasing n . 

Both c and c increase with t at a particular 

Jr 

height from the source in the stagnant layer . 

From fig. 3.1 and 3.2, it can be noted that c is 

always greater than c for k = 1, k, = 10 at a 

p a d 

particular height from the source for all times. 

For k & = k^= 1, from fig. 3.3 and 3.4, it is also 

observed that c is slightly greater than c^ for 

all values of n and t though the values of c and 

c are close. 

P 

For all values of n and r , it is noticed from 
fig. 3.5 and 3.6 that c^ is greater than c for 


At a given height, it is also noted from fig. 3.7 
that c increases with k^, but decreases as k & 
increases for fixed t . 

At a particular height from the source, it may also 
be seen from fig. 3.8 that c increases with k , 

L) <"•! 

but decreases with increasing k^ for a given instant. 


3 .4 00MCLUSI0H 


The dispersion of a pollutant from an area source 
emitting with constant flux in a stagnant fog layer has been 
discussed. It has been shown that both the concentration 
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of the pollutant in the gaseous phase and fog droplets 
increase with time and decrease as the height from the 
source increases. It has also been noted that the concen- 
tration of the pollutant in the fog droplets increases and 
the concentration of the pollutant in the gaseous phase 
decreases as the forward reaction rate increases, while 
reverse is the case when the backward reaction rate 
increases. It has been further seen that the concentration 
of the pollutant in the gaseous phase is greater or less than 
the concentration of the pollutant in the fog droplets, at 
a particular instant and height from the source, depends upon 
whether the backward reaction rate is greater or less than 
the forward reaction rate. 
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CHAPTER IV 


EFFECTS OF CONVECTIVE AMD DISPERSIVE MIGRATION Oh THE 
LINEAR STABILITY OF TWO SPECIES SYSTEM WITH 
MUTUAL ISTIC INTERACTIONS AND 
FUNCTIONAL RESPONSE 


4.1 INTRODUCTION 


In the earlier two chapters, the dispersion of 
pollutants in air causing the deterioration of the environ- 
ment and ecology, which might affect the living species in 
the habitat, has been studied. As pointed out in Section 1.1, 
such undesirable effects might induce the tendency in the 
living species to migrate from unfavourable to favourable 
regions for their survial in the natural habitats. Further, 
species can also migrate from one region to other due to over 
crowding, sources of food supply, topographical conditions 
and environmental and ecological gradients. In general, 
migration of the species can be studied by identifying it 
with convective and dispersive processes. Therefore, in 
Chapters IV to XI, the effects of convective and dispersive 
migration on the linear and nonlinear stability of systems 
of interacting species (mutualistic, prey-predator, competitive 
living in a linear or circular habitat are studied. In the 
present chapter, however, we investigate such effects on the 
mutualistic interactions of two species with functional 


response 
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The evolution and stability of the two interacting species 
systems with mutualistic interactions (symbiosis) have been 
studied by various investigators, Cunningham on, King et al. 

Q 6 n » Rescigno and Richardson [[ 113 , Travis and Post [[ 14 3 > 

V and acme e-r and Boucher £153* In particular, Travis and 
Post [[ 14 J have discussed the response of equilibrium population 
densities to cnanges in biological and environmental parameters 
and concluded that the equilibrium state is stable. It may 
however be noted here that though the effectsof functional 
response on prey-predator system have been investigated, no 
attempt has been made to study such effects in the case of 
symbiosis (Kazarinoff and Driessche C 53 , Bin and Kahn [[83 » 

May C 10 J ). The effects of convective and dispersive 
migration have also not been discussed on such interactions 
(Comins and Blatt [_l[[]» Hadeler et al. [[ 3 3 > Levin [[93s 
Segel and Levin [[ lf> 3 )• 

Keeping the above in view, in the following, the effects 
of dispersive migration on the mutualistic interactions of two 
species system with functional response have been studied in 
finite and semi-infinite habitats under flux and reservoir 
boundary conditions. The effect of convective migration on 
the stability of this system has also been discussed. 

Another aim of the chapter is to investigate the premise 
whether by making a system more complex its stability increases 
as suggested by Jacquez [[4 3 • 
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4 . 2 BASIC EQUATIONS 


Consider the mutual is tic interactions of two species in 
a finite or semi- infinite one dimensional linear habitat. By 
taking into account the effects of convective and dispersive 
migration of the species with functional response, the equations 
governing their evolution can be written as follows (Lin and 
Kahn j] 8~j , May [] 10]] , Segel and Levin [] 13 ]} , Travis and 
Post [] 14 ]J )■ 


9H-, 9N 1 p Up 9%- 

dt~ + P 1 Fx“ = a l N l" a ll N l + a 12 l+ccN 1 + D 1 7 2 

1 dx 


^2 

at 


+ v 



p nk 

a 2 N 2“ a 22 N 2 + a 21 1+aF^ D 2 



(4.2.1) 


where H^(x,t), N g (x,t) represent the population distributions 

of the two species at time tj D^,Dg their coefficients of 

dispersal and v their migration velocities. The 

interaction coefficients a q» a 2 ’ a n ,a i2 ,a 21 ,a 22 a are 

*1 

positive constants. The term (1+aN^)" determines the 
functional response in the model and the constant a may be 
interpreted as the measure of the strength of this response 
(Kazarinoff and Driessehe []5]])« 


Since the trivial equilibrium points are not of much 
interest in biological systems (Roberts [] 12 ]] ), the nontrivial 
positive equilibrium state (N*, Kg) of the system (4.2.1) 
may be obtained from, 
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a l + 


a 9 + 


fl2 !| 

1+aNj 


a 21 H l 

1+ctNj 


as 


N* 3 + ) N5 2 + 

1 v a a n 1 


: a nK 


■■ a 22 Hj . 


, we get the 

1 

O: 

2 

a ll a 22 a 


)n k*. 

clgcl. 



(4.2.2) 


Eliminating Hg in (4.2.2), we get the cubic equation for 


ll a 22~ a 12 a 21 


0. (4.2.3) 


a ll a 22 a 


After examining the nature of the coefficients of N* in 
(4.2.3), it is seen that there exists only one positive value 
for N* and therefore for N* . 

However, for a=0, solving (4.2.2) we get the following 
positive equilibrium point 


1 


a 2 a 12 + a l a 22 
a ll a 22“ a 12 a 21 


w * - a l a 21 + a 2 a ll 
J J-M o — 


a ll a 22** a 12 a 21 


(4.2.4) 


provided 


a ll a 22 > 


a 12 a 21 


(4.2.5) 


which implies that intraspecific interactions are more 
stronger than interspecific interactions. 
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To investigate tile local stability of the equilibrium 
state (H£, Hg) for a >_ 0, we linearise the system (4.2.1) by 
writing 


N^(x,t) = N* + u(x,t) 

Hg(x,t) = h* + v(x,t) . 

On substituting (4.2.6) in the system (4.2.1) 
(4.2.2), we get the following equations 


(4.2.6) 


and using 


du 
9 1 


dv 

at 


+ V 


+ V 


m . 0U + f V -+■ D 

1 9x 

2 lx = SU ~ bV + ^2 



u 

2 


(4.2.7) 


where e,f,g,h are positive constants given by, 


e = h* {3 


“ a 12 M 2 


g = 


a 2l H 2* 

( 1+aN *) 


(1+cdJf) 

2 » h = a 22 K 2 


3? f = 


a 12% 

(1+ocHj) 


(4.2.8a) 


and 

eh > fg (4.2.8b) 

for a >, 0. 

In absence of convective and dispersive migration, it 
can be easily seen from (4.2.7) that the positive equilibrium 
state (N*, Ng) is stable for a > 0. 



To investigate the stability of the equilibrium state 
for a > 0, the system (4.2.7) may be associated with the 
following conditions i 

(i) Flux conditions 

u(x,0) = r(x), v(x,0) = s(x) x > 0 

-0 1 —^4^- = U, -Dg ll (0,t) = V t > 0 (4.2.9) 

auCx, t) ^ 0 5v(x,t; 

5x * 9x 

as x - L for finite habitat or x - 00 for infinite habitat. 

(ii) Reservoir conditions 

u(x,0) = r(x), v(x,0) = s(x) x > 0 

u(0,t) = U 1 , v( 0, t) = V 1 t > 0 (4.2.10) 

u(x,t) - 0 , v(x,t) — 0 

as x - 1 for finite habitat or x - 00 for infinite habitat. 

4.3 STABILITY ANALYSIS WITH DISPERSAL IN A FINITE HABITAT 

In the following, we study the stability of the equilibrium 
state with dispersal in a finite habitat when flux or reservoir 
boundary conditions are prescribed (p-j_ = Vg = 0). 

Case (i) Flux conditions 

In this case, on solving the system (4.2.7) for 
V 1 = V 2 ~ 0 w ^- tl conditions (4.2.9) by using Laplace 



and finite Fourier cosine transforms, we obtain u(x,t) and 
v(x,t) as 

u(x,t) = [i exp(-p 10 2 ) + B o exp(-p 20 f O dT 

0 

oo "t» ... HTUX 

+ ^ S { / L A n e xp(-Pi n I) + B n exp(-p 2 n f ) 3 df } cos 3 - 

+ r exp ^ p 10 t) + F 20 exp (-£20^3 

OO 

+ I s {P ln exp(-p ln t) + F 2n exp(-p 2 n t)} cos 
n=l 

T(x,t) = -r / C C exp(-PiQ^) + E 0 exp(-p 20 T) 3 dT 
o 

00 "ij mcsm , 

+ r 2 {f 3 c n e xp(-Pin^) + B n exp(-p 2 n f) J 
n=l o 

+ % Z (} 1 0 ex 9^10 t ' ) + G 20 exp( “ p 20 t ^ 

00 JlTOC 

+ I X {G ln exp(-p ln t) + G gn exp(-p 2 n t)} cos -j- 
^ n=l 

(4.3.2) 


nxx 

“IT 


(4.3.1) 


, mi n-rcx 
df} cos 


where 
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(p-, „-b- 


In 


2 2 

% 

2 " 2 - 


21 71 D 0 ) U-fV 


^n 


h 


p ln “ p 2n 


B n = 


V h " p 2n> U+fV 

p ln ~ p 2n 


(P 


In' 


•e- D 1 ) V-gU 


0 


n 


p ln “ p 2n 


E n = 


2 2 

(—§“ + e-p 9 JV + gU 

L 


2n ' 


p In ** p 2n 


(4.3.3) 


y . 2 2 

( P ln - ti - 4 §- ft >) ^- £S . 


n n 


in" 


p = 
2n 


p ln “ p 2n 

(~j“2~“ D 2 +il “ p 2n^ R n + fS n 


p ln “ p 2n 


G ln = 


^ p ln“ e " ~f~ V VS R n 
p ln “ p 2n 


^ t 2 E l +e-p 2n^ ^n + ^ ^n 


G 2n= 

_u 

p ln * 

- ®2n 



p ln ^ 


2 2 


v, 2 2 


p 2n i 

1 

” 2 

( n k 
L 2 

D l + 

n % 
L 2 

V 


1 y,2 2 

+ IC (Hr E 


n 2 ^ 2 
1_ L 2 


(4.3.4) 


D ? +e-h) 2 + 4fg] 


1/2 


(4.3.5) 
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and ( ) , ( )iq, ( )gQ can be obtained from the 

corresponding expressions of ( ) n , ( ) ln , ( ) 2n by 

putting n = 0. The constants R , S n are given by 

R n = / r ( x ) c °s dx 
o 

S n = / s(x) cos dx, 
o 


Since p ln , p gn (n=l,2,3,...J are positive and increase 
with dispersal co- efficients (see appendix), the nonstationary 
solutions (4. 3,1) and (4.3.2) converge to the following forms 


as t — 


It u(x,t) 

t - oo 


It v(x,t) 
t 


- r 

L L. 


i 00 

0 + f2_>l z 


p 10 p 20 L n=l p ln p 2n 


r^+ Si-: cos n,tx 

U- n - n _ —l 


1 

L 


r °° + “] + 2 £ p £n_ 

Lp 10 P 20 J L n=l p ln 


E. 


n 

p 2n 




1 

(4.3.6) 

n-rcx 

L 

(4.3.7) 


and the rate of convergence increases with dispersal 
co-efficients, further, as (4.3.6) and (4.3.7) satisfy the 
steady state form of the system (4.2.7) with boundary conditions 
(4.2.9), the equilibrium state is asymptotically stable for 
a > 0 and the degree of stability (i.e. the rate of convergence 
of nonstationary solutions to stationary solutions) increases 
as dispersal co-efficients increase. 



Oase (ii) Reservoir Conditions 

On solving the system (4.2.7) for = Pg = 0 with 
conditions ( 4 * 2. 1C) by using Laplace and finite Fourier sine 
transforms, we get u(x,t) and v(x,t) as 


u(x, t) 


1 2 n ij exp(-p ln l) 

L^ n=l o 1 11 

+ exp(~p 2n l) } sin 


+ ^ n i 1 ^ I 'L exp( ' p ln t) + i ’2n exp( - p 2li t} 3 BiI1 TT 

(4.3.8 ) 

9 00 t 

V(x,t) = m 2 n {/ Cc> exp(-p ln T) 
n=l 0 n ln 

+ exp(-p 2n T)3dl} sin 


CO 

+ r n =i C s in e3 ®(-I>ln t ^ + & 2n exp (- p 2 n t >3 B “ TT 

(4.3.9) 

where • 2 g 

(Pm- -V- D lV f Vl 

t ±L . 

n Pm - p 2n 

(7# ¥^ a ) W f Vi 

n Pin ~ p 2n 

( p m“ D i” e ^ D 2 v i _g ^1^1 

c 

n p ln ~ p 2n 


B ' = 
n 


( l 2 D l +e ” p 2n /l D 2 Y 1 + S D 1 TJ 1 


p ln “ p 2n 


(4.3.10) 
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p ’ 

In 


( p ln” ^2 • D 2“ i1 ' ) R n~ f S n 

p ln “ p 2n 


p t 

2n 


V h - p 2n) R A + f S A 

p ln ” p 2n 


( p m" 7T V e) S A " s R n 

jj 


G 1 


i 


n 


p ln ” p 2n 


2n 


l 2 ' p i +e ~P 2n^ S n + s R n 
p ln “ p 2n 


(4.5. 11) 


p ln , p 2n are defined in (4.3.5) and R^, are gi^ren by 

= / L r (x) sin dx 
o 

= J s(x) sin dx . 
o 


Since P ln jP 2 n ( n=1 > 2,3 , . . . ) are positive, the ronstat ionary 
solutions (4.3.8) and (4.3.9) converge to the following forms as 
t - <=°, 


It u(x,t) 


It v(x,t) 

t->oo 



2 n 
n=l 



2 n 
n=l 




3 ' 

r_ii_ + JL.“| 

p In p 2n 


C 1 E* 

r - 2 — + -S- 3 

U n . n . — * 


p ln 


2n 


sin 


n-rtx 


T“ 


sin 


nixx 

L 


(4.3.12) 


(4.3.13) 



and the rate of convergence increases with dispersal 
co-efficients. As (4.3.12) and (4.3.13) satisfy the 
steady state form of the system (4.2.7) with boundary 
conditions (4.2.10), the equilibrium state is asymptoti- 
cally stable for a >, 0 and the degree of stability 
increases with dispersal co-efficients under reservoir 
conditions also. 


4.4 STABILITY ANALYSIS WITH DISPERSAL IN AN INFINITE 
HABITAT 


In the following, the stability analysis of the 
equilibrium state without convection (u ^ = u = 0) is 
discussed in a se mi- inf in it e habitat under flux and reservoir 
conditions (4.2.9), (4.2.10) respectively. 

Case (i) Flux Conditions 

In this case, on solving the system (4.2.7) for 
~ V 2 ~ 0 conditions (4.2.9), we obtain 


u(x,p) 


{ -' 1 7 a; e *p 


D iV ,f( P ) 


+ -ij — sxp G-y) /z (p )x 3 } 

y? (?) 

! “ expG-yi^te) 2 ! . . 

"** 2D D ^ ^ 1/p ~ > z+x ) 

2D 1 D 2 0 yV2( p) 


. exp G*yj!/^(P ) Z H 

+ F 1 (p,| z-x )3+ 1 7 §— 

y 2 7 (p ) 


”F 2 (p,z+x) + Fg(p,jz-x| )H > dz 


(4.4.1) 



v(x,p) = -j-- y exp £-y^ 2 (p )*3 

■ u l iJ 2 p y£ /<J (p) 1 

+ exp L^yi^CpjxJ > 

y 2 < p) 

i ,°° e xpG-y^ 2 (p)z] 


2D l D 2 o 


y ?%7 


L& 1 (p ,z+x) 


exp b-yg //2 (p ) C3 

+ &i(p » | z-xl ) 3+ i 7 §— 

yi 7 (p) 


_G g (p,z+x) + G g (p , 1 z-xj ) 3> dz 


(4. 


where u(x,p), v(x,p) are Laplace transforms of u(x,t) 
v(x,t) respectively and 


(D p y 1 -p-h) U-fV (p+h~D g y p )U+fV 

A(p) = - M— ; B(p) = • 82 


y l - y 8 


y l " t 


(D 1 y 1 -p-e)V-glT (p+e-L 1 y ? )Y+gU 

G (p ) = — ; B(p) = 


I’ 1 (p »x) 


F g (p,x) = 


g-lCp^x) 


y l " y 2 

(l> 2 yi-p-ii)r(x) - fs(x) 

~~2 

(p+h-H 2 y 2 ) r(x) + fs(x) 
y i-' y 2 

(Diyi-p-e) s(x)-g r(x) 


y i - y 2 


(4. 


G g (p,x) = 


(p+e-D 1 y g ) s(x) + g r(x) 


y l - y 2 


4.2J 


4.3) 


(4.4.4) 




Since eh > fg, it may be noted from (4.4.5) that y^ and y g 
are positive for all values of p >. 0 and this fact is used 
in obtaining tne solutions (4.4.1) and (4.4.2). 

To investigate the asymptotic behavior of u(x,t) and 
v(x,t) as t -» we use the following property of Laplace 
transform (see £ 7 , p. 315) 


It u(x,t) = It pu(x,p) , (4.4.6) 

t->-oo p-^0 

Using (4.4.6) in (4.4.1) and (4.4.2), we get 


It u(x,t) = 
t — *°° 


{ 


33 i d 2 r ± 


+ 


- ex P C-yi /8 ( 0 u3 


7 ^^ — exp C-y g /2 (0) (4.4.7) 


o) 


it v(x,t) = exp G-yJ / 2 (°)x^ 

t-oo 11 1 D 2 y7 /(i (o) 


+ exp C-y 2 /2 (°) xD> 

yg /2 (0) u 2 


(4.4.8) 


where y 1 (0), y g (0) are values of y 1 (p), y g (p) P=° 

respectively and A(0), B(G), 0(0), E(0) are obtained by 
putting p=0 in the corresponding expressions (4.4.3). 





It can be verified that the expressions (4,4.7) and 
(4, -4. 8) satisfy the steady state form of the system (4.2.7), 
(4,2.9). Thus it may he noted that the nonstationary 
solutions tend to the corresponding stationary solutions 
whezn t-*°=, snowing the asymptotic stability of the equilibrium 
state under flux boundary conditions. 


Case (ii) he servo ir conditions 

As before, on solving the system (4.2.7) for = P 2 = ^ 
with conditions (4.2.10), we obtain 


u(x,p ) 




+ 


1 

2D 1 D g 


{A'(p) exp [> y^ 2 (p) 3 
+ B*(p) exp [1 yg/ 2 (p) xj } 

°o expG-y^ 2 (p)z3 

{ — 


[)l' 1 (p, 2 +x) - P 1 (p,Z-x)3 


+ 


exp C-y| /8 (p )g| 
y) /8 (p) 


0 2 (p,zn-x)-i 1 2 (p,z-x) dz 


(4.4.9) 



v(x,p) = -C o * ( p ; exp G-y 1 1 // 2 (p)x J 

1 2 


+ E'(p) exp G-yJ/ 2 (p )x"2 > 


1 


2D 1 D 2 o 


/ { 


exp G-y^ //2 (p)z J 

yp^(p) 


exp C-y 9 //2 (p) z H 

[I &i(p ,z+x)-G 1 (p ,z-x) 3+ I7T7 


[G- p (p,z+x)-G ? (p,z-x)3 } dz 


(4.4, 10) 


where u(x,p), v(x,p) are Laplace transforms of u(x,t), 
v(x,t) respectively and 


A*(p) 

B‘(p) 

C f (p) 

B»(p) 


(D 2 yi-p-h) D 1 U 1 -fJ)gV 1 
~~2 “ 
(p+h^-Lgyg j D 1 U 1 +fD 2 V 1 

( D l y l“P“ e ^ ■ D 2 V r‘ s;D l lJ l 

— 

(p+e-D^g) DgV^gl)^ 

yl^yi 


(4.4.11) 


To investigate the asymptotic behavior of u(x,t) 
and v(x,t) as t — <», we use the property of Laplace 
transform given by (4.4.6) in (4.4.9) and (4.4.10) to get 
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It u(x,t) = exp L-y^ /2 (°)x J 

t-*°° l - ^ 2 

+ 3'(0) exp ^ (4.4.12) 

It v(x,t) = tAt- < C‘ (OJ exp C-yi /2 (°) €3 
t— °° 12 

+ E’(0) exp Q-Yg^CO) £J } (4.4.13) 

where A'(0), B*(0), C^u), E’(0) are obtained by putting 
p=0 in the corresponding expressions (4.4.11). 

It is seen that the expressions given by (4.4.12) and 
(4.4.13) satisfy the steady state form of (4.2.7) with boundary 
conditions (4.2.10) showing the asymptotic stability of the 
equilibrium state in this Case also. 

4.5 STABILITY ANALYSIS WITH CONVECTION AND DISPERSAL 
IN A PIN HE HABITAT UNDER RESERVOIR CONDITIONS 

To see the effect of convection on the stability of the 
equilibrium state of the system (4.2.7) under conditions 
(4.2.10), we consider the case for which 



i.e. the ratio of convective velocity to dispersal co-efficient 
is same for both the species. This problem is biologically 
reasonable because the species possessing more dispersal 
ability are likely to have faster convective migration velocity. 
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Using the following transformations 
u(x,t) = exp (^ I x/2D 1 ; u 1 (x,t) 

v(x,t) = exp (y 2 x/2Ug) v 1 (x,t; 

the system (4.2.7) can- be written as follows 
ou- 

rrr— = -e-,u- + fv^ + D-, 5?= 

0 v -i- -** -L jv 

ox 

. .2 

9v 1 d v 1 

FT = SU 1 “ h l V l + D 2 ^2“ 

where 

e 1 = e + y^/4D 1 j h 1 = h + y 2 /4Dg . 


(4.5.2) 


(4.5.3) 


(4.5.4) 


The initial and. boundary conditions (4.2.10) reduce 
to the following form 

u 1 (x,0) = r 1 (x), v 1 (x,0) = s 1 (x) x > 0 

( 0 , t ) = V , v 1 (0,t) = Y ± t > 0 (4.5.5) 

u^Xjt) — 0 , v 1 (x,t) - 0 

as x — L for finite habitat or x - 00 for infinite habitat. 
The functions r^(x) and s^(x) are defined by 

r 1 (x) = r(x) exp (- u 1 x/2D 1 ) 
s 1 (x) = s(x) exp (- ygx/2Ug) . 


(4.5.6 ) 



It is noted that the system (4.5.3) with conditions 
(4.5.5) is of the same form as of the system (4.2.7) for 


W 1 = U 2 = 0 wi'frk conditions (4.2.10) wnere e^,h^, r^(x), 
s^(x ) have been modified according to (4.5.4) and (4.5.6). 

Thus the solutions of the system (4.2.7) with (4.5.1) can be 
obtained from case (ii) of Sections 4.3 and 4.4 by making use 
of (4.5.2), (4.5.4) and (4.5.5) appropriately. 

For example, in the case of finite habitat, the solutions 
of the system (4.2.7 ) with reservoir conditions (4.2.10) for 
the case of (4.5.1) can be obtained from (4.3.8) and (4.3.9) 

U 1 X UgX 

by multiplying these equations with exp (g^~) a^d exp 
respectively and replacing e by e^, h by h^, r(x) by 
r 1 (x), s( x ) by s 1 (x) . Thus the corresponding eigenvalues 
C Pin* p 2n ( n =4> 2 »3,...)3 for the system (4.5.3) with (4.5.2) 
would be given by the expressions (4.3.5) by replacing e by 
e^ and h by h^ (see appendix). As shown in the appendix, 
it may.be noted that Pg n (n=l, 2 ,3, . . . ) are positive and 

increase with U 1 and Hence, in the case of finite 

habitat, it can be argued that the equilibrium state is 
stabilized further due to convective migration also and the 
degree of stability increases with convective velocities. 
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4.6 STABILITY ANALYSIS WITH CONVECTION AND DISPERSAL 

IN A TWO DIMENSIONAL FINITE HABITAT UNDER RESERVOIR 
CONDITIONS 


Event ho ugh migration of the species is considered only 
along a particular direction in the earlier sections,- in 
real situations, the species migrate in all possible directions. 
In such circumstances, the evolution of the species in a 
rectangular finite habitat £l,M being the length of the 
habitat along x,y-directions respectively]] is governed by 
the following dynamical equations 


5N 1 3N 1 6N 1 a 1p N 9 

F“ + y i ~ + U 3 ~ = N i( a i- a n H i + 

3 2 N, 


+ D 


) 2 N, 


+ JD 


3 ay 2 


S V U 


SN r 


p ^p i 

2 a x + y 4 Fy" = 1T 2^ a 2" a 22 N 2 + i+ocN 1 ^ 

d 2 N 0 3 2 N 

2~ + D 4 




(4.6.1) 


ax~ * ay* 

# 

where D^ Dg, D^, D^ are co-efficients of dispersal of the 
species ahd p 1 , u 3 , their convective velocities. 

As in Section 4.5, consider the case 


D-i 


D 2 ’ D 3 


D, 


(4.6.2) 


i.e. the ratio of convective velocity to dispersal co-efficient 
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is same for both tne species along the particular direction* 
Similar analysis, as in the previous articles, establishes 
that the equilibrium state is asymptotically stable in a two 
dimensional finite habitat under reservoir conditions. 

4.7 CONCLUSION 

Effects of dispersive and convective migration of the 
species, which arise due to environmental and ecological 
gradients, on the linear stability of the equilibrium state 
for two species system with mutualistic interactions have 
been investigated by taking into account the functional response 
in finite and semi-infinite habitats. In both these cases, 
it has been shown that the equilibrium state wnich is always 
stable without migration remains so with dispersal under flux 
and reservoir boundary conditions. 

In the case of a finite habitat, it has also been pointed 
out that the degree of stability increases with dispersal 
co-efficients of the species in absence of convective migration. 
It has been further shown that the effect of convective 
migration of the species is to stabilize the equilibrium state 
and the degree of stability increases as convective velocities 
increase in this case also. It has also been pointed out 
th a t the effects of convective and dispersive migration are 
to stabilize the equilibrium state further in two dimensional 
finite habitats under reservoir conditions. 
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iJPPSNDIX 

From (4.3.5), p^ n and. pg n can be rewritten as 

?r (D^+bg+e+h) 
tgZ (D^Dg+e-h) 2 + 4 fg-J V 2 (i) 

where 

„22 _ 2 2 

D 1 = - T -z- D 1> d 2 = s-r- d 2 • (“•) 

Jj JJ 

Since eh > fg, it is noted from (i) tnat P ln > P 2n 
(n = 1,2,3,...) are positive. To see the behavior of 
Pln> P 2n ( n " 1,2,3,...) w.r.t. D^,Dg and L, we find 

dp in _ 1 ^ 1 ~ + e " + 1 

d5 l 2 Q (S^Dg+e-h ,i ? V4fg j 1; ' 2 2 

Tbs = 1.1 (5 i ~ h 1 6 - 

d5 2 _5 " 2 [(vV»-« 2 + «s3 1/2 ’ 

d P 1n dp ln 

Since p. , — - — and — — — are positive for D. > 0, D 9 > 0, 

in db 1 dDg x ~ “ 

it is concluded tnat p ^ increases with D^ and Dg. Since 

D^, Dg are proportional to D-^Dg and inversely proportional 
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to L, it may be noted that p ln increases with Dp and D g and 
decreases as 1 increases. 


Similarly it can also be seen that p gn increases as 
Dp, D g increase or as 1 decreases. 


As pointed out earlier, the eigen values Pp n , p gn 
(n = 1,2,3,...) for the system (4.5.3) with (4.5.2) are defined 
by 



i - - y i 

w (D 1 +D 9 + e+h+ri+r^-) 

D 1 Dg 

±s Z (Dp-Dg+e-hf + 4fg] V 2 (iii) 

Dp Dg 


where 



2 2 
n -rt 

— £T“ 

41 



2 2 
n tc 




(iv) 


Since eh > fg, it is noted from (iii) that Pp n > Pg n 
(n = 1,2,3,...) are positive. Do see the behavior of Pp n » Pg n 

w.r.t. ^p and Pg, we find 


dp' 


in 


dy. 


1 

z5 i 

+ 

2D 


(D^ -D 9 +e~h + — - 
1 * D. 


u 2 

Do 


1 r (D^Dp+e-itf- ~ - ~) 2 + 4fg|] 1 / 2 
1 D. Do 
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a Pto 

d? s 


1 

2D 2 



(D 1 -D 9 +e-h+ rr 
1 * D. 



£ (5-,-D +e-h+~ - + ^fg^ 1 / 2 

U 1 2 D. Do 


Since p* , and are positive for D 1 > 0 and Dg > 0, 

in d “ d Ug 

it is concluded that p^_ n increases with and Ug. Since 
v y are proportional to v 1 , Ug respectively, it may be noted 
that p^ n increases with and v g. 

Similarly it can also be shown that P g n increases with 


and Ug» 



CHAPTER V 


EFFECTS OP CONVECTIVE AHD DISPERSIVE MICRaTIOI 
OR THE LINEAR STABILITY OP PREY-PREDaTOR 
SYSTEM WITH FUNCTIONAL RESPONSE 


5.1 INTRODUCTION 


In Chapter IV, effects of convective and dispersive 
migration on the linear stability of the equilibrium state 
for two species system with mutualistic interactions and 
functional response have been investigated. Similar 
effects on prey-predator system are studied in this chapter. 

Though the study related to the evolution of intera- 
cting species has been quite extensive (May £ll3» Rapport 
and Turner £ 13 3 , Rescigno and Richardson 3 14 3 , Rosen 
CiO. Scudo and Ziegler 3 16 3 , Smith £20^), the 
investigations related to the effects of environment on 
their co-existence are of recent origin, Chewing £l3* 

Comins and Blatt [] 2 3 » Gopalsamy [] 3 3 » Hadeler et al. 3 4 3 » 
Leung 3®U» 3 9 3 » Murtrie 3 12 3 > Segel and Jackson 

3 17 3 » Segel and Levin 3 18 H> Skellam 3 19 U » Vandermeer 
3213 * The effects of environmental and ecological 
gradients on the evolution of the species can be studied by 
taking into account their dispersive and convective migration 
in the habitat. Chewing 3 ^13 » Comins and Blatt 3^H» 
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Gopalsamy £ 3 ]j , Hadeler et al. Q 4 ] , Murtrie £ 12 ]] . 

It has been noted that migration in the habitat can lead to 
stability and coexistence of the species. Levin has 

presented models of systems showing that stable co-existence 
is more probable in a heterogeneous environment than in a 
homogeneous environment. However, dispersal can give rise 
to instability when an Allee effect is present in the prey 
dynamics and the predators are highly mobile while the preys 
are relatively sluggish (Levin £ 9 , Segel and Jackson £ 17 , 

Segel and Levin £ 18 ). It has also been noted that if the 
size of an isolated terrestrial habitat is less than a certain 
critical value, the dispersing populations are open to 
extinctions which would not occur in the absence of dispersal 
(Skellam £ 19 ~2 ) * 

It may be noted here that the effects of intra- 
specific interactions with functional response in prey- 
predator dynamics have also been investigated (Kazarinoff and 
Driessche £ 6 J , Holling £ 5^ » Lin and Kahn j~lCQ)» but hhe 
effects of convective and dispersive migration on the stability 
of such systems have not been studied. 

Keeping this in view, in the following, effects of 
dispersive and convective migration of the species on prey- 
predator interactions in finite and semi-infinite habitats 
have been investigated by considering the functional response. 



102 


5.2 BASIC EQUATIONS 


Consider the interactions of prey and predator species 
in a finite or semi-infinite (0 < x < L or 0 < x < °°) habitat 
with functional response. By taking into account the 
convective and dispersive migration of the species, the 
dynamical system describing the interactions (Segel and 
Levin £ 18 ^ ) can be written as 


an. an. n p a 2 N 1 

FT + U uF = n i (a 1 -a 11 l 1 -a 12 3^7) + \ —g 

1 OX 

0N p 8N p N 1 3 2 N p 

FT” + U 2 Fx“ = N 2^*‘ a 2 +a 2l l+a^ 1 ^ + D 2 ^2~ 


(5.2.1) 


where N^(x,t), N g (x,t) represent the spatial distributions 
of prey and predator populations at time t respectively; 

, I) g are dispersal coefficients; u 1 , y g are migration 
velocities. The interaction coefficients a^, ag, a a^g 
ag^, a are positive constants. Here the predator's 
functional response to prey as suggested by Holling [(5^ 
is characterised by the term (l+aN 1 )"" and a determines 
the strength of this response. 

The nontrivial equilibrium state (N^, N*) of the 
system (5.2.1) can be found as 


a. 


Nj 




N* = — 
2 a 


cl r 


a 2 a ll 


.-—21 “1 , ( 5 . 2 . 2 ) 

12^ a 21~ a 2 a ^ ^ ^ ( a 2i" a 2 a ' J V 



J.UO 


It may be noted from ( 5 . 2 . 2 ) that both H* and K* are 
positive provided (a^-agoc) > 3-2 a ll ^ 0 * 


To investigate the local stability of the equilibrium 
state (Kj% IT*), we linearise the system (5.2.1) by writing 

i^(s,t) = H* + u(x,t) 

Ng(x,t) = + v(x,t) 


(5.2.3) 


which on substituting in the system (5.2*1) and using 
( 5.2.2 ) gives 


u 

0 ^ 


+ 


y 5 u 
1 dx 


-eu - f v + I ) 1 ng 
dm * 4 


5v u 9 v 
at + 2 ax 


gu + Dg 



(5.2.4) 


where 


e = 


h* |_a 11 (ag 1 +aga) - a q a ( a 2i“ a 2 a ) 3 
(1+aN*) ( a 2i“ a 2 a ) 


^1 

( 1 +aN*) 


g = 


a 21^2 
( 1+olT* ) ' 


and u(x,t), v(x,t) are perturbed 
From (5.2.5), it may be noted that 
but e is positive if > a Q > 

0 < a ±1 < a Q , where 

( a 9 -i “apoc ) 

a = a^a > 7 ■ ■■■■ * t >0 for 

o 1 C a 21 +a 2 a ' *“ 


population sizes, 
f and g are positive, 
0 and negative if 


a > 0 . 


(5.2.6) 
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The linearised system (5.2.4) may be associated with the 
following initial and boundary conditions; 

( i) Flux conditions 


u(x,0) = r(x), v(x,0) = s(x) x > 0 


as 


_]) = u - d = v 

i ax J J 2 ax 

a u ( ; t ) /■> a v (x,t; r-i 

ax ’ ax 


x — L for finite nabitat or x 


t > u (5.2.7) 


for infinite nabitat. 


(ii) xieservoir conditions 

u(x,0) = r(x) , v(x,0) = s(x) x > 0 

u(0,t) = u 1 , v(0,t; = v 1 t > 0 (5.2.8) 

u(x,t) -* U , v(x,t) - 0 

as x - L for finite habitat or x - 00 for infinite habitat. 

5.3 ST'ldlLllY ANALYSIS ‘rfllCdDUID MIGRATION 

If tne species are nonmigrative, the system (5.2.4) 
reduces to 


du 

dT = “ 


eu - f v 


dv 


gu 


(5.3.1) 


Considering the solutions of (5.3.1) are of the form 
exp (- xt), the eigen values X satisfy the following 
equation 


x 2 -ex+fg=0 


(5.3.2) 
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which gives 



1 

2 


e + 


1 


Le‘ 


-p -i 1/2 

fgj y 


(5.3 .3 ) 


It may be noted from (5.3.3; tnat botn A^, A g or 
their real parts are positive for e > 0 i.e. a 11 > a , but 
X ^, X 2 or ibeir real parts are negative for e < 0 i.e. 
a H< a . Thus the following conclusions may be summarized; 

(i) For e = 0 i.e. a n= a Q , both A^ and A^ are purely 
imaginary and the equilibrium state is neutrally stable . 

This is possible even when a = 0 and a ^ = 0. 

( ii ) For e > 0 i.e. a^-^ > a Q , the equilibrium state is 
asymptotically stable. This condition is satisfied even 
for a = 0 because from (5.2.6), a Q = 0, giving a ^ > 0. 

(iii) For e < 0 i.e. a^ < a , the equilibrium state is 
unstable. This can happen even when a^ = 0, but a > 0 
[2 see ( 5 .2.5) 2| • 

Results (ii) and (iii) are identical to those obtained 
by Kazarinoff and Driessche £ 6 2} • 


5.4 STABILITY ANALYSIS Wifi! DISPERSAL IN A FINITE HABITAT 


In the following, we study the stability of the equili- 
brium state of tne system (5.2.4) with dispersal in a finite 
habitat when flux or reservoir boundary conditions are 
prescribed ( v q = = °)« 
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Oase (1) Flux conditions 

In tnis case, on solving tne system (5.2.4) for ^ y 2 = 0 

with flux conditions (5.2.7) by using Laplace and finite Fourier 
cosine transforms, we obtain u(x,t) and v(x,t; as 
1 t 

u(x,t) = ^ j L A 0 ex P ( ~P i o '-*- ) + p 0 exp(-p 2U r f ) ]dl 


+ f 1 { / L A n exp (-P ln T) 

n=l o 

+ p n exp (-P 2 n' I ^ dlI} cos TT 
+ £ L^IO exp ("PlO 1 ^ + F 20 exp 


+ ^ C, p in exp ("Pln^ + F 2n 


exp (-P 2n t) I] 


(5.4.1) 


v(x,tj = f 


l J L °o exp (- p lO T; + E o exp (- p 20 ,i: ^ dI! 


2 { / C G n exp (-Pin 1 ^ 


n=l o 


B n exp (-p^HdT) = os T 2 

£ C S 10 9X P (-p 10 t) + Ggo 9X P (-PgO^l 

°o 

r 2 E>m o^C-Pm*) + G 2n exp ^n 1 ^ 


(5.4.2) 



10? 


where 


■^n 


(Pin - s 2 

p ln “ p 2n 

( p 2 4 - Pan) U - fV 

p ln *" p 2n 


0. 


(P in ~ e - D 1 <)V - g U 


1 n- 


n 


E = 
n 


in y 2n 


^+e-PpjT + gU 


2n v 


Pin ” ^ 


in ^2n 


(5.4.3) 


E 


( p ln “ D 2 a n) R n + f S 


n 


In 


2n 


G ln " 


p ln “ p 2n 

( P 2 °n ~~ p 2n^ R n ~~ £ S n 
p ln “ p 2n 

^ p ln "" e ” ^1 a n^ ^ ^n 


p m “ p 2 


n 


G 2n - 


(D 1 a n + e - p 2n ) + g a* 


p ln 

p 2n 


\ 


p ln p 2n 


1 (D, o® + Do °l + e) 


= 5 ^1 “n + "2 n 

± C(D X “*-Dg °n + ®( 2 - 4 £ sl 


1/2 


(5.4.4) 


(5.4.5) 


a n = nn/L and ( ) Q , ( ) 10 » ( ) 20 can be Gained 

from the corresponding expressions of ( ) n , ( ) ln , ( ) 2n 

by putting n = 0. She constants R n , S n are giwen by 
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K 

n 

_u 

= / 

0 

r(x) 

cos 

nxx 

L 

dx 

S 

n 

L 

= / 

s(x) 

cos 

nnx 

L 

dx 


o 


Prom (5.3,3), (5.4,5) it may be seen mat p 


10 


p go = A g 5X1(1 ^e^ce p 10 , P£0 or ^heir rea -l parts are positive 


* - and hence 

or negative depending upon whether e>0 or e<0. It 
may also be seen from (5,4.5) that when e = 0, p^ 0 , p 2 q 
are purely imaginary numbers, but p ln » P 2n (n = 1,2,3,...) 
or their real parts are positive. For e > 0, p.j_ n , p 2n 
(n = 0,1, 2, 3,...) or their real parts are always positive. 
However, when e < 0, p 1Q , p 2Q or their real parts are 
negative, but p ln , p 2n (n = 1,2,3,,..) or their real parts 

'TT ^ 

are positive whenever (— h 1 + e) > 0. In this case, if 
2 L 

(“ D^+ e) < 0, then there may exist a value of n, say n = n Q , 
L 

for which p 2n is negative. 

It may also be noted from (5.4.5) that p ln , p 2n 
(n = 1,2,3,...) always increase with dispersal coefficients 


provided (D^ ±> 2 ) is fixed. 

P or e > 0, as p , p 2 ^ (n = 0,1,2,.. •) or their 
real parts are positive, the terms in (5.4.1) and (5.4.2) 
which arise due to initial conditions of (5.2.7) tend to zero 
as t - Hence (5.4.1) and (5.4.2) take the following 

forms as t -* °° provided e > 0, 
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It u(x,t ) 

"fc -> OO 


It v(x,t) 

"t -*■ °o 


^ A^ B 
1 r- C 0 


c 50 A B 

f v , ^ "H i 2 v [—* B B — t Brcx 

L L J?10 £ 2 cP L n=l £ ) 2n“^ 


I r + fo~“J + 2 ” r- £n fn_ 

L L PlO P 2 0 J L n=l L Pln P 2 n 


(5.4.6) 


H cos 


nxx 


(5.4.7) 


It may be verified that (5.4.6) and (5.4.7) also 
satisfy the steady state form of the system (5.2.4) with 
boundary conditions (5.2.7). Thus it is seen that for 
e > 0, the nonstationary solutions (5.4.1) and (5.4.2) tend 
to the stationary solutions (5.4.6) and (5.4.7) respectively 
as t — 00 showing that the equilibrium state of the system 
is asymptotically stable. further, as P ln > p 2;n (n=l, 2,3, . . . ) 
increase with dispersal coefficients for fixed (D^~ Dg) 
including the case D^ = Dg, the degree of stability of the 
equilibrium state increases with dispersal [see (5.4.1), 
(5.4.2) 2] co-efficients. It is also noted from (5.4.6) and 
(5.4.7) that the steady state solutions form spatial pattern 
which arise due to nonhomogene ous boundary conditions* 


The following conclusions may also be drawn from 
(5.4.1) and (5.4.2): 

(i) for e = 0, the equilibrium state is neutrally stable, 
since p 1Q , p 20 are purely imaginary and p ln , p 2n (n=l,2,3, .. . ) 
or their real parts are positive. 


(ii) for e < 0, p 1Q , p 20 or their real parts are negative 
irrespective of whether Pq^> Pg J:i (n^ : l,2 > • . • ) or their real 
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parts are positive or negative and hence the equilibrium 
state is unstable. 

pa.se (ii) Reservoir conditions 

On solving the system (5.2.4) for v 1 = = 0 with 

conditions (5.2.8), we obtain u(x,t) and v(x,t) as 


u(x,t) = -| Z n {/ A) exp (-p- I) 

L n= 1 o 

+ exp (-Pg n ®) 13 di } sin ^2 


+ ^ n =l exp( ~ p ln t} + F 2n exp(-p 2n t) ^ sin £g x 


(5.4.8) 


9 °° "C _ 

V(x,t) = ^ B {/ C°n e xp(-Pi n lJ + 4 exp(-p gn l)3dl} 


sm 


n-rcx 


+ r n=l exp( - p m t} + G 2n exp(-p 2n t)3sin 


(5.4.9) 


where 




B’ 

n 


0* 

n 


E' 

n 


(Pin" 

Vn> : 

D 1 U 1 + 

f 

D 2 V 1 


P in - 

P 2 n 




- P 2 0 

D l u l 

-f 

Vi 


Pm - 

P 2 n 



(Pin' 6 

-»!<& 

D 2 V 

g 



p ln " 

P 2n 



( D l°n 

+ e -Psn> V 

1 + 

g %% 


p ln 

p 2n 




(5.4.10) 



a n = Tj^ Pin* P?n are 3- e Pl ne( l i n (5.4.5). The constants 

B4 and are given hy 

= f r(x) sin 2^2, dx 
o 

S’ = S s(x) sin dx. 

1 o 

Since p ln , p gn (n=l, 2,3, . . . ) or their real parts are 
2 

positive for (£- D^e) > 0, the terms arising due to initial 

L 

conditions in each of the expressions (5.4.8) and (5.4.9) 
tend to zero as t -* °°. Hence (5.4.8) and (5*4*9) take the 


2 

% 


following 

forms 

> as t 

-*■ OO 

provided 

. (% I^+e) > 0, 

L 


f 

( 

t 

f 

It u(x,t) 

_ 2tc 

— o 

OO 

Z n | 

r-^ 

B* 

+ p 2 3 


(5.4 

| 

.12) j 

•fc OO 

L 2 

n=l 

Pm 

P2n 



i 

| 



OO 


E ' 

S iB £g£ . 


1 

It v(x,t) 

_ 2tc 
= 

Z n j 

-°n 

L *-y 

- r-3 

(5.4 

.13) 

t OO 

L 2 

n=l 

Pin 

p Sn 






It is noted that (5.4.12) and (5.4.13) satisfy the 

steady state form of (5.2.4) with boundary conditions (5.2.8) 

snowing again the asymptotic stability of the equilibrium 

_2 

state for a _> 0, provided (~ „a,+ e ) > 0, under reservoir 

li 

boundary conditions also and the degree of stability increases 
as dispersal coefficients increase for fixed (D^-Dg) including 
the case = Dg. In this case also, spatial pattern is 

formed by steady state solutions due to nonhomogeneous boundary 
c ondit ions . 


The following conclusions may also be drawn from (5.4.8) 
and (5.4.9): 

(i) For e = 0, the equilibrium state is asymptotically 
stable as p 1r| , p 2n (n=l,2,3, . .. ) or their real parts are 
positive. Thus it may be concluded that the dispersal has 
the effect to stabilize tne otnerwise neutrally stable 
equilibrium state • 


(ii) For e > 0, the equilibrium state is stable. 

(iii) For e < 0, the equilibrium state is stable provided 
2 

(£- d + e) > 0. Thus it is noted that dispersal has the 

Ii 1 

effect to stabilize even an otherwise unstable equilibrium 

2 

state. However, when (— D 1 + e) < 0, the equilibrium state 

I* x 

may be unstable . 

It may be pointed out here that in all these cases the 
degree of stability increases with dispersal co-efficients 
provided (D^-Fg) is fixed. 


113 


'■Thus it may be remarked that the stable equilibrium 
state without dispersal remains stable with dispersal in the 
case ox finite habitat both under flux and reservoir boundary 
conditions. In the case of flux conditions, the instability 
of the equilibrium state is not affected by dispersal, while 
under reservoir conditions the dispersal has the effect -to 
stabilize an otherwise unstable equilibrium state provided 
certain condition holds. 


5.5 STABILITY ANALYSIS WITH DISPERSAL TS AH IliFBillE 
HABITAT 


In the following, the stability of the equilibrium 
state of the system (5.2.4) without convection is 

discussed in a semi-infinite habitat under flux and .reservoir 
conditions . 


Case (i) Elux conditions 

In this case, on solving the system (5.2.4) for 
with conditions (5.2.7), we obtain 


u (x,p) - p^p L -^72 exp + y l/2" e2 p(“Y2 / ' x ^-l 


£ — „ f { 

2D lVo 

exp(-y?/ 2 z) 


exp(-y^ 2 z) 


. [^(p, z+x) + E 1 (p,iz-x|)3 

^i 


z+x) + F g (p, | z-x| ) H>dz (5.5.1) 
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^’P- 1 = 575V exp(-yV 8 x) + «xp(-yg /e x)3 

1 d , y l y 2 

^ oo ©ucp ( — y ~ 2 3 

j { — [[G-^p, z+x) + G-j_(p j | z-xj ) 3 

y l 


2;D 1 D 2 o 


+ 


©ixp ^ ^ ^ Zr 3 

-T7§“ L G 2^ Pj Z+X ^ + G 2^ p ’ I z-x l ^ J }dz (5.5.2) 


y c 


where u(x,p), v(x,p ) are Laplace transforms of u(x,t), 
v(x,t) respectively and 

(p—r ) 2 y 2 ^ U - fV 
~2 

(5.5.3) 


(DpYl-p) U+fV 

a(p ) = - ~r^rr — ? B(p) 

y l y 2 


(D y -p-e) Y-gU 

C(p) = — — ; E(p) = 


y l " y 2 

(p+e-D 1 y 2 ) Y+gU 


p -,(p»x) 


P 2 (P »x) 


G 1 (p ,x) = 


G g (p ,x) 


y l “ y 2 y l - y 2 

(f 2 y i- p ) r (x) + f s(x) 


y 1 (p ) 
y 2 (p) 


\ 


/ 


y l - 

' y 2 



(p-f 2 y g ) r(x) - 

f 

s(x) 

y l - 

' y 2 



(D^-p-e ) 

s(x) 

- 

g r(x ) 

y i - 

’ y 2 



(p+e-J^yg) 

s(x) 

+ 

g r(x) 

y l " 

y 2 




(5.5.4) 


= I (£±e 

2 '•D, 


V 1 1 c 


£ ) 2 _ 
D ' 
v 2 


2 4 f g — 1 1/2 

D 1 D 2 J * 


(5.5.5) 
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It may be noted from (5,5.5) that y^ and yg or rneir real 
parts are positive for all values of p > 0 provided e is 
positive and this fact is used in obtaining the solutions 
(5.5.1) and (5.5.2). If e is negative, tnese solutions are 
not consistent with the boundary conditions (5.2.7). 

To investigate the asymptotic benavior of u(x,t) and 
v(x,t ), we use the following property of Laplace transform 
(see H 7 2 , p. 315), 


It u(x,t) = It p u(x,p). 

t — oo p -* 0 


Using (5.5.6) in (5.5.1), (5.5.2), 


1 A( 0) 

t 1 ! „ u(x,t) = exp 


B(0) r 

+ ~ ytp exp [_-y 

y<p( 0) 


we get , 

Z-i\ /z W2 

2 /2 (0)^) 


(5.5.6 ) 


(5.5.7) 


It v(x,t) 

— > OO 


1 . _0 ( 0 2 _ 

D 1 D 8 y^7 5 (0) 


exp :-y^ /z (0)3 


_E(°) 

yg /Si ( o) 


exp 


C-y^ s (0): } 


(5.5.8) 


provided e > 0, where A (0), B(0), 0(0), E(0) are obtained by 
putting p = 0 in the corresponding expressions (5.5.3) and 
y^O), y 2 (°) are the values of y 1 (p), y g (p) for p = 0. 

It can also be verified that the expressions (5.5.7) 
and (5.5.8) satisfy the steady state form of the system (5.2.4) 



116 


with boundary conditions (5.2.7). Thus it may he noted that 
for e > 0, the nonstationary solutions tend to the corresponding 
stationary solutions as t “ snowing that the equilibrium 
state remains asymptotically stable with dispersal under flux 
boundary conditions. It may also be seen from (5.5.7) and 
(5.5.8) that the steady state solutions form spatial pattern. 

Case (ii) Reservoir conditions 

As before, on solving tne system (5.2.4) for y 1 = u 2 = 0 
with conditions (5.2.8), we obtain 


u(x,p) =Tnr~r n O’(p) exp(-y^/ 2 x) + B 1 (p ) exp ( -l]/ X )U 
12 p 


+ 


f 

2D 1 D 2 J o 


{ — l~Pi (p , z+x (p , z-% ) 1 


yf* 


— p gg(p , )-g g (p . a-aOl 1 iz 

** (5.5.9) 


v(x,p) = 5 


-4-“ r o * Cp ) exp(-yy 2 x)+E'(p)exp(-yg /2 x)D 
1 D 2 p 


1 / l 


exp(-y^ 2 z) 


+ ^?2o 


jyg — Z & 1 (P ,z+x)-Gi(p,z-x)3 


— c& 2 ( P ,wx) - g 2 (p,z-x) 3) a z 




(5.5.10) 
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where u(x,p), v(x,p) are Laplace transforms of u(x,t), v(x,t) 
respectively ana 



F 1 (p,x), S'gCpjx), G- 1 (p,x), Grg(p,x) are defined in (5.5.4) and 
y^(p), y 2 (p ) are defined in (5.5.5). 

Here it should he noted that the expressions given by 
(5.5.9) and (5.5.10) are consistent with the boundary conditions 
provided e > 0. 

To investigate the asymptotic benavior of u(x,t) and 1 

v(x,t) as t — we use the property of Laplace transform given , 
by (5.5.6) in (5.5.9) and (5.5.10), for e > 0, to get 

i 

f 

It u(x,t) = { A*(0) exp C-y^ 2 (0) x] j 

t -*• °° 1 2 ; 

+ B* (0) exp C-yJ/ 2 (0) x^ } (5.5.12) | 

f 

It v(x,t) = { 0 *(0) exp C-y^CO) x"2 j 

■fc -♦ °o 2 | 

i 

+ e'( 0) exp C-yJ/ 2 (0) x] } (5.5.13)1 
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where a'(0), b'( 0), C f (0), E*(0) are obtained by patting p = 0 
in the corresponding expressions given by (5.5.11). 

Ir is seen that the expressions (5.5.12) and (5.5.13) 
satisfy the steady stare form of the system (5.2.4) with 
boundary conditions (5.2,8) showing that the equilibrium state 
remains stable for e > 0. In this case also, spatial pattern 
arises due to nonhomogeneous boundary conditions which may be 
noticed from (5.5.12) and (5.5.13). 


5.6 STABILITY ANALYSIS rfTTH CONVECTION AND DISPERSAL IN A 
FINITE HABITAT UNDER RESERVOIR CONDITIONS 

To study the effect of convection on the stability of 
the equilibrium state of the system (5.2.4) with conditions 
(5.2.8), we consider the case 



( 5 . 6 . 1 ) 


which means that the ratio of convective velocity to dispersal 
co-efficient is same for both the species. 

By using the transformations 

y l x 

u(x,t) = exp (g^-) u^x, t) 

v(x,t) = exp (g|p-) v 1 (x,t) (5. 6. &) 

2 


i 


[ 


\ 


the system (5.2.4) can he written as follows 
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0U 

aT" = - e i u i - 


2 

d u. 


f V 1 + D 1 


dx 


T 


av l 

d^T" 


g u-i 


h l V 1 + % 


0 2 T. 


Sx 


wnere 


'1 


: e + 


M 1 

4D, 


h 


1 “ 4IL 


(5.6.3) 


(5.6.4) 


The initial and boundary conditions 
the following forms 

u 1 (x,0) = r 1 (x), v 1 (x,0) = s 1 (x) 

u 1 (0,t) = U 1 , v 1 (0,t) = Y 1 

u 1 (x,t) - 0 , v 1 (x,t) - 0 


(5.2.8) reduce to 

x > 0 

t > 0 

(5.6.5) 


as x — L for finite habitat or x -* 00 for infinite habitat. 

The functions r^(x) and s^(x) are given by 

r^(x) = r(x) exp (~v^ x/2D^) 

B^(x) = s(x) exp (-Vg x / 2D 2^* (5.6.6) 

It is noted that the system (5.6.3) with conditions 
(5.6.5) is almost the same as that of system (5.2.4) for 
= Ug = 0 with conditions (5,2.8), where e^, h^, r^(x), 
s 1 (x) are defined in (5.6.4) and (5.6.6). Hence the solutions 
of the system (5.2.4) with conditions (5,2.8) are obtained as 
in case (ii) of Sections 5.4 and 5.5 by making use of (5.6.2), 
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(5.6.4) and (5.6.5) appropriately and the stability of the 
equilibrium state is discussed as before. 


For example, in tne case of finite habitat, the stability 
of the equilibrium state is determined by the nature of eigen 
values p | n , Pg n (n = 1,2,3,...) wnere 



1 ,n 2 n 2 

= S’ F“ D : 

li 


v, 2 2 

, n % 

+ rr~ 

b 


D 2 + e + W. 


+ 


2 

4D r 


2 2 


+ — \~ c£-JL_ d _ 
“ 2 L- ^ "l 


n 2 2 

-4- D r 

L 2 ‘ 


+ e + 


4D. 




(5.6.7) 


It may be noted from (5.6.7) that , pg n (n = 1,2,3,...) or 
their real parts are positive provided 



2 

D 1 + e > 0 

1 


(5.6.8) 


and hence the equilibrium state is asymptotically stable under 
the condition (5.6.8). 

From (5.6.8), we can draw the following conclusions : 

(i) For e > 0, the equilibrium state is asymptotically 
stable with dispersal and convection. 

(ii) For e < 0, the convective migration has the effect to 

stabilize an otherwise unstable equilibrium state with 

dispersal because condition (5.6.8) may be satisfied, 

2 

but y + e) < 0. 
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It may be noted from the above analysis that the 
convective migration in tne direction of some environmental 
gradient can enhance the stability of 'ere equilibrium state 
for prey-predator system with functional response. 

5.7 STABILITY ANALYSIS V.TTH COBVECIIOE AND DISPERSAL IN A 

T.*0 D IlhEES I Oh An PUTTIE HASH AT UDDER RESERVOIR COED III OKS 

In the earlier sections, we considered the migration of 
the species only along x-direction for simplicity. In real 
situations, tne species migrate in all possible directions in 
the habitat. The dynamical system governing the evolution of 
the species in a rectangular finite habitat []L, M being the 
lengths of the habitat along x,y-directions respectively]] 
can be written as 


dE 1 6E 1 dE 1 

_ + ^ — 


— E 1 (a 1 - a 11 E 


a l2 N 2 


11 1 1+ccEL 


^-i p "b p 

Sx^ ay* 2 


1 (5.7.1) 


aE 0 dEp SEp a P1 E 1 d E p d E p 

W- + av + n w- ■ + 


where D^, Dg> Dg, D^ denote the dispersal co-efficients of the 
species; u 3> y 4 denote their convective velocities* 

As in the previous section, consider the case in which 



(5.7.2) 
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i.e. the ratio of convective velocity to dispersal co-efficient 
is same for both the species along a particular direction. 

Analysing the stability of the equilibrium state 
as in tne earlier sections, it may be seen that the equilibrium 
state is asymptotically stable provided 

2 2 

+ db + 72 + rz d 3 + e ^ 0# (5.7 .3 ) 

1 o Jj M 

Comparing (5.7.3 ) with (5.6.8), it is noted that the equilibrium 
state which is unstable in one dimensional habitat can become 
stable in two dimensional finite habitat provided the inequality 
(5.7.3) holds. 

5.8 CCfCLUSIOH 

Effects of dispersive and convective migration of the 
species on the linear stability of the equilibrium state of 
prey-predator system have been investigated by taking into 
account the predator’s functional response to prey in finite 
and semi-infinite habitats. It has been shown that the 
functional response has destabilizing effect. In the case of 
a finite habitat, it has been noted that an otherwise stable 
equilibrium state remains stable with dispersal both under 
flux and reservoir conditions and the degree of stability 
increases with dispersal co— efficients when their difference 
is fixed. In this case, it is also shown that an unstable 
equilibrium state remains unstable with dispersal under flux 
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conditionsi bat dispersal has the effect to stabilize an 
otherwise unstable equilibrium state under reservoir boundary 
conditions. Even in the case of semi- infinite habitat, an 
otherwise stable equilibrium state remains so with dispersal. 

It has also been seen in the above cases that the steady state 
population of the species depenas only upon nonhomogene ous 
boundary conditions and form a nonuniform spatial pattern. 

It has been shown that the effect of convective migration 
is not only to stabilize the already stable equilibrium state 
with dispersal, but also to stabilize evenan otherwise unstable 
equilibrium state with dispersal in the case of finite habitat. 

Further, it has been pointed out that the equilibrium- 
state which is unstable in one dimensional habitat may become 
stable in two dimensional finite habitat under reservoir 


conditions 
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CHAPiSK. VI 


EFFEOfS 
L j~u slR. 


OF COuVECflVE jilTD DISPERSIVE aUGEAIIOK OR 
STABILITY OF £'a 0 OOlPEilhG SPECIES SxSIEH 
wild F LUM'O IIOITaL KCSPCIiSE 


SHE 


6.1 IE iRODUG 1'IOH 

In Chapters IV and V, effects of convective and 
dispersive migration on the linear stability of systems of 
two interacting species with functional response in cases 
of mutual istic, prey-predator interactions respectively have 
been discussed. In this chapter, similar investigations 
are carried out for two competing species model. 

The evolution, stability and continued co-existence 
of population over time of two competing species have been 
extensively studied, Pr&uer |™l]], Gomatham and Mac Donald 
[] 2 3 t Leou and lumpson []8]], May |Tll]], Maynard Smith 

Rescigno and Richardson [] 13 J . By considering a 
model of two competing preys and one predator having switching 
property of predation, Teramoto et al. [1 15]] have demonstrated 
that the predation pressure has a stabilizing effect on the 
community of competitive species. A model has been presented 
by Hallam []5]] showing that proto co-operative species can 
drive a stable competitive community to extinction. It may 
however be noted here that the effect of functional response 
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on prey-predator system has been studied, but such effects 
have not been investigated in the case of competition model 
(Kazarin off and Driessche L^J, Lin and Kahn £ 10 ^ , May M )• 

It is also seen that the effect of dispersal on the 
dynamics of biological populations has been investigated by 
several workers, Gopalsamy £ 3 , Hadeler et al. 

Harada and Pukao L 6 H j Levin |_ 9 ]] , Segel and Levin £ 14 ]] . 
Effect of boundary conditions on the dynamics of dispersive 
species in a linear habitat has been investigated by Gopalsamy 
Q 3 J and Hadeler et al. L^I] in cases of competition and 
prey-predator models respectively. Harada and Pukao £6]] 
have investigated the co-existence of two competing species 
in a linear habitat of finite length by considering Volterra- 
Yerhulst model. It should be pointed out here that the 
effects of convective and dispersive migration on the stability 
of the equilibrium state for system of two competing species 
*with functional response have not yet been investigated. 

Keeping the above in view, in the following, the effects 
of convective and dispersive migration on the linear stability 
of the equilibrium state for two competing species system with 
functional response have been studied in finite and semi- 
infinite habitats under flux and reservoir conditions. 
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6.2 BASIC EQUaIIOHS 


Consider the evolution of two competing species with 
intraspecific interactions and functional response in a linear 
one dimensional finite (0 < x < 1 ) or semi-infinite 
( 0 <_ x < °°) habitat which are capable of migration along 
x-direction. By talcing into account the effects of convective 
and dispersive migration of the species, the dynamical system 
describing their interactions (Kazarinoff and Driessche £ 7 3} » 
Segel and levin [] 14]] ) can be written as 


9E 5Bi 

r-r~ + v 

8t 1 3x 


9N 2 +U 9N 2 
FT + U 2 ” 


a. pi' d% t 

l 1 (a 1 -a 11 l 1 - 

a 9l l 1 9 2 N p 

N 2^ a 2“ a 22 N 2” + D 2 


( 6 . 2 . 1 ) 


where l 1 (x,t), Eg(x,t) represent the population distributions 
of two species at time t; , Eg their co-efficients of 
dispersal; ^ their convective migration velocities. 

The interaction co-efficients a^,a g » a n» a i2 » a 2i» a 22 a 

are positive constants. The term (1+ocE^)” stands for 
the functional response in the model (see []?]], []] 11]]] ) and 
a determines the strength of this response. 


The nontrivial positive equilibrium state (El",!*) 
of the system (6.2.1) is given by 
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a 2“ a 22 IT 2 


a 12 Xl 2 
1+aiJ j 


a 2l^l 

1+aWf 


= 0 


= 0 


( 6 . 2 . 2 ) 


Eliminat ing ]j| from (6.2.2), we get the following cubic 


equation for IT* 
, 2- T *3 


a H a 22 a 1 ' f i^“ a 22 a ^ a i a ~ 2a H^i 2 
+ l_a 11 a g2 ~2a 1 a a 22~ a i2^ a 21“ a 2 a ^ N l 

+ a 2 a i2“ a l a 22 = °* 


(6.2.3) 


It may be noted here that for a = 0, the nontrivial positive 
equilibrium state of the system (6.2.1) is stable provided 


a 22 > ^2 > fgl 

a 12 a l a ll 


(6.2.4) 


In the case a > 0, for U* and IT* to be positive, from 
(6.2.2) and keeping in view (6.2.4), we have 


a 21"’ a 2 0C 


a., 

> —■ > E* 

a ll 1 


(6.2.5) 


giving 


ag y a gOC 


( 6 . 2 . 6 ) 


It is also noted from (6,2.3) that there exists only one 
positive value for Ej and IT* under (6.2.4) and (6.2.6). 



13 0 


To analyse the local stability of the equilibrium 
state (hjj% if*;, we linearise the system (6.2.1) by writing 


ILj^Xjt) = i\ T |’ + u(x,t) 
ii g (;c,t) = 11 * + v(x,t) 


(6.2.7) 


in the system (6,2,1) and we neglect the higher order terms 
t o get 


II + y i tl - - ea - fv + D i y? 

II + v 2 li = -g^-hr+Dg 


( 6 . 2 . 8 ) 


where u,v are small perturbed population distributions and 


.* 


Li. 


e = 


f = 


g = 


(1+aIJ) 


1+aNj 


(1+aHp 2 


T~ ( a n"* a i a+ 2 a na H^) 


(6.2.9) 


* 


» h a 22 N 2 


which are positive in view of (6.2.4) and (6.2.6). 

The system (6.2,8) is associated with the following 


initial and boundary conditions: 
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(i) 


as x 
(ii) 


as x 


Flux conditions 


u(x,0) = r(x), v(x , 0) = s(x-) 


-D 


a u.( o,-t 

1 OX 




U, 


X, 22(0, t) 

J 2 ox 


V 


£u(x,t) Q 6v(x,t) ^ Q 
ox ’ 5x 


- L for finite habitat or x -* °° 


Reservoir conditions 

u(x,0) = r(x), v(x,0) = s(x) 
u(0,t) = U 1 , v(0,t) = Y 1 
u(x,t) - 0 , v(x,t) - 0 

- L for finite habitat or x — °° 


x > 0 

t > 0 (6.2.10) 


for infinite habitat. 


x > 0 

t > 0 (6.2.11 ) 


for infinite habitat. 


6 .3 


STABILITY ANALYSIS WITHOUT MIGRATION 


If the species are nonmigrative, the system (6.2.8) 
reduces to 

du = -eu-fv 

(6.3.1) 


6J 


= -gu-hv . 

Considering the solutions of (6 .3. 1) to be of the form 
exp (-It), the eigen values X^x^ satisfying the equation 

X 2 - (e+h ) X + eh-fg = 0 (6.3.2) 


are given by 



132 



= | (e+h) ± | L (e ~ h ) 2 + ^fgj 1 / 2 . 


(6.3.3) 


Since e,£,g and h are positive, ^ is always positive 
and ig is positive iff 


eh > fg 


(6.3.4) 


which from (6.2.9) gives 


a 2 2&ll-( a l- a ll^ + a ll a ^l-] ^ 1+a ®i) 2 > a i2 a gl* (6.3,5) 

Thus it is concluded that the equilibrium state (N* N*) 
is asymptotically stable provided (6.3.5) is satisfied. 


It may be remarked here that if the set of parametric 
values consistent with (6.2.4) and (6.2.6) are such that 


Ik* > 


2a 


11 


(6.3.6) 


then (6.3.5) holds automatically. 

Sl-i 

As — > if [see (6.2.5)[J, it is also remarked that 

a ll 1 

the inequality (6.3.5) may not be valid for a particular set of 
parametric values, though the conditions (6.2.4) and (6.2.6) 
are satisfied showing that the functional response has 
destabilizing effect. 


6 .4 STABILITY ANALYSIS WITH DISPERSAL IN A FINITE HABITAT 

In the following, we examine the stability of the 
equilibrium state with dispersal in a finite habitat under 
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flux and reservoir boundary conditions ( y i =V 2 =( ^* 

Case (i) Blux conditions 

On solving the system (6.2.8) for ^=1^=0 with flux 
conditions (6.2.10), we obtain u(x,tj and v(x,t) as 
t 


u 


(x,t) = £ J [A o exp(-p 10 i 1 ) + B o exp(-p 20 T) H d2 


+ ~ Z /f J lA n exp(-p tn I) + B n exp(-p 2n r i)I]dl} cos 


n=l o 


In' 


2n 


nigc 

T 


+ j [l 10 exp(-p 10 tj + ? 20 exp(-p 20 t)d 


+ I £ . t? ln exp(-P ln t) + P Sn exp(-p 2ll t)3oo s SjS (6.4.1) 


nrac 


n=l 


in'- 


t 

v(x , t ) = ~ f [p o exp(-p 10 l) + B o exp(-p 20 l)3d'Il 
o 

o° ~fc 

+ r 2 { J lP n e xp(-Pi n T ) + 3 n ex P(-p 2 n ,j: ^ dT} cos 
n=l o 

+ r P 1 0 ex P(-P 10 ' fc ) + 8 20 exp(-p 80 t:) ^ 

oo __ 

+ f S + S 2 n exp (- p 2n t) 3° 0s TT 

n=l 


nux 

L 


( 6 . 4 . 2 ) 


where 
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A = 
n 


B = 
n 


(Pm-n-Vn) u + 


c 


p ln“ p 2n 

(°8°n +h -Pgn) V - fT 
p ln~ p 2n 

(Pm- e '¥l) v + en 


n 


p ln“ p 2n 


E 


^i a n +e p 2n^ ^ ^ 


n 


E = 
In 


E 


G 


2n 


In’ 


S 2iT 

/ \ 

I p ln 

P 



p ln~ p 2n 




(p 

In-b-Dg^n^ R n 

+ 

f 

S n 


p ln“ p 2n 




(D 

2 CT n + h_p 2n ^ R n 

- 

f 

S n 


p ln~ p 2n 




(P 

ln“ e ”' D i a n^ S n 

+ 

g 

E n 


p ln“ p 2n 




CD 

!°n +e " p 2n) S n 

- 

g 

R n 


p ln“ p 2n 


= I ( D 1 a n + D 2 a n + e + 

1*2 nj 

+ i r (D. a 2 _ D 0 a 2 + e - h) 2 + 

~ 2 ^ 1 n 2 n J 


°n ~ IT 3X1(1 ( 'o’ ^ ^10’ ( ^20 3X6 

corresponding expressions of ( ) , ( 

putting n = 0. The constants R , S. 


(6.4.3) 


(6.4.4) 


4fg^ 1 / 2 (6.4.5) 

obtained from the 

bn’ < bn ^ 

are given by 
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R, = J r(x) cos dx 


S = J s(x) cos -g- dx. 


It may be nored from (6.3,3) and (6.4.5) that p 10 = 

PgQ = It may also be seen from (6.4.5) that p ln (n=l,2,3, 

are always positive and p (n=l,2,3 , . . . ) are positive if 


D 1 tc D 9 it 

( — 5 — + e) (— 5 — + h) > fg. 


(6.4.6) 


It has also been shown that p ln , Pg n (n = 1,2,3,...) increase 
with dispersal co-efficients (see appendix, Chapter IV). 

If (6.3 .4 ) ho Ids i.e. p p , p 2 n ^ 6 , 1,2, 3,... ) ar e 

positive, the nonstationary solutions (6.4.1) and (6.4.2) 
converge to the following forms as t — 

A T 5 00 A "R 

t 4 _ / \ 1 r" 0 , 0 — 1,2 « r— n , n — « _ xitex / r A n * 


C 

it T(x,t) = £ 
t - 00 ^ ‘-Pk 


B 

0 - 

1+ - 

OO 

r 4_ 

p 20" 

-l + 1 

h \ 

n=l 

Lp m 

f°_- 

1 , 2 
J + L 

00 

2 I 

n=l 

r-4_ 

p 20“ 



Since Pp n » pg n (n-l,2,3, . . . ) are positive and increase with 
dispersal co-efficients, the rate of convergence increases 
with dispersal co-efficients. Further, as (6.4.7) and (6.4.8) 
satisfy the steady state form of the system ( 6 . 2 . 8 ) with 
boundary conditions ( 6 . 2 . 10 ), the equilibrium state is 
asymptotically stable for a. > 0 and the degree of stability 
increases with dispersal co-efficients. It can also be seen 
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from (6.4.7) and (6.4.8) that the steady state solutions form 
spatial pattern which arise due to non homogeneous boundary 
conditions . 

Sven if the inequality (6.3,4) does not hold, the 
equilibrium state c&n become stable provided = Dg and 
the initial and boundary conditions satisfy the following 
relations 

( h- X 9 ) U = f V 

(6.4.9) 

or (e~* 2 ) V =• gU 

(h- x 2 ) r(x) = f s(x) 

(6.4. 10) 

or (e- A g) s(x) = g r(x ) 

which are likely to be satisfied under very special ecological 
circumstances, where x g is defined in (6.3.3)'. In general, 
such conditions may not be applicable to realistic habitats. 

Case (ii) .Reservoir conditions 

Again in this case, on solving the system (6.2.8) 
with conditions (6.2.11), we obtain 

r, 00 t 

u(x,t) = ^| 2 n {/ £ a ' exp(-p 1;n 'J:) 

L n=l o 

+ exp(-p 2n !T)^dT} sin ££2 

oo 

+ ~ n ?. 1 ^'ln exp ^“ p ln t ^ + ■ F 2n exp ^“ p 2n t ^ sin L" 


(6.4.11) 
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Jj 


J r(r) sin dx 


n 


1 

S ’ = J s(x) sin 
o 


n-nx 


dx . 


If (6.4.6) is satisfied, as p^, Pg n ( n =l5 2,3 , . . . ) 
are positive, the nonstationary solutions (6.4.11) and (6.4.12J 
converge to the following forms as t — °° 


It u(x,t) 

± -> CO 


2 % v 
P Zu 

3j n=l 


n 




r -- + 

u p 


B' 

n 


In 


2n 


■ J sin 


nmx 

L 


It u(x,t) = ^4j- Z n 
t - oo L n=l 


O' E 

r_iL. + J2_j 

Lp ln p 2n J 


sin 


n-rcx 

L 


( 6 .4 . 15 ) 

(6.4.16) 


and the r a te of convergence increases with dispersal 
co-efficients. Since (6.4.15) and (6.4.16) are steady state 
solutions of the system (6.2.8) with boundary conditions (6,2.11), 
the equilibrium state is asymptotically stable for a > 0 and 
the degree of stability increases with dispersal co-efficients 
under reservoir conditions also. It is further noted that the 
steady state solutions give rise to spatial pattern. 

Keeping in view the inequalities (6.3.4) and (6.4.6), 
it is observed that the stable equilibrium state without 
dispersal remains so with dispersal and the effect of dispersal 
is also to stabilize even an otherwise unstable equilibrium 
state without dispersal under (6.4.6). 



139 


Even if the inequality (6.^.6) does not hold, the 
equilibrium state can become stable provided = D p , tne 
initial and boundary conditions satisfy vhe relations such 
as (6.4.10) and 


(°-V u i = 

or (e- J g J = giq 

which are not applicable to realistic situations. 


(6.4.17) 


6.5 STABILITY ANALYSIS WITZ DISPOSAL H .All ILEIhlTE 
HABITAT 


In the following, the stability of the equilibrium 
state without convection ( p q = v g = °) is analysed in a 
semi-infinite habitat under flux and reservoir conditions. 

Case (i) Blux conditions 

On solving the system (6.2.8) for = y g = 0 with 
conditions (6.2.10), we obtain 


= d^p { ^ ( p 7 exp 

+ -I^fr e x pC-yo 1/2 (p)xn} 
yg 7 (p) 2 

exp E-yi /2 (p)€l 


rV 2 / 


2D 1 D 2 o 


I i 


exp [ 2 -y p 72 (p ) z] 


yg(p) 


7 172 ' [| 1 (p,z+x)+P 1 (p,jz-x| )3 

yf (p) 


Jl? p (p,5*x) + P p (p,i z-x| )3>iz 


(6.5.1) 
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v(x,p) 


_!_{ -2ie2_ 

TyKp) 


exp 


Lry^tpJx J 




20 ,D 


/ { 


12 o 


+ - ffe '} — exp E-yg //2 (p)x J} 

Yg 7 (p) 
exp [^-y^ //2 (p )z] 

~F^PT 


C g i(p >^x) + G- (p , J Z-X| ) ] 


+ 


expEI-y 2 /2 (p)i3 _ . . 

T/g — — [&g(p ,z+x) + Gp(p , z-x ) 3 > dz 

yi 7 (p) 2 2 


(6.5*2) 


where u(x,p), v(x,p) are Laplace transforms of u(x,t), 
v(x,t) respectively and 

(0 2 y r p-h)U + fV (p+h-0 2 y 2 )U _ fV 


A(p) = 


y l- y 2 


B(p) = 


y l- y 2 


(D^-p-e)? + gU (p+e-D y )V 

0(p) = - - 1 - 5 E(P) LJS ~ 


- gu 


^(PjX) = 


p 9 (p,x) = 


Gi(p»x) = 


&g(p»x) = 


y l~ y 2 

( OgYp— p~h) r (x ) + f s(x) 


y l” y 2 


y l“ y 2 


(p+h-0 2 y 2 ) r(x) - f s(x) 

(■ D l y l“ p " e ^ sW + s r ( x ) 


y l“ y 2 


(p+e-D^yg) s(x) - g r(x) 


y l- y 2 


(6.5.3) 


\ 


(6.5.4) 
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lx may be noted from (6.5.5) that y^ and y 9 are 
positive for all values of p >_ 0 provided eh > fg and 
this fact is used in obtaining the solutions (6.5.1) and 
(6.5.2). 


To know the asymptotic behavior of u(x,t ) and v(x,t ) 
as t — the following property of Laplace transform 

It u(x,t) = It pu(x,p) (6.5.6) 

t - °° p - 0 

is used in (6.5.1) and (6.5.2) to get 


It u(x,t) 
t — °° 


D 1 D 2 


{ 


A.( 0) 

y^(oj 


exp 


c-yi /a (o)xn 


It v(x,t ) 

"t oo 




+ 


B(0) 

y^*(o) 


exp 


C-yg /8 (°) =0 > 


.. Slfi.1 - 

y5p 8 (o) 


exp 


L-y;[ /2 (&.) 


exp C-y| /2 (o) O i 


(6.5.7 ) 


(6.5.8) 


where y 1 (0), y g (0) are values of y 1 (p), y g (p) for p = 0 
and A(0), B(0), 0(0), E(0) are obtained by putting p = 0 
in the corresponding expressions (6.5.3). 


It can be verified that the expressions (6,5.7) and 
(6.5,8) satisfy the steady state form of the system (6.2.8) 
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with boundary conditions (6.2,10). Thus it is seen that 
the nonstationary solutions tend to the corresponding 
stationary solutions as t - , showing the asymptotic 

stability of the equilibrium state under flux boundary conditions 
for eh > fg. It is also observed from (6.5.7) and (6.5.8) 
that spatial pattern is formed by steady state solutions. 


Case (ii) Keservoir conditions 


In this case, on solving the system (6.2.8) for 2 =0 

with conditions (6.2.11), we obtain 


u(x,p) = *A'(p) exp G-yJ //2 (P ) sQ 


'1^2- 


+ B'(p) exp Lr-yJ/ 2 (p)x_| } 


+ 


x “ expL-yJ /2 (p)0 


2 % D 2 { 1 

. 1 / 

2 

TTZ 7; 


— L?l(p >z+x)-I' 1 (p,z-x) 13 


exp L~yp //2 (p )d _ 

+ — — — T V S [j 2 (p,z+x)-P 2 (p,z-x)3}(iz (6.5.9) 


y . 7 *(p ) 


v(x,p) = ^C'(p)exp t-y^ 2 (p) xU 

1 2 


+ S’(p)exp£-yy 2 (p) £j} 


+ 




1 expH-yJ^tp)!] 


/ { 


2 - D l I) 2 o 


expE-y£ /2 (p)H 


L g i(p » z+x ) _G i (p » z - x )U 

y/ (p) 


~-T7V OMp ,z+x)-G p (p ,z-x) ~j }dz (6.5.10) 

yl /2 (p) 2 2 
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where u(x,p), v(x,p ) are Laplace transforms of u(x,t), 
v(x,t) respectively and 


(P ) 

B'(p J 

C'(p) 

E‘(p) 


( Bgy -j-d— h )d^u^+ 

f 

D 2 V 1 

y r y 2 



(p+hr-Bgyg)!)^- 

f 

DgVf 

y 2 - y 2 



(D iy i-p-e)D 2 v 1+ 

g 

3 1 U 1 

y x - y g 



(p+e-D^gjDgV^ 

g 

D l u l 


yi- y 2 


(6.5.11) 


To investigate the asymptotic behavior of u(x,t) and 
v(x,t) as t using (6.5.6) in (6.5.9) and (6.5.10), we get 

It u(x,t) = <A*(0) exp G-yi /2 (°) X H 

t -* oo ^±2 

+ B’(0) exp Q-y;j/ 2 (0) x]]} (6.5.12) 

It v(x,t) = “j- {C’(0)expL'"yi /2 (0) x 3 

t - 00 12 

+ E*(0) exp [}*yg //2 (0) x^ } (6.5.13) 

where A* (0), B’(0), G’(0), E'(0) are obtained by putting 
p = 0 in the corresponding expressions (6.5.11). 

Since the expressions given by (6.5.12) and (6.5.13) 
satisfy the steady state form of the system (6.2.8) with 
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reservoir conditions (6.2.11), the asymptotic stability of 
the equilibrium state is established in this case also 
provided eh > fg. It can also be seen that the steady state 
solutions give rise to spatial pattern due to nouhomogeneous 
boundary conditions. 

Thus it is concluded that an otherwise stable equilibrium 
state without dispersal remains so with dispersal in the case 
of infinite habitat. 

6.6 SUABILITY ANALYSIS itfllH CONVECTION AND DISPERSAL UNDER 

RESERVOIR CONDITIONS 

To investigate the effect of convective migration on the 
linear stability of the equilibrium state of the system (6.2.8) 
under reservoir conditions (6.2.11), consider the case for 
which 



i.e. the ratio of convective velocity to dispersal co- 
efficient is same for both species. This is biologically 
reasonable because the species having more dispersal ability 
are likely to have faster convective migration velocity. 

Using the transformations 


u(x,t) = exp (y 1 x/2D 1 ) u^(x,t) 
v(x,t) = exp (y 2 x/2D g ) v 1 (x,t) 


( 6 . 6 . 2 ) 
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the system (6,2.8) is reduced to the following form 


dir, a 

FT" = “ e i u i “ fT i + J i 

3n 

a 2 v^ 


(6.6.5) 


6v 


1 


= -gu. - h 1 v 1 + D c 


ot ' 11 2.2 

ax 


where 


6^=3+ y 2 /4D 1 ; h 1 = h + y|/4D 9 . 


(6.6.4) 


the initial and boundary conditions (6.2.11) are modified 
as follows s 

u 1 (xjO) = r 1 (xjj v 1 (x,0) = s 1 (x) x > 0 

U;] _(0,t) = u x , v^Ojt) = v 1 t > 0 

u 1 (x,, t) - 0 , v^(x,t) - 0 (6.6.5) 

are x -* L for finite habitat or x — °° for infinite habitat. 
The functions r^(x) and s^(x ) are given by 


r 1 (x) = r(x) exp (-Ujx/2 i 1 ) 
s 1 (x) = s(x) exp (-VgX/^ig) * 


( 6 . 6.6 ) 


It is noted that the system (6.6.3) with conditions 
(6.6.5) is of the same form as that of the system (6.2.8) for 
V 1 ~ y 2 = 0 w ith conditions (6.2.11) where e^ f h^,r^(x), 
s^(x ) have been modified according to (6.6.4) and (6.6.6). 
Therefore, the solutions of the system (6.2.8) with (6.6.1) 
are obtained from case (ii) of Sections 6.4 and 6.5 by making 
use of (6.6.2), (6.6.4) and (6.6.5) appropriately. 
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Eor example, in the case of finite habitat, the solutions 
of the system (6.2.8) with conditions (6.2.11) for the case 
of (6.6.1) can be obtained xrom (6.4.11) and (6.4.12) by- 
multiplying tnese equations with exp( v ^x/2D^) and e:xp ( Vgx/2Dg ) 
respectively and replacing e by e^, h by h^, r(x) by 
r 1 (x), s(x) by s 1 (x). Ihus the corresponding eigen values 
p^ n , p^ n (n=l, 2,3 , . . . ) for the system (6.6.3) with (6.6.2) 
are given by 


Pin 

Pin 


1 n 2 -^ 2 

± f£L2L_ b + 

2 1 t 2 D l + 

h 


2 2 

TT" D 2 + a l + h l^ 


9 9 

1 1 C (- *- 


2 2 


”2“ D i“ £ rf~ D 2 + e l~ V + 4f Sj 1/2 * (6.6.7 ) 
LI 


It is seen from (6.6.7) that 

p ositive provided 
2 

(jT D i + 0 1 + e) ( £2 V 




P in (t 1 — 1 , 2 , 3 , • . • ) are 


+ h) > fg. (6.6.8) 


It has also been shown that (see appendix, Chapter IV) 

Pin* P^n ( n= l» 2 >3,.. . ) increase with y ^ and 

Keeping in view the conditions (6.4.6) and (6.6.8), it 
is concluded that even an otherwise unstable equilibrium state 
with dispersal can become stable due to convective migration 
and the degree of stability increases with convective velocities. 
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6 .7 STABILITY ANALYSIS WITH COHVEC1IOE AED DISPERSAL Hi A 

lVt r O DIITES SIOL AL FINITE HABITAT UKHESR RESERVOIR CONDITIONS 


Though «7 e considered the migration of the species only 
along x-dir action in the earlier sections, in real situations, 
migration tales place in all possible directions. In such 
cases, the equations governing their evolution in rectangular 
habitats *{[[1,1.1 being the lengths of the habitat along 
x,y-dir actions respectively [} are written as 


5N 1 oil 9EL 

__A + y _ — £ + p„ — -£ 
3 ^ + 1 5x 3 ay 


3N dEp QNp 

sr + u e ax + v 4 ay 


3 i 2 H, 


6 


a 1 pNp \j juv jy -Lv *4 

( a i“ a n N i“ T +a E ^ + D i ” 2 + ^3 ~ v 2 

1 3x 3y 

(6.7.1) 

a p1 N 1 3 2 N p 3 2 Ep 

Eg ( a 2 “ a 22 H 2“ l+aE 1 ^ + D 2 ~^2~ + ^4 ^2“ 


where D^, Dg, D^, D^ are co-efficients of dispersal of the 
species and y^, y 3 , y^ their convective velocities. 

Analogous to Section 6*6, consider the case 



i ,e . the ratio of convective velocity to dispersal co-efficient 

is same for both the species along the particular direction. 
Analysing as before, it may be observed that the equilibrium 
state is asymptotically stable under reservoir conditions 
provided 
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2 2 K 2 2 

( n + 5— 7 ) + + iVlL- q + 2 

\2 J l + tt 2 "V £ET + 4D„ ; S 2 J 2 , t 2 

1 1.1 1 3 Jj ii 


+ + -li 

T 4Dg + 4D 4 


) > fg . 

(6.7.3) 


Comparing (6.7.3) with the condition (6.6.8), it is 
concluded tnat the equilibrium state which is unstaole in 
one dimensional finite habitat may become stable in two 
dimensional finite habitat due to convective and dispersive 
migr at ion . 


6.8 CONCLUSION 

Effects of convective and dispersive migration on the 
linear stability of the equilibrium state for two competing 
species system with intraspecific interactions and functional 
response have been examined. It has been shown that the 
functional response has destabilizing effect. In both finite 
and semi-infinite habitats, it has also been shown that the 
stable equilibrium state without dispersal remains so with 
dispersal under flux and reservoir conditions. In all these 
cases, it has been pointed out that the steady state populations 
form spatial pattern due to nonhomogene ous boundary conditions. 

In the case of finite habitat, it has been noticed that 
the degree of stability increases with dispersal co-efficients 
in absence of convective migration and dispersal may stabili z e 
even an otherwise unstable equilibrium state without migration 
under reservoir conditions. It has also been pointed out that 
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the effect of convective migration of the species is to 
stabilize even an otnerwise unstable equilibrium state 
with dispersal and the decree of stability increases with 
convective velocities* It has been nored tnat the 
equilibrium state wnich is unstable in one dimensional 
finite habitat may become stable in two dimensional finite 
habitat under reservoir conditions. 
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CHAPTER VII 


EFFECTS OF^DISPERSAL OI'I THE LINEAR STABILITY OF PREY- 
PREDATOR SYSTEM v,ITH FUNCTION Al RESPONSE III A 
CIRCULAR FINITE HABITAT 


7.1 INTRODUCTION 

In Chapter V, effects of convective and dispersive 
migration on tne linear stability of the equilibrium state 
for prey-predator system with functional response have been 
examined in one and two dimensional linear habitats. As a 
continuation, in the present chapter, similar study is 
conducted on prey-predator interactions with dispersal in a 
circular finite habitat. 

Extensive research has been done on the evolution and 
co-existence of interacting species by considering temporal 
variations of the species, Holling [] 3]] » Kazarinoff and 
Driessche [] 4 ]] » H* 1 ahd Kahn []?]] , Rapport and Turner co- 
It should, however, be noted that the effects of dispersal on 
the stability of the equilibrium state of systemsof interacting 
species have also been investigated, Comins and Blatt j“l]], 
Hadeler et al. []2]], Leung Levin [] 6 ]} , Murtrie [] 8]] , 

Segel and Levin £ 10]] , Skellam [] 11]] . In particular, 

Comins and Blatt CO have pointed out that dispersal can lead, 
to stability and continued co-existence of the species. It may 
also be noted here that the equilibrium state of systems of 
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intar acting species with functional response are not always 
stable (Kazarin of f and Driessche kin and Kahn [] 7 ) 

and the effects of dispersal on the stability of such systems 
have not been investigated. 

Keeping all these in view, in the following, we investi- 
gate the effects of dispersal on the linear stability of the 
equilibrium stare for prey-predator system in a circular 
habitat by considering intraspecific interactions and functional 
response. 


7.2 BASIC EQUATIONS 


Consider the evolution of two interacting and dispersing 
species forming a prey-predator system in a circular (0 < r < R) 
habitat. By taking into account the effects of functional 
response, the equations governing the evolution of the species 
can be written as, 


T 1 


a ^2 




) + 


D ( x 

r ar ^ 


, T / _ , a 21 ij l 

- Ng(-a 2 + 




d_ 

ar 


(r 


an, 

c 

dr 


an. 

Fr“ 


(7.2.1) 


where H^(r,t), N g (r,t) represent the population distributions 
of prey and predator species respectively at time t; D^, Dg are 
the co-efficients of dispersal? and the interaction co-efficients 
a l , ag » a n > a ig » ^ i > a ar e positive constants. 


The term 
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— 1 

(1 + ccli^) characterizes the functional response in the 
system end a determines its strength. 


She nontrivial equilibrium state (H*, I'T*J of the system 


(7.2.1) is given by 


lC = 


3-r 


1 •"•ctgOC * 


o I'J”^ — 

9 (TO 


a 21 La-i(a 2 i- a 2 a) - a 2 a lx J 


a 


12 


( e* 2 -^""agCX ) 


(7.2.2) 


It is noted from (7.2.2) that both N* and IT* are positive 


provided 


ag-^^agOC ) > a 2 a ll * (7.2.3) 


To examine the local stability of the equilibrium state 
(N*, N*), we linearise the system (7.2.1) by writing 

U 1 (r,t) = IT* + u(r,t) 

h 2 (r,t) = H* + v(r,t) 

which on substituting in the system (7.2.1) and using (7.2.2) 
gives 


(7.2.4) 


5 u 

FF 

9v 

FE 


-eu - fv + r g r 


19 /„ du 


(r 22) 
' dr J 


= gu + B 2 I fj (r 0) 


(7.2.5) 


where 


N* &il( a 21 +a 2 cx ) “ a l a ( a 21 _a 2 a ) H 
(1+aH*) ^ a 2l a 2 a '* 


f = 


a 12 N 1 


a 2l H 2 


(1+ocN*) (1+a N*)‘ 


(7.2.6) 
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From (7.2.6), it may be noted that f and g are positive 
under (7.2.3). ihe constant e is positive for > a Q 

and negative for < a c where 


( Q*p i ap£% ) 

a = a 1 a — — — > 0 . 


(7.2.7) 


It should be here noted that the equilibrium state 
(H*, E*) is asymptotically stable in absence of dispersal 
(refer Chapter V ) for a > u provided 

e > 0 . (7.2.8) 

She following initial and boundary conditions are 
associated with the system (7.2.5). 

(i) Reservoir conditions 

u(r, 0) = u Q (r), v(r , 0) = v Q (r) 0 < r < R 

u(R,t ) = U 1} v(R, t ) = V 1 t > 0 (7.2.9) 

u(r,t) and v(r,t) are finite at r = 0, 


In physical situations, this kind of problem arises when 
the population density of the dispersive species is prescribed 
at the boundary . 


(ii) 


Efflux conditions 


u(r,0) = u n (r), v(r ,0) = v n (r) 

V, 


du 

aF 


o^"'’ ’ ' ~ ' o' 

5v 


+ ku = Ug , + kv 


0 < r < R 

t > 0 
at r = R 


u(r,t) and v(r,t) are finite at r = 0 . 


(7.2.10) 
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She boundary conditions (7.2.10) indicate that the 
flux of the population density at the periphery of the 
circular habitat depends upon density of tne population at the 
bo undeny. 

In the following, the linear stability of the equilibrium 
state (R*, h* ) is studied under (7.2.3). 


7.3 SSABILISY .ANALYSIS WISH DISPERSAL OH HSR RESERVOIR 
COL'DISIOKS 


By using Laplace and finite Hankel transforms, the 
system (7.2.5) is solved with conditions (7.2.9) in a circular 
habitat (0 < r <_ R) to obtain u(r,t) and v(r,t) as 


u(r , t ) 



x. 


J o< r 
n J 1 (R x n ) 


t 

f H A n exp (-p ln S ) 
o 

+ B n exp (-p 2n S)^dS 
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~2 1 
k n 


J cV *0 
[?!<* x n)3 2 


CP ln exp (-p ln t) 

+ P gn exp (~P Sn *)J 


v(r , t) 


2 

E 


2 

n 


x. 


n 


J p( r *n> 

q(R * n ) 


t 

/ C c n e ^P (~ p ln T ^ 
o 

+ 3 n exp (-p gn I)]ai 


a . 

R S n U 1 (H x n ) 2 2 


C S ln exp (-p ln t) 

+ S gn exp (-PgPOQ 


(7.3.1) 


(7.3.2) 


where J ffi (x) 


is Bessel function of order m and the summation 
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is taken over all the positive roots, x n , of the transcendental 


equation 

J 0 (R x) = 0. (7.5.3) 

The constant co-efficients in (7.3.1) and (7.3.2) are defined 
as follows: 
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(7.3.4) 


(7.3.5) 


(7.3.6) 
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where u Qn and v are finite Ranked transforms of the 

first kind defined as 

R 

u on = ■>' E u o< r) J o (r ^ 
o 

(7.3.7) 

R 

v = J r v (r ) J Q (r x n ) dr. 
o 


From (7.3.6) it may he noted that p ln and p 2n or 
their real parts are positive for all values of n if 

min (D-x? + e) > 0. (7.3.8) 

n 

It is also seen that p ln and p gn or their real parts increase 
with dispersal co-efficients for fixed (D^- Dg). 

If the inequality (7.3.8) holds, i.e. p^ n , Pg n (n = ^>2,3,...) 
or their real parts are positive, u(r,t) and v(r,t) given by 
(7.3.1) and (7.3.2) take the following forms as t -*■ 


lim 

"fc oo 


lim 

t -* OO 


u(r, t) = 


8 „ J oU 


J o (r 





(7.3.9) 


(7.3.10) 


It can be easily verified that the expressions (7.3.9) and 
(7.3.10) represent the steady state solutions of the system 
(7.2.5) with boundary conditions (7.2.9). Since the 
nonstationary solutions (7.3.1) and (7.3.2) tend to the 
corresponding stationary solutions (7.3.9) and (7.3.10) as 
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t - the equilibrium state is asymptotically stable under 
(7.3.8). It may also be pointed out that the degree of 
stability of the equilibrium state increases with dispersal 
co-efficients provided (D^- Dg) is fixed. 

Keeping in view of the condition (7.3.8), the following 
conclusions can also be drawn : 

(i) When e = 0, the equilibrium state is asymptotically 
stable. Hence dispersal stabilizes an otherwise neutrally 
stable equilibrium state. 

(ii) When e > 0, the equilibrium state is stable. Hence the 
stable equilibrium state without dispersal remains stable with 
dispersal. 

(iii) When e < 0 i.e. a ±1 < a Q , the equilibrium state is 
stable if (7.3.8) is satisfied. 

Finally, if (7.3.8) does not hold, then there exists 
least one positive value of n, say n = n Q , for which pg n is 
negative and in that case the equilibrium state is unstable. 

7.4 SUABILITY ANALYSIS WITH DISPERSAL HINDER EFFLUX 
CONDITIONS 

As before, on solving the system (7.2.5) with conditions 
(7.2.10) in a circular habitat (0 < r < R), u(r,t) and 
v(r,t) are obtained as 
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u(r,t) = | E I n (r) / exp (-p^ 1 ) + exp(-p^ n T)I]ai 


P X n( r ) 

+ 2 — — — [> ln exp (-P x t) 

R n J 0 (R y n ) ln ln 


+ Rg n ex P (~P2n C 7 * 4 * 1 ) 


v(r,t) = | S I n (r) / C c n exp ("Pin ^ + E n ex P(“ p 2n 


y -- fi, 

+ ,2 L 
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In( r ) 


n J q (R y n ) 


C G i n ex P (“Pin ^ 


m 


+ Gg n exp (-Pg n t)]] (7.4.2) 

where the summation is taken over all positive roots, y n , of 
the transcendental equation 

k J 0 (R y) = y J 1 (E. y). (7.4.3) 

The function I n (r) is defined as 


T / \ J p( r y n^ y n 

” " J 0 ( a y c > ' 


(7.4.4) 


The constant co-efficients in (7.4.1) and (7.4.2) are given as 


follows : 
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where u^, are 

B 

u on = { r “o<*> J P r y P dr 

(7.4.8) 

R 

v on = { r v o< rJ J o (r y n> dr - 

From (7.4.7), it may be noted that pj_ n , Pg n or their real 
parts are positive for all values of n provided 


min (D 1 yj? + e ) > 0. (7.4.9) 

n -*• “ 

If (7.4.9) does not hold, then there may exist least positive 
value of n, say n = n Q , for which pg n is negative. 

, Pg n or their real parts 
are positive for all values of n, u(r,t) and v(r,t) given by 
(7.4.1) and (7.4.2) respectively converge to the following 
forms as t co, 


If (7.4.9) is satisfied i.e. p 1 


It u(r,t) = | Z I (r ) 
t - °° * n n 



LPln P2n-i 


It v(r,t ) = | 2 I n (r) 
t -» °° n 



P2n- 


(7.4.10) 


(7.4.11) 


As the expressions (7.4.10) and (7.4.11) satisfy the steady 
state form of the system (7.2,5) with boundary conditions 
(7.2.10), the equilibrium state is asymptotically stable under 
(7.4.9). It may also be noted that the degree of stability of 
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the equilibrium state increases with dispersal co-efficients 
for fixed (D 1 - Dgj. 

Keeping in -view of the inequality (7.4.9), the following 
conclusions are also noted ; 

(i) When e = 0, the equilibrium state is asymptotically 
stable. Hence dispersal has the effect to stabilize 
an otherwise neutrally stable equilibrium state. 

(ii) When e > 0, the equilibrium state is always stable. 

Hence dispersal stabilizes further the already stable 
equilibrium state without dispersal. 

(iii) When e < 0, the equilibrium state is asymptotically 
stable provided (7.4.9) is satisfied. If the inequality 
(7.4.9) does not hold, the equilibrium state may be 
unstable. 

7.5 CONCLUSION 

In a circular habitat, effects of dispersal on the 
linear stability of the equilibrium state for prey-predator 
system with intraspecific interactions and functional response 
have been investigated. Both under reservoir and efflux 
boundary conditions, it has been shown that the equilibrium 
state is asymptotically stable for e > 0. Thus it is concluded 
that the equilibrium state which is stable without dispersal- 
remains so with dispersal. However, it has also been pointed 
out that even an otherwise unstable equilibrium state can 
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become stable with dispersal provided certain inequalities 
such ms (7.3.8), (7.4.9) involving dispersal co-efficients 
hold. it has also been observed that the degree of stability 
of the equilibrium state increases with dispersal co-efficients 
if tn.eir difference is fixed. 
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CHAPTER VIII 


EFFECTS OF DISPERSAL OH THE LIREaR STABILITY CP TtfO 
COx£PETIHG SPECIES SYSTEM IH A CIRCULAR 
FIR I IE HADIIAI 


8.1 lH IRODUC T IOR 

The linear stability of the equilibrium state of system 
of two competing species with functional response have been 
studied in Chapter VI by considering the effects of convective 
and dispersive migration of the species. In the current 
chapter, similar analysis is carried out on dispersive compe- 
titive species in a circular habitat. 

Since the days of Volterra, the evolution and co- 
existence of interacting species have been studied by ignoring 
the effects of migration, Brauer HID* Gomatham and Mac Donald 
C 23 . Leou and lump son [2 6]]}, Lin and Kahn [] 8]] • The effects 
of dispersal on the dynamics of biological species have also 
been investigated more extensively, Gopalsamy []] 3 ]]} , Hadeler 
et al. [] 4 ]] , Harada and Pukao [] 5]]] , Levin £] 7 ]]} , Segel and 
Levin []9]]}, Skellam [] 10]]] . In particular, Hadeler et al. 
£43 Lave examined the effect of boundary conditions on the 
dynamics of dispersive prey-predator species in a linear 
habitat. Similar analysis has also been carried out by 
Gopalsamy []3]]} on the dynamics of two competing and dispersive 
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species in a linear habitat and found that reservoir type 
boundary conditions lead to stable co-existence only if a 
delicate relation on the boundary holds. Harada and 
Bukao l 5 J have also investigated the co-existence of two 
competing species over a linear habitat of finite length by 
considering the nonlinear effect. 

Keeping the above in view, in the following, we investigate 
the effects of dispersal on the linear stability of the equili- 
brium state of two competing species system with intraspecific 
interactions and functional response (Lin and Kahn M) in 
a circular habitat. 


8.2 BASIC EQUATIONS 

Consider the evolution of two competing and dispersing 
species with intraspecific interactions and functional response 
in a circular (0 <_ r < R) habitat. The evolution of the 
species in such a habitat is described by the following system 
of dynamical equations (Lin and Kahn £ 8^ , Skellam £ ICQ ) 


sN l 

at 


= I^a^a, -,R 






) + 


D i ft- 

"1 r ar 


3H- 


( 8 . 2 . 1 ) 


SN p a p1 R 1 

bT" = II 2^ a 2'" a 22 1I 2“ TZaSq 


) + D 


1 

2 r dr 


(r 


«8 

dr 


where K 1 (r,t), Ng(r,t) represent the population distributions 
of two species at time t ; Dg are the co-efficients of 
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dispersal. She inter act ioii co-efficients , ag , » a j_2 5 a 2 1 ’ 

— 1 

agg and a are positive constants. ihe term (1 + aN^)“ 
which arises due to rhe functional response in the model is 
characterised by a. 

'Ihe non trivial positive equilibrium state (II*, Ng) of 
the system (8.2.1) is obtained from 


a 1 pHp 

a 1 -a 11 N 1 * - = 0 

1 11 1 1+aN* 


a 91 N? 

a p -a pp Mj - -Z-L-± 
* ^ * 1+aN? 


(8.2.2J 


= 0 


It has been verified from (8.2.2) that there exists only one 
positive value for H* and N* under the following conditions 
(see Chapter VI) 



a 22 > 

fg > a 21 

(8.2.3) 


a 12 

a l a ll 

and 

a 21 > 

ag oc • 

(8.2.4 ) 


The system (8.2.1) is linearised to study the local 
stability of the equilibrium state (N*, N*) by writing 

' N 1 (r,t) = N* + u(r,t ) 


Ng(r,t) = N* + v(r,t) 


(8 ? 2.5) 


which on substituting in (8.2.1) and making use of (^8,2,2) 
gives 
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A L_ ( r 1Z) 
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( 8 . 2 . 6 ) 


where u,v are perturbed population distributions and 

HJ 

e — (a^-r-a^a + 2a., -a H?) 

(1+aNj) 11 1 11 1 

f = a 12 fl 

(1+ ah'*) 
a ?1 iTp 

s = uTS*)® ; h = as2l 2 (8- s - 7 > 

which are positive in view of (8.2*3) and (8.2.4). 

If the species are nondispersive, it has been shown that 
the nontrivial positive equilibrium state (N£, Rg) is 
asymptotically stable for a > 0 (see Chapter VI) provided 

eh > fg. (8.2.8) 

In the forthcoming sections, the stability of the equili- 
brium state (R*, R*) of the system (8.2,1) is analysed by 
associating the linearised system (8.2,6) with the following 
initial and boundary conditions. 

(i) Reservoir conditions 

u(r,0) = u Q (r ), v(r,0) = v Q (r ) 0 < r < R 

u(R,t) = VS , v(R,t) = V 1 t > 0 (8.2,9) 

u(r,t) and v(r,t) are finite at r=0 . 



169 


(ii) Efflux conditions 

u(r,G) = u Q (rJ, v(r,0) = v Q (r) 0 < r < R 

+ ku = Ug, ||-+ lev = V 2 at r = R (8.2.10) 

u(r,t) ana v(r,t) are finite at r=G . 


8.3 S1A3ILI1Y MaLYSIS Jim DISPERSAL IJHEER RESERVOIR 
COHDITIOUS 


In this case, on solving the system (8,2,6) with conditions 
(8.2,9) in a circular habitat (0£r<R) by using Laplace and 
finite Hanbel transforms, we obtain u(r,t) and v(r,t) as 


2 v- x n J o^ rx n^ ^ 


= I l + B^xpC-p^DldT 

f £ ^ 

0 11 2 0’in ex P^Pln i:) + F 2n ex P^2n t} I3 


+ ^2 2 


R* n [J i(Rx n ) 3 


(8.3.1) 


v 


0 J (rx ) ^t 

(r,t) = ^ E x n J 1 (R x n ) q E J n e3p ^" f> ln 1 ^ + B n exp ^" p 2n' I ^ <iT 


J ( J?X ) 

+ 2 _ 2 —2 o C G in exp (-Pm 17 ) + G 2 n exp ( -~£ , 2n' t ) 3 


n [Ji(Rx n ):] 


(8.3.2) 


where J m (x) is the Vessel function of order m and the 
summation is taken over all positive, roots, x n , of the 
transcendental equation 


J 0 (Rx) = 0. 


(8.3.3) 
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The constant coefficients used in (8.3.1) and (8.3.2) ■ 

are given by 


A = 

n 


■^(P ia~ i) 2 :x n" il ^ 7J l +f ' D 2 Y 1 
Pin" P 2 n ~~ 


B n = 


\ fr* n ~ p 2n l u l~ f ' -V 7 ! 
p ln“ p 2n 


C n = 


B 2^ p In-^j-^n" 6 [ { 
p ln“ p 2n 


E = 
n 


E 2 (DlX n+ e-pg n )Yi-gEiI J i 

p ln“ p 2n 


(8.3.4) 


E = 
In 


(Pi n -E 2 xJ-ii)u oI1 +xv 


on 


p In p 2n 


S' = 

2n 


(D g y^p 2n ) u on-fv 0 n 
Pin” p 2n 


G 


In” 


(p ln~ ;D l :K n- e ) v ~on 4 'g u on 
p ln“ p 2n 


G 2I r 


(D lV e-p an )Yo D -gu on 

Pin" p 2n 


(8.3.5) 



I < I) l 3 4 +D S x ^ +e+h > 


i I C (D^-DgX^+e-h)^ 4fgl] 1 / 2 


(8.3.6) 
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where u on , v on are defined as 

R 

v. = s ru 0 n) dr 

0 


v or: = ■> < - I) J 0 ( rx n j 41 


(8.3.7) 


?rom (8.3.6), it is noted that p ln , p 2n (n = 1,2,3,...) 
satisfy the equation 

p 2 -(D 1 x^+D 2 x^+e+h)p + (D 1 x 2 +e)(D 2 x^+h)-fg = 0 
and hence P ln » Pg n (n = 1,2,3,...) are positive if 


min (D 1 x 2 +e)(D 2 x 2 +h) > fg. 


(8.3.8) 


If (8.3.8) is satisfied i.e. P ln »Pg n (n = 1,2,3,...) are 
positive, u(r,t) and v(r,t) given hy (8.3.1) and (8.3.2) take 
the following forms respectively as t - 


P J (rx ) A B 

It u(r,t) = # I * -2—2- r— + — 

t — n J^(Hx n ) p m p 2n' 


p J n (rx ) C E 

It v(r,t) = | Z x --S— 2- r-2- + JL.~\ 

, ^ 9 J R n t f "D-rr- \ n . n _ — * 

{jr —*■ oo 


(8.3.9) 


(8.3.10) 


R n n J 1 (R^ n ) Up ln ' p 2n' 

It may be verified that the expressions (8.3.9) and (8.3.10) 
satisfy the steady state form of the system (8.2.6) with boundary 
conditions (8.2.9). fhus it is seen that the nonstationary 
solutions (8.3.1) and (8.3.2) tend to the stationary solutions 
(8.3.9) and (8.3.10) respectively as t - showing the 
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asymptotic stability of the equilibrium state provided (8.3.8 ) 
holds. 

Keeping in view trie inequalities (6.2.8 ) and (8.3.8), 
it is noted that the stable equilibrium state without dispersal 
remains so with dispersal. It may also be noted that an 
otherwise unstable equilibrium state without dispersal can 
become stable with dispersal under (8.3.8). 

Even if the inequality (8.3.8) does not hold, the 
equilibrium state may become stable provided = Dg and the 
initial and boundary conditions satisfy the following relations 


(h - Xg) JJ 1 = £V ± 
or (e - Xg) V 1 = gU ± 

(h - x £ ) u Q (r) = fv Q (r) 
or (e - x g ) v Q (r) = gu 0 (r) 


(8.3.11) 


(8.3.12) 


where X g i s given by (6.3.3). 

It should be here noted that the above conditions are not 
applicable to realistic situations. 


8 .4 SIABUM ANALYSIS WITH DISPERSAL UNDER EPFLUX COEDITIONS 

On solving the system (8.2.6) with conditions (8.2.10) 
in a circular habitat (0 < r < R), we get u(r,t) and v(r,t) as 
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u(r,t ) 


a 2 I n ( r ) / ex P("Pin 2 ^ + ex P dT 

Jl r\ 


o 


Ir,(r) 


+ 


jj2 £ J 0 (fty n 7 ex P(-Pm t - 1 + ? 2n exp (~ p 2n t - > 3 C 8 - 4 - 1 ) 


v(r,t) = | £ i n ( r ) / [_0^ exp(-p^) + 3^ exp(-p^i;3ai 


+ -- £ 31 


2 

n 

*nOO 


r 2 n J 0 (Ry n ) 


[Z ex P ( —P ^ + %n ex P ^”P2n"^ ^ -3 (8.4.2 ) 


where the summation is taken over all positive roots, y n , of 


the transcendental equation 

k J 0 (Ry) = yJqCB-y) 


(8.4.3) 


and the function I (r ) is defined as 


J ( ry ) y 
t ( r ) _ ^n 

n v J 




(8.4.4) 


The constant coefficients in (8.4.1) and (8.4.2) are given by 




, D l^ln - Vn “ h ) U 2 + f ¥2 


B. 


n 


C T = 
n 


1 * 

n 



Pin - 

p 2n 


■ 

■h h - J 

p 2iP U 2 

- fD 2 V 2 


^ ! 

p ln ~ 

P 2n 


D 2^ p ln ' 

2 

- D„V 
l J n 

- e)V 2 

+ gDy 2 


^ t 

p ln “ 

fan 


Wn 

+ e - 

P2n> 7 2 



Pin ~ 

p 2n 



(8.4.5) 
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P' 


in 


'2xi 


G ln 


G 2n 


^ln ' 

.. 2 
“ Vn 

- h) 

u on 

4- 

fv 1 
on 


P 

m “ j 

?2n 



2 

'■Vn 

+ h - 


u on 

- 

fv f 
on 


p 

In “ p 2n 



(Pin ’ 

- 

- e) 

v on 

+ 

§ u on 


P 

in “ p 2n 




+ e - 

p 2n^ 

V f 

on 

- 

§ u on 


(8.4.6) 


p ln ” p 2n 



I (D l y n +B 2 y n +e+1 ^ ±5 CtBxy^gy^e-h^fgjVS 

(8.4.7 ) 


where u’, v’ are defined as follows : 
on 7 oji 


R 

u on = f *%<■*> J 0 ( ry P dr 
o 

S 

v on = / ^0^) J 0 (^ n ) ar. 
0 


(8.4.8) 


As in the previous section, it is noted from (8.4.7) 
that p^, pg n (n = 1,2,3,...) are positive provided 


min (D^+e) (D 2 y^+h) > fg . (8.4.9) 

If (8.4.9) holds, i.e. pjj^, p^ n (n = 1,2,3,...) are 
positive, (8.4.1) and (8.4.2) converge to the following forms 
respectively as t -* 
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a! B* 

It u(r,t) = | Z L(r) + -~H (8.4.10; 

t - ~ n p ln g 2n 

Q f g * 

It v(r,t) = | ;yr) [iti-] . (8.4.11) 

* -* " n Pin P2n 

It may "be notel t.iat the expressions (8*4.10) and (8. -±,11) 
satisfy the steady state form of the system (8.2.6) with 
boundary conditions (8.2.10), showing the asymptotic stability 
of the equilibrium state. It is concluded from (8.2.8) and 
(8.4.9) that the equilibrium state which is stable without 
dispersal remains so with dispersal and also the effect of 
dispersal is to stabilize even an otnerwise unstable equili- 
brium state under (8.4.9). 

Even if (8.4.9) is not satisfied, the equilibrium 
state can become stable provided E^ = Eg, the initial conditions 
satisfy the relations (8.3.12) and the boundary conditions 
satisfy the following relations 


(h- V U 2 = f7g 
or (e - X 2 ) V 2 = gUg 


(8,4.12) 


where X^ j_s defined in (6.3.3) 
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8.5 COL CLUSIOI'T 

In circular habitats, effects of dispersal on the 
linear stability of the equilibrium stare of two competing 
species system ,vith intraspecific interactions and functional 
response have been investigated under reservoir and flux 
conditions. It has been asserted from (8.2.8;, (8.3.8) and 
(8.4.9) that the stable equilibrium state without dispersal 
remains so with dispersal. However, an otherwise unstable 
equilibrium state without dispersal can become stable with 
dispersal provided certain inequalities such as (8.3.8), 

(8.4.9) involving dispersal co-efficients are satisfied. Even 
if the conditions such as (8.3.8), (8.4.9) are not satisfied, 
an otherwise unstable equilibrium state may become stable 
under the following restrictions j 

(i) both the species should have equal dispersal ability 

i .e . = Dg . 

(ii) initial conditions should satisfy the relation (8.3,12). 

(iii) the boundary conditions should satisfy the relations 
such as (8.3.H), (8.4.12). 

Thus it is concluded that dispersal can stabilize 
even an otherwise unstable equilibrium state provided either 
certain inequalities involving dispersal co-efficients or 
certain relations involving initial and boundary conditions 
are satisfied. 
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CHAPTER IX 


SPRIG IS 01 DISPERSAL GL SIAEIuIiT OP IwO SPECIES 
SYSIELI PIIH I.UIUALISIIC ILiSRAC HCuS AXIL 
PULCTIGLAL RESPONSE 

9.1 IUIRODUC LICIT 

In Chapter IV, effects of convective and dispersive 
migration on the linear stability of the equilibrium state 
for two species system with mutual istic interactions and 
functional response have been discussed for constant convective 
and dispesive co-efficients under non- homogeneous boundary 
conditions. In this chapter, the effects of variable 
dispersal co— efficients on both the linear and nonlinear 
stability of the equilibrium state of the above mentioned 
system under homogeneous boundary conditions, using Liapunov's 
Direct Kethod , have been investigated. 

Ihe study related to the nonlinear stability of the 
equilibrium state for systems of interacting species are of 
recent origin, Alan Hastings Cob C 2 I3 > dome and 

Carmi []3^|> Segel and Levin Q 6 • Effects of dispersal 

on the nonlinear stability of the equilibrium state have also 
been investigated, Alan Hastings C i] > Jorne 811(1 Garnii c 33 , 
Segel and Levin [2^13 • Eventhough G-oh [[2 3 established 
global stability in a class of prey-predator models including 
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functional response (Kazarin off and Driessche C4j, Lin and 
Kahn [6]. Segel and Levin 6 3 ) » no attempt has been made 
to study such effects in the case of mutualistic interactions 
of species (‘Travis and Post £ 7 t Vandermeer and 3oucher 
dsn ). It should be noted nere that the effects of dispersive 
migration have also not been investigated on such interactions. 

Keeping the above in view, in the following, we investigate 
the effects of variable dispersal co-efficients on the linear 
and nonlinear stability of the equilibrium state for two 
species system with mutualistic interactions and functional 
response under homogeneous boundary conditions. 


9.2 BASIC EQUATIONS 

Consider the mutualistic interactions of two species in a 
finite (O^x^L, L being the length) one dimensional linear 
habitat. The equations governing their evolution with 
functional response and dispersal can be written as follows 
(Lin and Kahn d > Se gel and Levin d ^ 3 > Travis and Post d^d* 
Vandermeer and Boucher d 8 d ) : 


Q1L 

aT 


3N 2 

W 


a 1 pKp 3N1 

1I l^ a l" a ll K l + T+aN^ + dx ^ D 1 dx ^ 


ag^i 


m e 


( 9 . 2 . 1 ) 


V a 2~ a 22V * fe (D 2 d^> 


where N 1 (x,t), Ng(x,t) represent the population distributions 
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of two species at time tj D^Dg their variable co-efficients 
of dispersal and are positive. The interaction co-efficients 
a 1 ,ap,a 11 ,a 1 g,a 91 ,a gg and a are positive constants. The 
term (l+all^)“ corresponds tc the functional response in the 
model and the constant a may be interpreted as the measure 
of the strength of this response (Kazarinoff and Driessche 

The nontrivial positive equilibrium state (1»£, 1 i*j of 
the system (9.2,1) is given by 


ci-j 4- 


a 12 M 2 

1+cdSj 


a ll^l 


a 2 + 


a 21^1 

1+ocN* 


a N* 

a 22 2 • 


( 9 . 2 . 2 ) 


It can be easily verified from (9.2.2) that there exists only 
one positive value for N* and H* (see Chapter IV). 

Writing 

N 1 (x,t) = N* + n 1 (x,t) > 0 

(9*2*5 

H g (x,t ) = N* + n 2 (x,t) >0 

in the system (9.2.1) and using (9.2.2), we get the following 
system of nonlinear equations 
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an. 

FT 


(llj+n^ _ 

— g S+a^gaNg )n^ + a^g ( 1+ah* ) n g3 


- 6ll 1 

+ ~ (I), T~) 


an 2 

FT 


(ii*+n 9 ) 


? an 

1“ a 21 n l~ a 22 Sll 2 -3 + ” ^ 


ax ^ i ax 
an, 

ax ( J 2 W 


(9.2.4) 


wnere 


S = ( l+ah|' ) ( l+aEj*+an^ ) . 


(9.2.5) 


The linearised version of the system (9.2.4) is 

3 n i 3n 1 

W = - en 1 + f 2 ig+ q^(D 1 2 ^- ) 


8n . an 

gn-j-hn^ ~ (3, 


at 


2 ax v 2 ax 


) 


(9.2.6 ) 


where e,f,g,h are positive constants given by 

^ aa l2^2 _ a l?^l 

e = N* [a 1+ — - * P D , f = 12 1 


1 U11 (l + aht) 2 ' 


(1+aNJ) 


g = 


a 21^2 


( 1+ctN*) 


*\2 ’ h a 22 N 2 


(9.2.7 ) 


The systems (9.2.4) and (9.2.6) are associated with the 
following homogeneous boundary conditions: 

(i) Flux conditions 


an 1 (o,t) 

ax 


an^ (I j t ) 


3x 


= 0 


an (0,t) an 9 (L , t) 


(9.2.8a) 


3x 


ax 


o 
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These conditions imply that there is no migration of the 
species across the boundary of the habitat. 

(ii) Reservoir conditions 


^(Gjt) 


= n 1 (L ,t ) = 0 


n 2 ( 0 j t ) = n<p(L , t ) = 0 . 


(9.2.8b) 


These conditions suggest that the species populations remain 
at equilibrium level on the boundary of the habitat. 

In the following, we study the linear and nonlinear 
stability of the equilibrium state (If*, N*) of the system 
(9.2.1) using Liapunov’s Direct Method. 


9.3 LINEAR STABILITY ANALYSIS 

To study the local stability of the equilibrium state in 
absence of dispersal, the following positive definite function 
is considered 

E(n 1 ,n 2 ) = | (n^ + cn^) (9.3.1) 

where c is a positive constant to be chosen suitably. The 
time derivative of (9.3.1), which on using (9.2*6), becomes 

|| = -en^ + (f+cg) n ± xi 2 - chn^ . (9.3.2) 

Since eh > fg, by choosing 


c = f/g 


(9.3.3) 
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it may be noted from (9.3.2) that < 0 for n^ ^ 0, n 2 ^ 0 
and dE/dt = 0 only when n^= 0 and n^ = 0 ensuring the 
asymptotic stability of the eq ailibrium state (h*, II*) . 

fo see the effects of dispersal on the linear stability 
of the equilibrium state, consider the positive definite 
function 

E(n 1 ,n 2 ) = | f (n| + cn 2 )dx 

where c is chosen as in (9.3.3). 

ihe time derivative of (9.3.4) is 

j-m L Qn 1 dn p 

dT = { (n i W + cn 2 dx 

which on substituting (9.2.6) gives 


(9.3.4) 


(9.3.5) 


£§=-“/ C ^^-(f+cg^^g+chngj dx 


a dn i 

"1 3x^1 dx" 


+ J Ol Av ) + cn c> A -v 


dn 

2 3x v "2 5x 


3 


(9.3.6) 

Integrating the second integral in (9.3.6) and making use of 
homogeneous boundary conditions (9.2.8), we get 

Xi 

H = - / Z en^-(f+cg)n 1 n 2 +chng ^ dx 
o 

L 3n 1 2 L an 9 2 

" { D l ( — • ) dx - C { D 2 ( — } dx * 


(9.3.7) 
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Since the first integral of (9.3.7) is nonpositive, it is 
seen that <_ 0 and ~ = 0 only when n^= 0 and n g = 0, 
showing tnat the equilibrium state is asymptotically stable 
with dispersal also. 

Thus it is concluded that the stable equilibrium state 
without dispersal remains stable with dispersal as well. 


9.4 H OKI 31 BAR STABILITY ANALYSIS 

To investigate the nonlinear stability of the equilibrium 
state in absence of dispersal, consider the positive definite 
function 


n t ftp 

-L \ . _ -to-tt -I «. / -I . C ' 


E(n 1 ,n 2 ) = n 1 -llj ln(l+ — ) + c[> 2 -N| ln(l+ -&) 2 


h; 


n 1 > -if, n p > - Nj 


(9.4.1) 


where c is a positive constant to be chosen appropriately. 
The time derivative of (9.4.1) is 


dE n l ^1 A cn 2 ^2 

anr + 




N i +n i 


H i +n 2 


(9.4.2) 


which on using (9.2.4) becomes 


d^ " S ( a ii s+a i2 ctN 2 ^ n l“ ca 22 n 2 
+ &l2^ +a ^l^ + ca 21-^ ”S 


(9.4.3) 
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By choosing 


c = 


a 12 

a 2l 


(1+ctu*) 


(9.4.4) 


it is noted from (9,4.3) that ~| < 0 and = 0 only when 
n^ = 0 and = C. hence the equilibrium state (!>!*, U*) 
is asymptotically stable in the entire positive quadrant 
P = { (1I 1 ,1T 2 )|H 1 > 0, l\ig > 0} of the p base -p lane . 

To examine the effects of dispersal on the nonlinear 
stability of the equilibrium state, the following positive 
definite function may be considered 

E(n 1 ,n ? ) = f f n-h* In (1+ — ) 

1 d o 1 1 N* 


n„ 

+ cCn p -3J* ln(l+ -£)3>dx (9.4.5) 

c * 


where c is chosen as in (9.4.4). 

The time derivative of (9.4.5) is 


# = / r-tn !f± + 2f|-i ta 

^ o L Hpn 1 ^ »|+n g VEJ x 

which on using (9.2.4) gives 


(9.4.6) 


dE _ r L / / _ „ . . n ! 


dJ 


/ {(a^S + a 12 aK 2 ) + 0 a 22 n 2 


n„n 


- [a 12 ( 1+aH* )+ ca 2qIl ~jr-^ ) dx 


+ / c 

o 


n. 


H^+n^ 


— (D 
dx ^ 


9n„ 


cn. 


1 9x ^ + m* 


1 2 +n 2 


a 9n c 

£_ fry £ 

9x ^2 9x 


^)3dx. 


(9.4.7) 
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Evaluating the second integral in (9.4*7), keeping 
the homogeneous boundary conditions (9.2.8), we get 

o 


dS 

dt 


/ {(a 


llS+ a i 2 cai*) 


3 

S 


+ c 


a 22 n S 


in view 


_ inn 

~ &12( + ca giJ ~ s~~ ~ ^ 


1 


- / 


o 


pi s r 

(H*+n 1 ) 2 


dn 1 

> 


2 L 

dx-c / 
o 


yg 

(N*+n g ) 2 


an p p 

<53T> ta ' 


(9.4.8) 


Since the first integral is nonpositive, it may be noted 
from (9.4.8) that ~ < 0 for n^ ^ 0, n g ^ 0; = 0 only 

when n^= 0 and n g = 0. Hence the equilibrium state 

(H*,1I*) is asymptotically stable in the entire positive 
quadrant P of the phase-plane. 


Thus it is concluded that the nonline a^ly asymptotically 
stable equilibrium state without dispersal remains so with 
dispersal. 


9.5 CONCLUSION 


Using Liapunov’s Direct Method, effects of variable 
dispersal co-efficients on the linear and nonlinear stability 
of the equilibrium state for system of two interacting species 
with mutualistic interactions and functional response have been 
investigated. It has been shown that the equilibrium state is 
asymptotically stable in the entire positive quadrant of the 
phase-plane with or without dispersal. 
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10.1 IL-RCli'CrlcK 

In Chapter V, the linear stability of the equilibrium 
state for prey-predator system with functional response has 
been studied by considering the effects of convective and 
dispersive migration and solving the linearised equations 
under non- homogeneous boundary conditions. In the present 
chapter, we study the linear and nonlinear stability of the 
above system under homogeneous boundary conditions, using 
Liapunov's Direct Method, as in the case of mutualistic 
interactions studied in Chapter IX. Ihe domain of attraction 
of the equilibrium state is also estimated. 

In recent years, a great deal of attention has been 
devoted to the study of nonlinear stability of systems of 
interacting species and sufficient conditions have been 
established for the nonlinear stability of the equilibrium 
state by considering temporal variations of the species, 

Alan Hastings £ 1”| , Gatto and Rinaldi £ 53 , Goh £ 6,8,93 > 
Goh and Agnew £ 7 3 , Hsu £ H3 . In particular, Gatto and 
Rinaldi £ 53 ; considered a generalized prey-predator Lotka- 
Yolterra model and studied the stability properties of its 
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non- trivial equilibrium state by means of an energy function 
first proposed by Volterra in the context of conservative 
ecosystems. Hsu j_llj presented certain criteria for 
nonlinear stability of the locally stable equilibrium state. 
Effects of dispersal on the nonlinear stability of the 
equilibrium state have also been investigated, alan Hastings 
J orne and 0 armi C12J . Hosen £ 15j . In particular, 

Jorne and Carmi [] 12 3 considered the diffusive Lotka-Vol terra 
system of equations with a number of interacting species and 
concluded that the diffusive system with zero flux boundary 
conditions is stable for different positive dispersal co-efficients 
and the role of diffusion is to damp out all spatial variations. 
Event ho ugh G-oh 9 H has established nonlinear stability in 
prey-predator system with functional response (Hoi ling £ ICQ , 

Lin and Kahn £ 14 J ) , the effects of dispersal on the nonlinear 
stability of such systems have not been investigated. 

Keeping these in view, in the following, we examine the 
effects of variable dispersal co-efficients on the linear and 
nonlinear stability of prey-predator system with functional 
response. 

10.2 basic equations 

Consider the interactions of two dispersive species 
constituting prey-predator system in a linear one dimensj(onal 
finite habitat (0£ x < L, L being the length) with intraspecific 
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interactions of prey and functional response. The dynamical 
equations governing their evolution with dispersal along 
x-direction can be written as follows (Kazarinoff and Driessche 
[3 13]] » H n and Kahn 3 14 J , Jorne ana Card £ 12 J , Segel 
and bevin £ 16 3 ) 


5l 1 

W 

as, 

at" 


& A 

, lV a f a irr ICT 


i.’l ^ CL *"" , cL -j < — i . 


12 “ 2 ) + 8_ (3 flj 
3x 1 3x 


Qlv 


h 2 (-a 2 + ) + g x (r> 2 ax^ 


( 10 . 2 . 1 ) 


where II^(x,t), Ng(x,t) are population distributions of prey 
and predator species, respectively, at time tj D^jDg their 
positive variable dispersal co-efficients. The interaction 

co-efficients c^l ? 3-g J a ll ,a l2* a 2l 3X1(1 a 3X8 P os i^i Te constants 
The functional response in the model is characterized by the 
term (l+(xW 1 )“ where a measures the strength of the response 
(Kazarinoff and Driessche £13 3)* 

The nontrivial positive equilibrium state (II*, N*) of 
the system (10.2.1) is given by 


a l“ a ll N l 


a 12 N 2 

1+aN? 


= 0 


3- 0 


a 2l^l 

1+aN* 


= 0 


( 10 . 2 . 2 ) 
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which on solving gives 


a r 


i: 


a 2l~ a 2 a 


i]* 

a 2 




a 


z± 

12 


a 2 a ll ■ 


( a 21 a 2°^ 


2 


(10.2.3) 


It is noted from (10.2.3) that both E* and ITg are positive 
provided 


a l^ a 21“* a 2 a ^ ^ a 2 a H * 


(10.2.4) 


To investigate the stability of the equilibrium state 


(N£,N|), writing 


N 1 (x,t) = E* + n 1 (x,t) > 0 
N 2 (x,t) = N* + n 2 (x,t) > 0 


(10.2.5) 


in the system (10.2.1) and keeping in view (10.2.2), we get 
the following system of non-linear equations 


dn 1 (N*+n 1 ) 

W = 


an. 


an 0 (N*+n 2 ) 


1 1 y a 

S C (a 12 aN 2 -a 11 S)n 1 -a 12 (l+alJ)n 2 (]+ j£~) 

(10.2.6) 


2 

FV 


where 


a 21 n l + 8x 


an, 

2 3x 


■) 


S = ( l+oN*)( l+ocN*+an^) 


(10.2.7) 


The linear system corresponding to (10.2.6) is 
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where 


e = 


a 6n, 

fn Q + $— (D 1 -r— = 
2 ox v 1 3x 

0n o 

2_ ( o 2 x 

1 T ox lJ 2 3x ' 


[][a^^ ( Ug j_+ agoc ) ™ a^a ( ag ^"~UgOC ) 


0n 1 


8V = 

-en^ - 

8n 2 

gr. ± + 

aTT - 


( 10 . 2 , 8 ) 


f = 


(1+aKf) 


a i2 ji r 


( ag '^"■“ci.gCX ) 


(10.2.9) 


g = 


a 21 H 2 


( 1+ccN*) ' (1+cxNj)' 

and n (x,t), n p (x,t) are small perturbed population 


*1 v ~ y y 2 
distributions. 


From (10.2.9), it may be noted that f,g are positive and 
e is positive for a ±± > a Q > 0 and negative for 0 < a 1;L < a Q , 


where 


a„ = a., a 


( Ug -j-aga ) 


1 ( ag aigOC^ ) “ 


> 0 for a > 0 . 


(10.2.10) 


As considered in the previous chapter, the following boundary 
conditions ahe associated with the systems (10.2.6) and (10.2.8); 

(i) Flux conditions 

an 1 (0,t) 3n 1 (L,t) 


5x 


9x 


= 0 


an (0,t) _ 9n 2^ L>t ^ 

— dx 


0 . 


(10.2.11a) 
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The boundary conditions (lu.2. 11a) represent the fact that 
the population is confined within the region and there is 
no migration across the boundary. 

( ii ) Reservoir conditions 


n 1 (0,t) = n 1 (l,tj = 0 
n g( 0 , t) = n g (L,t) = 0 . 


(10.2.11b) 


These conditions imply that the species populations at the 
boundary are eq.ua! to their equilibrium Values. 

In the forthcoming sections, using Liapunov's Direct 
Met nod, we investigate the stability of the equilibrium state 
(HJ, N*) under (10.2.4). 

10.3 LINEAR STABILITY ANALYSIS 


To examine the linear stability of the equilibrium state 
in absence of dispersal, consider the positive definite function 

B(n 1 ,n 2 ) = | (n^ + cn g ) (10.3.1) 


where c is a positive constant to be chosen appropriately. 
The time derivative of (10,3.1) is 



dn 1 

n l cHT + 


dn 

cn 2 dT~ 


(10.3.2) 


which on using (10.2.8) gives 

= -en^ + (f-cg) xi ± xx 2 . 


(10.3.3) 
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Choosing 


f/i 


( 10.3.3; becomes 


dB 


dJ = " en 


(10.3.4) 


(10.3.5) 


Brom (10.3.5), It is noted that ~ < u for e > C 
and n^ ^ 0. Defining the set R = { (n^,n 2 )| ^ = CJ , 

It can ha easily verified that the origin is the only 
invariant subset of R and, therefore, on using the theorem 
on extent of asymptotic stability (see CO. p 224 ) , it 
may be noted that the equilibrium state is asymptotically 
stable provided e > 0. 

Jo see the effects of dispersal on the linear stability 
of the equilibrium state, consider the positive definite 
function 


E(n x ,n 2 ) = | f (nf + cn 2 )dx 

where c is chosen as in (10.3.4). 

ihe time derivative of (10.3.6) is 


(10.3.6) 


dE 1 an l dn 2 

d^ = { (n iFT + cn 2 a^ dx • 


(10.3.7 ) 

Using (10.2.8), (10.3.4) in (10,3.7) and rearranging, we get 


i - - { L 4°* + { k i ^ + k ^ 


(10.3.8) 
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Evaluating the second integral and making use 'of the homo- 
geneous boundary conditions (10.2.11), we have 


1 © 1 on 1 2 3 n Q 

dB _ „ r a.. r t- -rs / 1 \ . _-n / 2 


dT 


-e / n^ dx - / r D i( + cIi 2^3x ^ 3 dx ’ (10.3.9) 


It is seen from (10.3.9) that gj < 0 for e > Oj = 0 
only when n^= 0 and n^= 0. Hence e > 0 is the sufficient 
condition for asymptotic stability of the equilibrium state. 
Thus it is concluded that an otherwise stable equilibrium 
state without dispersal remains so with dispersal. 


Even if e < 0, we can find certain condition for 
asymptotic stability of the equilibrium state under homogeneous 
reservoir conditions (10.2.11b) provided and Dg are 

positive constants. In such a case, using Poincare’s 
inequality (see £ 4 , p 225) in (10.3.9), we get 

-,-p D 1 n; 2 L « D 9 % 2 1 P 

H < -(e + S njdx-c / n~ dx . (10.3.10) 

L 0 L o 

It may be observed from (10.3.10) that ^ <_ 0 for 

D-jir 2 

+ e > 0 (10.3.11) 

n 

d.E 

and ^ = 0 only when n^ = 0 and ng = 0. 

Hence (10,3.11) gives the sufficient condition for asymptotic 
stability of the equilibrium state showing that dispersal has 
stabilizing effect. 



196 


10.4 MLIiiEAR STA3ILIH ANALYSIS 


Consider the following positive definite function to 
study the nonlinear stability cf the equilibrium state in 
absence of dispersal 


n i n 9 _ 

3(n.,n 9 ) = xi.-J* ln(l+ -±) + c[n 0 -17* Oja(l+ — )~\ 

1 * 1 1 N* ^ * Ng - 


n l > ~ i ' J i» n 2 > ” ^2 


(10.4.1) 


where c is a positive constant to be chosen suitably. 

On substituting (10.2,6) in the time derivative of 
(10.4.1), we have 


n: 


a xl a (l+aNp(I 0 -N 1 ) 


+ C? a g^-a^2(l +a ®*) 


(10.4.2) 


where 


IT = 

o 


a^a(ag^-aga) - a 1 - ) a 


H 21 


a ll a ( a 21-V> 


(10.4.5) 


By choosing 


a 


c = 


a, 


12 

21 


(l+ocN*) 


(10.4.4) 


the following results are noted from (10.4.2) regarding the 
regions of stability: 
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(i) If 


1 < 0 
o — 


i .e , 


a 


11 


a^a 

- a^7 Ca 21 -a 2 a) 


dE 


(10.4.5) 


it may be observed that <_ 0 in the positive quadrant 
P = { (3^,IJg)ji^> 0, 1 j 2 > 0} of the original phase-plane. As 
the origin is the only invariant subset of R = {(n^ngjj ~| = 0} , 
the equilibrium state (21*, U|) is asymp totically stable in 
the entire positive quadrant P under (10.4.5). 


It is also remarked from (10,4.5) that it may not hold 
for a particular set of parameters value constrained by 
(10.2.4) showing the destabilizing nature of a. 


(ii) If 21 > 0 

o 


l .e 


a 


11 


a^a 

< ( a 81- a 2 a) 


(10.4.6) 


we 


can find a subregion A = { (N^,Ng) |N^> H q > 0, Ng> 0^ 


of the positive quadrant containing the equilibrium state, 
dE 


whsTs Tpjr 0 ? provided 

v (agl-aga) 

a ll > a l a , s 

( u.g j|4" a 2 a ) 


> 0 


(10.4.7) 


which is also the condition for linear stability i.e. e > 0. 
Combining (10.4.6) and (10.4.7), the criteria for asymptotic 
stability of the equilibrium state in A is 


anx (a ?1 -a ? a) 

— — (ag^-aga) > a^-j_> a^a * (10.4.8) 

21 (ag^+agcc) 
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Thus it is concluded that if the equilibrium state is linearly 
stable, then there exists a subregion A in the phase-plane 
where it is nonlineariy asymptotically stable under (10.4.8). 

further, it is seen from (10.4.8) that it may not valid 
for a particular set of parameters value, though the inequality 
(10.2.4) is consistent, showing that the functional response 
has destabilizing effect. 

The effects of dispersal on the nonlinear stability of 
the equilibrium state (N*,N*) can be studied by considering 
the positive definite function 


1 n. 

B(n 1 ,n ) = / (n.,-1* In (1+ -i) 

r s o - 1 x 

_ B o 

+ cfn -Nj ln( 1+ — — ) ”1 } dx 
^ c w* 

■“2 

where c is chosen as in (10.4.4). 


(10.4.9) 


The time derivative of (10.4.9) is 
n„ 9 n 


dB 

d^ 


L "**^1 

o ^ N J+n 1 9 1 " N|+n g d t 


1 Cn 2 Sn 2 n , 

+ 5T- J bX . 


(10.4.10) 


Substituting (10.2,6), (10.4.4) in (10,4.10) and rearranging, 
we get 

L n 2 

|| = / a 11 a ( 1 + aBj) (K o «H 1 ) gi dx 


n. 


dn. 


- f l ~— fe <J>1 * o 

o (llj+n 1 ) 9x 1 dx 


n , dn c 

2 - - - A— ( t ) — , ) } 

dx ^2 dx ; ax * 


(H*+n g ) 


(10.4.11) 
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Evaluating the second integral of (10.4.11) and keeping in 
view the homogeneous boundary conditions (10.2.11), we have 


dE 

dt 


L 


= J a. 1± a (1+aN*) (IT^- II 1 ) 


n 


1 


o - 1' S 


dx 


- J 


L 


dn 1 2 


1 Vi 


o (H*+n 1 ; 


q ^ ~ c J n 

dx O (K|+n 2 )^ 


an 9 2 
^“) dx * 


(10.4.12) 


The following possibilities arise, which are noted from 
(10.4.12), regarding the regions of stability; 

(i) If H <_ 0 i.e. (10.4.5) holds, the equilibrium state 
is asymptotically stable in the entire positive quadrant P 
of the phase-plane. 

(ii) If N > 0 and < N* i.e. (10.4.8) holds, the 
equilibrium state is asymptotically stable in the subregion 
A = { (H^Ug)! U.j> N o > 0, Kg > 0} of the positive quadrant. 

(iii) If K o > > 0, there exists a subregion 

B = {(I^,lTg)|o < < U o , Kg > 0} of the positive quadrant 

where the equilibrium state is asymptotically stable under 
reservoir conditions (for D^Dg positive constants) provided 
certain inequality involving dispersal coefficients is satisfied 
as shown below; 

In this case, the inequality (N* + n^) < N o is used 
in the second integral of (10.4.12) to get 



200 


§< - / a 13 i a(l + «r*) (%-K 0 >f 


ll 


n: 


dx 


>,H? L dn, 2 L IT* 

-cf J (q^) dx - CD 9 / ^ 

iJ o 0 


2 o (i«f-n 2 ) 2 Bx 


3n 0 2 

dx. 

> AV * 


(10.4.13 ) 


Using Poinc a re’s inequality (see CO. p 225) in the 
second integral of (10.4. 13), we get 

.2 


H < - / a 11 a(l+aU*)(S 1 -h o ) — dx 


n: 


D.tc 2 N* L o 

-L 1 r „<o 


/ nr dx-cD 0 / 


H 2* 


P J W-L/p J O Vi T 

L N^ o 1 ^ o (N|+n g r ° 


3n 2 
(t~) dx. 


(10.4.14) 


dB 


It is seen from (10.4.14) th a t ^ < 0 in the subregion B 
provided 


D^ 2 

p 2l a -i i a “7 
1 N* 


(10.4.15) 


ensuring the stabilizing effect of dispersal. 


10.5 G OH Cl US ION 

Effects of variable dispersal co-efficients on the 
linear and nonlinear stability of the equilibrium state for 
prey-predator system with functional response have been 
examined using Liapunov's Direct Method. It has been shown 
that the functional response has destabilizing effect on the 
equilibrium state . 


Ui 


The following conclusions have also been drawn: 

(i) In otherwise linearly or nonlinearly stable equilibrium 
state witnout dispersal remains so with dispersal. 

(ii) If the equilibrium state is linearly stable -witnout 
dispersal, then there exists a subregion of the positive 
quadrant in the phase-plane where it is nonlinearly stable 
with or without dispersal under (10.4.8). 

(iii) Dispersal has stabilizing effect on the linear and 
nonlinear stability of the equilibrium state as can be 
noted from conditions (10.3*11) and (10,4.15) respectively. 

(iv) Even if (10.4,8) is not satisfied, still we can find 
a subregion of the positive quadrant as the region of 
stability provided certain condition such as (10.4.15) 
involving dispersal co-efficients holds. 
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CHXPIEtt XI 


KB'?3013 03 DISP3RS-1 01. SI-3ILIIY CP IVC 0CI.IP31ILG 
SPI0I33 SYSIS.I rfTIH FULCUIQIul EESPOLSS 

11.1 il'zrgdugiic:: 

In Chapter VI, effects of convective and dispersive 
migration on the linear stability of the equilibrium state 
for two competing species system with functional response 
have been studied by solving the system of equations 
analytically under non- homogeneous boundary conditions. 

It has been assumed that the convective and dispersive 
co-efficients are positive constants. In this chapter, 
using Liapunov’s Direct Method, effects of variable dispersal 
co-efficients on the linear and nonlinear stability of the 
equilibrium state for two competing species system with 
functional response are studied under homogeneous boundary 
conditions . 

The nonlinear stability of systems of interacting 
species have been studied by many investigators with or 
without dispersal, Alan Hastings £l,2]j, Gatto and Einaldi 
[3 4 3 » Gilpin £ 5£ > Goh £6,8,9 £, Goh and itgnew [£ *7 £] , 

J orne and Carmi £ 11]] . Seramoto et al. £l5£ have 
demonstrated that the predation p res sure has a stabilizing 
effect on the community of competitive species by considering 
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a model of two preys and one predator having switching 
property of predation. Hall am []lO]]j has presented a 

model showing that a protoco-operative soecies can drive 
a stable competitive co ununity to extinction. It should 
be here noted that though Goh £ 9 j has studied the nonlinear 
stability of prey -predator system with functional response 
(Kazarinoff and Driessche £ 12 ]]] , Lin and Kahn £l3J), such 
investigations have not been carried out in the case of 
competing species system with or without dispersal. 

Keeping these in view, in the f ollowimg, we investigate 
the effects of -variable dispersal co-efficients on the linear 
and nonlinear stability of the equilibrium state for two 
competing species system with functional response under 
homogeneous boundary conditions. 

11.2 B*SIC EQUATIONS 

Consider the evolution of two dispersive species 
competing for common ecological niche in a linear one 
dimensional finite (0 <_ x< L, L being the length) habitat 
with intraspecific interactions and functional response. 

The equations governing their evolution with dispersal along 
x-direction can be written as follows (Kazarinoff and 
Driessche []]12]]J , Lin and Kahn £ 13 » Segel and Levin [] 14]]] ) : 
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oi; 


•i Q'-j P L ^ 0J.1 

^ ~r^) + fj g~0 


i v 1 u irr i+aiT 


( 11 . 2 . 1 ) 


fT" L 2^ a 2" a 22^2 ifel^ + ax ^2 9x * 


£ 


where Ij^Xjt), l!g(x,t) are population distributions of 
two species at time t; D^, Dg are variable dispersal 
co-efficients of the species and are assumed to be positive. 
The interaction co-efficients a i> a 2» a H ,a i2 ,a 2l ,a 22 3X1(1 a 

are positive constants. The strength of the functional 

— 1 

response (1+aE^)” in the system is determined by a 
(Kazarinoff and hriessche H 12 !])* 

The nontrivial positive equilibrium state (H*,N*) of 
the system (11.2.1) is obtained from 


* a 12 K 2 

a 1 -a 11 N* - = 0 

1 11 1 1+aHj 

a 01 L T f 

a ? -a ? H? - = 0 . 

^ * 1+aN* 


( 11 . 2 . 2 ) 


When a = 0, it should be here noted that the nontrivial 
positive equilibrium state of the system (11.2.1) is stable 
provided 


a 22 

a 12 



^21 ' 
a ll 


(11.2.3) 


In the case 


a > 0 , for N*, Kg to be positive, from 
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(11.2.2) and keeping in view (11,2.3), we have 


a r 


a 


- — > — > K* 

a 21 a 2 a a H 1 


(11.2.4) 


giving 


a 91 /* • 


(11.2.5) 


It is also noted from (11.2.2) that there exists only one 
positive value for II* and IT* under (11.2.3) and (11.2.5). 

In this chapter, therefore, we investigate the linear 
and nonlinear stability of the equilibrium state (N£,IT*) of 
the system (11.2,1) under (11.2.3) and (11.2.5) using 
Liapunov* s Direct Method . 


Using the transformations 

N 1 (x,t) = Nj + n 1 (x,t) > 0 

Ng(x,t) = N* + n 2 (x,t) > 0 

in the system (11.2.1) and keeping in view (11.2.2) 
the following system of nonlinear equations 


( 11 . 2 . 6 ) 
, we get 


an, 

FT 


(U*+n 1 ) 

g — - Q-(a 11 S-a 1 gOcHg)n 1 -a 1 g(l+aN*) n g 3 


a dn, 

+ ax ^i 


(11.2.7) 


an a p1 n 8n 

FT" = ( K 2 +]a 2^ S a 22 n 2^ + ax ^ D 2 8 x ^ 
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where 


S - ( 1+ouT^) ( l+ah'*+an^ ) . 


( 11 . 2.8 ) 


The linearised form of the system (11.2.7) is 


dn 1 dn 1 

rrrcr— = — en 1 -fn o + t— (i)., t — — ) 
ot 1 2 ox v 1 3x 

an p dn 

W = ~ gn l- to 2 + dx ^ 2 5#) 


(11.2.9) 


where n^, Hg represent the small perturbed population 
distributions and 


n: 


e = 


(1+aN*) 


(a^-a-joc + 2a 11 aH*) 


a 12 K l 

(1+ccHj) 


a 2l H 2 


(l+ccN*)‘ 


h = 


a 22^2 


( 11 . 2 . 10 ) 


which are positive in view of (11.2.3) and (11.2.5). 

The systems (11.2,7) and (11.2,9) are associated with 
the following boundary conditions: 


(i) Flux conditions 


Sn 1 (0,t) 

an^ 

(L,t) 

dx 

= dx 

= 0 

dn (0,t) 

dn 


dx 

— 

~ dx 

(lift) = 0 


(11.2.11a) 
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These conditions imply that there is no migration of tne 
species across the boundary of the habitat. 

(ii) Heserv'oir conditions 


n^(C,t) = n^(L,t) = 0 
n 2 (0,t) = n 2 (L,t) = 0 


(11.2.11b) 


These conditions suggest that the species populations remain 
at equilibrium level on the boundary of the habitat. 


11.3 LINEAR STABILITY ANALYSIS 

To study the linear stability of the equilibrium state 
(N*, N*) in absence of dispersal, consider the positive 
definite function 

S(n 15 n 2 ) = | (n| + cn^) (11.3.1) 

where c is a positive constant to be chosen suitably. 

Using (11.2.9), the time derivative of (11.3.1) can be 
written as 

1^1 = -en^-(f+cg)n 1 n 2 - chn| . (11.3.2) 

By choosing 

c = f/g 

it may be noted from (11.3.2) that ~| < 0 provided 


(11.3.3) 
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eli > fg 


(11.3.4) 


1 .e. 


a g2&ii“ a ( a i“ a ii®*) +a n a ^*I3 ( > a i2 a 21 (11*3*5; 


and ^ = 0 onl;' when n^=G and n g =G. 

Hence (11.3.5) gives the sufficient condition for asymptotic 
stability of the equilibrium state which is satisfied, in 
particular, if 


N* > 


a- 


2a 


11 


(11.3.6) 


It may be here remarked that the inequality (11.3.5) for 
a > 0 may not be satisfied, in general, for a particular set 
of parameters value, eventhough the inequalities (11.2.3) and 
(11.2.5) are consistent showing the destablizing nature of a. 


The effect of dispersal on the linear stability of the 
equilibrium state is investigated by considering the following 
positive definite function 

B(n 1 ,n g ) = \ S (n| + cn|)dx (U.3.?) 


where c is chosen as in (11.3.3). 

The time derivative of (11.3.7) is 
dE r L an l an 2 

aT = { (“i FT + on 2 TT* 6 * 


(11.3.8) 


which on using (11.2.9) becomes 
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dE 

d^ = 


1 

- J 

o 

L 

+ J 

o 


£ en 2 + ch n| + (f+cg) n^Ug^dx 

__ a fcn 1 3 8n„ 

LTl l 3x ^1 3x“^ + ^2 ex ^2 


dx . (11.3.9) 


Evaluating the second integral in (11.3.9) and keeping 
in view The homogeneous boundary conditions (11.2,11), 
we get 

H = - / [en^ + chn| + (f+cg) n^g^dx 

L dn 1 2 1 an 9 2 

- / ^(g—) dx - c / f> 2 ^dx ^ dx * (11.3.10) 

If the inequality (11.3.5) holds, it may he noted from 
(11.3.10) that < 0 for n 1 ^ 0, n g ^ 0; ^ = 0 only 

when n^= 0 and ng= 0. Hence the asymptotically stable 
equilibrium state without dispersal remains so with dispersal. 

Even if (11.3.5) does not hold, certain condition 
involving dispersal co-efficients can be found for asymptotic 
stability of the equilibrium state under homogeneous reservoir 
conditions (11.2.11b) when Dg are positive constants. 

In such a case, Poincare's inequality (see £ 3^ » P 225) may 
be used in (11.3.10) to get 

< - / [en 2 + (f+cg) n-jUg + chng^dx 

D.% 2 L „ cD p 7t 2 L p 

- -nr J ' n i - —2" S dx . 

L o Jj o 


(11.3.11 ) 
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By choosing 


a 


c = 


12 


l 21 


( 1+ocliJ) 


(11.4.1) 


it is observed from (11.4.3) that < 0 in the entire 
positive quadrant P = { (K^,hg ) jli ^ > u, ii_, > of the 
p hase-pl&ne provided 

a.gg ( l+alj£) &qq“ a ( ) H 2. a i2 a 9l (11.4.5) 


giving the sufficient condition for asymptotic stability of 
the equilibrium state (If*, N*) in the entire positive 
quadrant P . 

a i _ 

Since — > N* [see (11.2.4)2]} it is observed that the 
inequality (11.4.5) may not be satisfied for a particular set 
of parameters value constrained by (11.2.3) and (11.2.5) 
showing again the destabilizing nature of a. 

Even if the inequality (11.4.5) does not hold, there 
exists a subregion A = {(H^,Ng)| 0 < < H^, 0 < Ng > of 

the positive quadrant in the original phase-plane, in which 
< 0, containing the equilibrium state provided 


N* > If > 0 
1 o 


(11.4.6) 


where 


n = ■= 
o a 


a 12 a 2l 


a(a 1 -a 11 Ip 


^ _j_ . -*-r 

a 22 a n( 1 +aNp a ll 


- 13 (11.4.7) 
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ensuring the asymptotic stability of the equilibrium state 
(H*,17*) in k under (ll.4t.6j. 

To investigate tne effects of dispersal on the nonlinear 
stability of the equilibrium state (17*, II*), the following 
positive definite function may be considered 

E(n,,n ) = f { n -17* ln(l+ -;) 
x^o 1 x IT* 

_ ; n 

+ c [no-H* ln( 1+ — )]]} dx (11.4.8) 

where c is chosen as in (11.4.4). 

She time derivative of (11.4.8) is 

dE r 1 t- *1 dni i Cn2 

^ = { + ^r 2 

which on using (11.2.7) gives 
L 

■p = - J { (a 11 S-a l2 a!7g) + c 

_ n- n 

+ [a^g(l+aH*) + ca-giH ' g" ~} dx 
L 

+ / 
o 

Evaluating the 
of homogeneous 


n. 


H* + n 1 


(B- 


an. 


■) + 


cn. 


* 9n, 

i_ (j) t 

3x ^ 2 dx 


03 dx . 


9x v i ax ' N * 

2 + 2 

(11.4.10) 

second integral in (11.4.10) and making use 
boundary conditions (11.2.11), we get 


g^^dx (11.4.9) 
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dB 

dt 


L 


n 2 

J { (a 11 S-a 1 -aIJ|j + ca 99 nf 
o 


22 2 


- * . n l n P 

+ j_a 1 o ( 1+aL - ) + ca p1 j — • — } dx 


L D 


12 

T * on. 2 


11 


J 

o (il*+n 1 ) 


\2 ^dx 


- ~ L D il* on p 2 

CstT - ) dx - c J - - -S (tx— ) dx . 


-r* \2 '"ox 

o (^ 2 +n.g J 


(11.4.11) 


From (11.4.H), the following possibilities regarding the 
nonlinear stability can be asserted: 


(i) If (11.4.5) holds i.e. N q < 0, the equilibrium state 
which is asymptotically stable in the entire positive quadrant 
P of the phase-plane without dispersal remains so with 

disp er sal. 

(ii) If N* > N q > 0 i.e. (11.4.6) is satisfied, the 
equilibrium state is asymptotically stable in the subregion 
A of the positive quadrant even with dispersal. 

(iii) If H > K* > 0, we can find a subregion 

B = { (N^jhg) | 0 < < H q , 0 < K g } of the positive quadrant 

in which the equilibrium state is asymptotically stable under 
reservoir conditions (for D^,D^ positive constants) provided 
certain inequality involving dispersal co-efficients holds as 
shown below: 

In this case, using the inequality (N*+n 1 ) < H q in the 
second integral of (11 .4. 11), we have 
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dE . 
dt - 


2 

n -, 

* \ 1 


- J {(a 11 S-a 1 g <xH 2 ^ S + Ca 22 n 2 


- _ _- n l n o 

+ l_ai 2 (l+alif) + ca 2iJ“ 


D-IT^ L 6n_ 2 
- / (~+) dx-cDgiig J 


T 2 
"‘o 


O ( ii g+Hg ) 


dn 2 

dx 


•,2 v dx 


(11.4.12) 


Eurtner, Poincare’s inequality (see [ 3 [] , p 225) is used in the 
second integral of (11.4.12) to get 

L n 2 

dE x r r/_ 1 


^ / {(a 11 S-a 12 aNj) 3- + caggng 

_ n.n ? 

+ j_a l2 ( 1+aN*) + ca 2 lH " s ' " ^ 


D^ 2 Mf L 


1 ' U , 9 ^ > 

— p— — ^ J n-dx - cD ? Up J 
L o ^ ^ o (Hg+n 2 ) 


*. „ >2 v ox 


an„ 2 

to. 


(11.4.13) 


✓a-p 

Erom (11.4.13), it may he noted that < 0 in the 
subregion B provided 


B 1 tc 2 IT* 

ag 2 2 " 2 ) ( 1+ocN-j^ ) — a^gqEg ~~j ^ a 12 a 2l ( 11.4.14 ) 

L E\ 
o 


showing that dispersal has stabilizing effect. 

11.5 CONCLUSION 

Using Liapunov's Direct Method, effects of dispersal on 
the linear and nonlinear stability of the equilibrium state 
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for system of two competing species with functional response 
have been investigated. It has been noted that the functional 
response has destabilising affect on the equilibrium state. 

The following conclusions are also drawn: 

(i) An otherwise linearly or nonlinearly stable equilibrium 
stare without dispersal remains stable with dispersal as well. 

(ii) There exists a subregion of tie positive quadrant in 
the phase-plane where the equilibrium state is nonlinearly 
asymptotically stable under (11.4.6). 

(iii) Dispersal has stabilizing effect on the linear and 
nonlinear stability of the equilibrium state of the system 
as can be seen from the conditions (H.3.12) and (11.4.14) 
respectively . 

(iv) Even if (11.4.6) does not hold, still we can find the 
region of asymptotic stability which is also the subregion 

of the positive quadrant provided condition such as (11.4.14) 
involving dispersal co-efficients is satisfied. 
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3YY30 IS GI DISPSiSAL Cl. SEA3HI2Y C? JCi.CBI.HEA 
J3?ID3.;iC I-IGDEL 

1 2 . 1 Ii. x Rv. C uo ixu-> 

In Chapters IV to XI, effects of convective and dispersive 
migration on the systems of interacting species have been studied 
by considering functional response. It may be noted that 
migration can also affect the spread of infectious diseases 
under various environmental and geographical conditions. One 
such disease is gonorrhea caused by sexual contacts between 
males and females. It can be spread from one region to another 
by migration of either of these species and therefore this 
aspect is investigated in this chapter. 

The spread of infectious diseases has been extensively 
studied by various investigators considering temporal variations, 
Aronsson and Mellander []l”[» Bailey C 2 > 3 > 4 II» Gooic e and 
Yorke 3 7 3 , Cooke £ 8^ » Hethccte £ 9,10,113 » lajmanovich and 
Yorke 3 12 3 > Waltman 3 153 » Wichmann 3 16 3 • In particular, 
a mathematical model for gonorrhea has been first developed by 

73* Hethcote 

model to study the control of gonorrhea. Lajmanovich and Yorke 
3123 have presented a generalized model for gonorrhea in a 


Cooke and Yorke 3 


3113 used a host-vector 
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non homogeneous population and studied asymptotic properties of 
its equilibrium. Yorle et al. L M 3 have pointed out that 
saturation in a sexually active core population limits the 
incidence of gonorrhea and this core causes gonorrhea to remain 
endemic. 

Attempts have also been made to study geograpnical spread 
of infectious diseases by considering dispersive migration of 
susceptibles and infectives, Bailey Q2_j, Capasso on, 

Capasso and Portunato £ 6 » Mar cat i and Pozio [^13^, Radcliffe 

|~ 14 3 • In particular, Capasso and Portunato Q6 J studied 
the asymptotic behavior of a class of reaction diffusion 
equations related to models for the spatial spread of infectious 
diseases. Recently, Harcati and Pozio Z3 investigated the 

global behaviour of a vector disease model by considering spatial 
spread and hereditary effects which may be applicable to growth 
and spread of malaria. a theoretical expression giving the 
velocity of propagation for geographical spread of host-vector 
and carrier borne epidemics has been developed by Radcliffe £143 • 

Keeping the above in view, in the following, effects of 
dispersal on the stability of the endemic equilibrium state of 
the system governing the spread of gonorrhea are investigated 
using host-vector model. 
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12.2 BASIC EQUATIONS 

Consider that the total number of sexually active and 
promiscuous portion of the population, living in a linear one 
dimensional finite habitat (0 < x < 1) of length L, consists 
of two groups viz. suscepti'oles and infectives. Let c^,Cg 
be the total number of promiscuous males and females respectively 
and N 1 (x,t), Hg(x,t) be the total number of the corresponding 
infectives at the point x -and time t in the habitat. She 
dynamical equations governing the population distribution of the 
infectives are given by (Bailey [2,3,4], Capasso [[ 5[] , Capasso 
and F or tun at o [[ 6 [] , Mar cat i and Pozio [[ 13 [] , Radcliffe [[ 14 [] ) 

as. 8 2 h h a 2 !! 

FV = -ap 1+ b 1 (o 1 -H 1 )(lI 2 +a 1 — j-) + % ^5“ 

( 12 . 2 . 1 ) ' 

a 2 ®-, a 2 N 2 

FT = -a 2 H 2 +b 2 (c 2 -N 2 )(N 1+ a 2 ^-) + D g 

where D-^JDg are dispersal co-efficients of the infectives due 
to self -diffusion and assumed to be positive constants. The 
constants a^, a g are cross dispersal co-efficients of the 
inf ectives, a^ag denote the rate at which male, female 
infectives are cured and b lf b g represent the rate of infection 
of male, female susceptibles respectively. 

The endemic equilibrium state (N£, Eg) of the system 

(12.2.1) is given by 


I 
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,j„ = b l b 2 C l C 2- a l a 2 
1 a l b 2 +b l b 2 C 2 


b, b 0 c., c n -a, a 0 
n* = 12 12 12 

i4 2 agbi+b'ibg c " 


( 12 . 2 . 2 ) 


which are positive provided 

b^bpC^Cg > a^ag . (12.2.3) 

io investigate the stability of the equilibrium state 
(W*,!!*) of the system (12.2.1), using the following 
t r an sf ormar i ons 


U^(x,t) = h*+n^(x,t) > u 
K g (x,t) = N*+n g (x,t) > 0 
the system (12.2.1) can be written in the form 


(12.2.4) 


dn^ 

’at” 


an, 


6 2 n r 


•( a i +b i H 2 +b l n 2 ) n i +fn 2 +l) i q ^2 + fcc l dx 2 

S 2 n 0 


3n f 

£ 

at 


b l a l n l . 2 
3x 


3 2 n 2 d 2 n x 

-( a 2 +b 2 hJ+b 2 ni)n 2 +gni+D 2 + ga g — -g- 


3 2 n, 




2^ 2 2 s 

9x 


( 12 . 2 . 5 ) 


where 


f = b 1 (c 1 -N*), g = b 2 (c 2 -hj) . 


( 12 . 2 . 6 ) 


It may be noted from (12.2.2) and (12.2.6) that f and 
g are positive constants as c^ > and c 2 > N*. 
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The linearised version of the system (12.2.5) is 


where 


arm 

c t' = “ ell l + ^ n 2 + ^l 


a 2„ 
a n . 


ox 


+ fee 


p 

a n, 
< 

1 . 2 
dx 


an. 


p 

an. 


3x 


*9 u j-'p o n- 

FT = gri 1 _m 2 +B 2 712” + ga 2 ^2“ 


e = a^+b^L r *, h = a 0 +b 0 H* 


(12.2.7) 


( 12 . 2 . 8 ) 


**1* ^ l* 2 " * S ^ Cl»{p * 

and n^,ng are the perturbed population distributions in 
the linearised model. 

The systems (12.2.5) and (12.2.7) are associated with 
the following flux boundary conditions 


3n an n 

ST (O.tj = (l,t) = 0 

(12.2.9) 

an p an p 

T 5 * (o,t) - j;2 <L,t) - o • 

The boundary conditions (12.2.9 ) represent the fact that the 
population is confined within the region and there is no 
migration across the boundary. 

In the following, we study the linear and nonlinear 
stability of the equilibrium state (n£,N*) under (12.2.3) 
using Liapunov’s Direct Method. 
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12.3 LINEAR STABILITY ANALYSIS 


To examine the linear stability of the endemic equilibrium 
state (II*, H*) in absence of dispersal, consider the following 
positive definite function 


S(n 1 ,n 2 ) = | + pn|) 

where p is a positive constant to be chosen. 
The time derivative of (12.3.1) is 
dE _ n dn l + R n dn 2 

nr - n i ar + jr 


(12.3.1) 


(12.3.2) 


and which on substituting (12.2.7) becomes 

|| = -en^+(f+pg)n 1 n 2 -pim^ . (12.3.3) 


Since eh > fg, by choosing 

P = f/g (12.3.4) 

it may be noted from (12.3.3) that dE/dt < 0 for 
n^ / 0, n 2 / 0| dE/dt = 0 only when n^ = 0 and n 2 = 0. 
Hence the equilibrium state (0,0) [or (B'£,N*) of the 
original system (12.2.1) [] is always asyrnpt otic ally stable. 

The effects of dispersal on the linear stability of the 
equilibrium state (N*,N 2 ) is studied by considering the 
positive definite function 

E(n 1 ,n g ) = | / (n|fpn|)dx (12.3.5) 

where p is chosen as in (12.3.4). 
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'file time derivative of (12.3.5) is 

L an- an 

af = J o ( n i aT + ? n 2 rr> dx 

and which, on substitution ox (12,2.7), gives 

Xl 

= - / jen^-Cf+pg^n^g+phnl J dx 
o 


L 


3 2 n- 


d 2 n c 


+ J UV 1 ! o + $ I) 2“2 ...2 




dx 

2 

3 ua r 


dx 


3 n 


1 — 

+ fa 1 n 1 — ~ + P ga 2 n 2 — 2 ~ -i dx * 

dx dx^ 


(12.3.6 ) 


(12.3.7) 


'The second integral in (12.3.7) is evaluated and making 
use of flux boundary conditions (12.2.9), we get 

||= - J [en|-(f+pg)n 1 n g +phn^dx 


dn 1 2 


dn g 2 


- { L%(ar> + 

dn- dn p 

+ (fa^+pgag) 3 dx * 


(12.3.8) 


Since the first integral is nonpositive, it may be seen from 

(12.3.8) that dE/dt < 0 provided 

iXLjhg > fg(a 1 +a 2 ) 2 (12.3.9) 

and dE/dt = 0 only when n^ = 0 and n g = 0. Hence 

(12.3.9) gives the sufficient condition for asymptotic 
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stability of the equilibrium state in presence of both self 
and cross dispersal co-efficients of the infectives. 

However, when = 6 and ag = 0 i.e. cross dispersal 
co-efficients of tne infectives are sere, the inequality (12.3.9) 
holds automatically showing that the endemic equilibrium state 
which is stable without dispersal remains so with self dispersal. 
Thus, keeping in view of (12.3.9 ), it may be remarked that cross 
dispersal of the infectives have destabilizing effect on the 
equilibrium state of the system and ultimately the entire 
population may be infected. 

12.4 NONLINEAR. STABILITY ANALYSIS 

* .jt. -X- 

The nonlinear stability of the equilibrium state (N^,N g ) 
in absence of dispersal is analysed by considering the positive 
definite function (12.3.1). Using (12.2.5) in (12.3.2), we 
get 

dTf = ~ &i +b i^ T g +n 2 ^ 3 n l &2 +t 2^1 +n l^ -3 n 2 

+ (f+pg^ng (12.4.1) 

where p is given in (12.3.4). 

Since a^ag =fg, it may be seen from (12.4.1) that 
dE/dt < 0 in the feasible region A - {(N^,Ng)| 0 < < c^, 

0 < Eg < Cg > of the phase-plane showing that the equilibrium 
state (N*,N*) is nonlinear ly asymptotically stable in the 

JL <, 

entire region A. 
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lo see the effects of dispersal on the nonlinear 
stability of the eiuiiibrium state (IT*, Eg), consider the 
positive definite function (12.3.5). From (12.3.6; and 
(12.2.5), ve have 

= - J C (a 1 +b 1 E|+b 1 ng;n^+p(ag+bgn*+b g n 1 )n^ 


- (f+pgKn ]dx 


12 

9 0 0 

3 B 1 9 n 2 9 n 2 


+ J 2 + , 2 + fa i n l 2 

o ox 3x dx 


9 2 n 1 9 2 n P a 2 n 1 __ 

+ pga 2 n 2 ^2 b i cc i n 1 —g ^ b 2 a 2 n 2 “H-T 3 dx 


(12.4.3) 


where p is chosen as in (12.3.4). Evaluating the second 
integral and making use of flux boundary conditions (12.2.9), 
we get 

fl = -/Da 1 +tiH2+t 1 n 2 )n2 + p (a2+b2N . +b2 n 1 ) n | 

- (f+pg^ng^dx 


L Sn 1 2 an Q 2 

- { < “hair) + 

an 1 an„ 

+ ff oc ^“*2b^cx ^n^+p (goCg— 2bgtXgng ) ^ ~ ^ dx.( 12.4.4) 

As the first integral is nonpositive, it is noted from 
(12.4.4) that dE/dt < 0 provided 
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4p D^Dg > [a 1 (f-2'b 1 n 1 )+pag(g~21)gn p ) J ~ (12.4.5) 

wnicn represents the region bounded by pair of straight 
lines given by 

2a^b :L n i + 2p«2' i:3 2 n 2 = ^ a i + PS 0 ig- 2 (pD^I)g (12.4.6a) 

2a 1 c 1 n 1 +2pa 2 b 2 n 2 = fa 1 +pga 2 -f2 (plLjDg) 1 / 2 (12.4.6b) 

containing the equilibrium state (0,0) when the inequality 
(12.3.9) holds. If we denote the region bounded by the 

straight lines n^= -ll*,n 2 = -N*, n i =c i -II l» n 2 = c 2“^2 811(1 

(12.4.6a) in the phase-plane by B, then the equilibrium 
state (0,0) [or (N£,Hg) of the original system [] is 
asymptotically stable in the region B provided (12.3.9) and 
(12.4,5) are satisfied. 

Thus it is concluded that if the equilibrium state 
is linearly stable with dispersal, then a subregion B of A 
can be found in the phase-plane where it is nonlinearly stable 
under (12.3.9) and (12,4.5). 

If the parameters are such that 

2(pD 1 D 2 ) 1 / 2 > b 1 a 1 (c 1 +h*)+pb 2 a 2 (c 2 +Hp (12.4.7) 

the straight line (12.4.6a) does not intersect the feasible 
region A and in that case the region B corresponds to A 
showing that the equilibrium state is asymptotically stable 
in the entire feasible region under (12.4.7), As before, 
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when a ± ~ u ccg = 0 i.e. cross dispersal 

co-e if relents of tne inf ec fives are zero, the inequality 
(12.4.7) is always valid snowing that cross dispersal may 
destabilize tne equilibrium state. 

12.5 COHCiUSICl'j 

Using Liapunov's Direct Llethod, effect of dispersal 
on the linear and nonlinear stability of the endemic equilibrium 
state for gonorrhea epidemic model has been investigated. It 
has been noted that the equilibrium state is nonlinearly 
asymptotically stable in the entire feasible region of the 
phase-plane without dispersal. 

The following conclusions may also be drawn with dispersal: 

(i) The endemic equilibrium state is asy up to tic ally stable 
in the feasible region when cross dispersal co-efficients of 
the infect ives are zero. fhis implies that in absence of 
cross dispersal co-efficients, the disease is confined to the 
equilibrium population in the long run as envisaged by Yorhe 

et al • £ 17 33 • 

(ii) In presence of both self and cross dispersal co-eff icients 
of the infectives, if the equilibrium state is linearly stable 
with dispersal, then there exists a subregion in the phase-plane 
where it is nonlinearly stable under (12.3.9) and (12.4.5). 
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(iii) In presence of both self and cross dispersal- 
co-efficients of the infect ives, the equilibrium state 
is asymptotically stable in the entire feasiolo region of 
the phase-plane provided (12.4.7) holds. 

It has also been pointed out from the conditions 
(12.3.9), (12.4.5) and (12.4.7) that the increase in cross 
dispersal co-efficients of the infectives may lead to 
infection of the entire population considered. 
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